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r~| I Abstract Mass loss is a key process in the evolution of massive stars, and must 

I' be understood quantitatively if it is to be successfully included in broader as- 

JL ' trophysical applications such as galactic and cosmic evolution and ionization. 

^ . In this review, we discuss various aspects of radiation driven mass loss, both 

• from the theoretical and the observational side. We focus on developments in 

! the past decade, concentrating on the winds from OB-stars, with some excur- 

sions to the winds from Luminous Blue Variables (including super-Eddington, 
continuum-driven winds), winds from Wolf-Rayet stars, A-supergiants and Cen- 
^ . tral Stars of Planetary Nebulae. After recapitulating the 1-D, stationary standard 

' model of line-driven winds, extensions accounting for rotation and magnetic fields 

00 . are discussed. Stationary wind models are presented that provide theoretical pre- 

■"sj" ■ dictions for the mass-loss rates as a function of spectral type, metallicity, and 

I the proximity to the Eddington limit. The relevance of the so-called bi-stability 

jump is outlined. We summarize diagnostical methods to infer wind properties 
from observations, and compare the results from corresponding campaigns (in- 
OO . eluding the VLT-FLAMES survey of massive stars) with theoretical predictions, 

' featuring the mass loss-metallicity dependence. Subsequently, we concentrate on 

two urgent problems, weak winds and wind-clumping, that have been identified 
from various diagnostics and that challenge our present understanding of radia- 
tion driven winds. We discuss the problems of "measuring" mass-loss rates from 
^ . weak winds and the potential of the NIR Bra-line as a tool to enable a more pre- 
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cise quantification, and comment on physical explanations for mass-loss rates that 
are much lower than predicted by the standard model. Wind-clumping, conven- 
tionally interpreted as the consequence of a strong instability inherent to radiative 
line-driving, has severe implications for the interpretation of observational diag- 
nostics, since derived mass-loss rates are usually overestimated when clumping is 
present but ignored in the analyses. Depending on the specific diagnostics, such 
overestimates can amount to factors of 2 to 10, and we describe ongoing attempts 
to allow for more uniform results. We point out that independent arguments from 
stellar evolution favor a moderate reduction of present-day mass-loss rates. We 
also consider larger scale wind structure, interpreted in terms of co-rotating inter- 
acting regions, and complete this review with a discussion of recent progress on 
the X-ray line emission from massive stars. Such emission is thought to originate 
both from magnetically confined winds and from non-magnetic winds, in the lat- 
ter case related to the line-driven instability and/or clump-clump collisions. We 
highlight as to how far the analysis of such X-ray line emission can give further 
clues regarding an adequate description of wind clumping. 

Keywords hydrodynamics • stars: atmospheres • stars: early-type • stars: mass 
loss • stars: winds, outflows 



1 Introduction 

Within the last decade, our understanding of the physics of massive stars has sig- 
nificantly improved. It has been realized that massive stars are critical agents in 
galactic evolution, during both the present ep och as well as in the early Universe, 
and they are regarded as "cosmic engines" dBresolin et all l2008h . For instance, 
a population of very massive. First Stars (for a recent review, see Ijohnson et all 
12008.') is thought to play a dominant role in the reionization of the Universe and its 
first enrichment with metals. Already with the next generation of telescopes, the 
integrated light from these First Generations of Stars might become observable 
in the near-infrared (NIR), by means of Ly^ (Barton et al 2004) and He II A 1640 
(iKudritzki and The GSMT Science Working Groupll2003h emission from the sur- 
rounding interstell ar med ium, if the c orresponding init ial mass function (IMF) is 
indeed top-heavy dAbel et al 2000, Br omm et aill200ll) . A proper knowledge of 
the mass-loss mechanisms during such early stages is crucial to understand the 
interplay of these First Stars with their environments. Another highlight concerns 
the long gamma-ray bursters (GRBs), which are likely the result of the terminal 
collapse of massive stars at low metallicity. GRBs may also become important 
tracers of the star formation history of the Universe at high redshift. Again, mass 
loss plays a dominant role in the evolution of angular-momentum loss, thus con- 
trolling whether the star will become such a burster. 

There are three factors which have contributed to the aforementioned progress: 

(i) new observational facilities such as ground-based 10-m class telescopes equip- 
ped with multi-object spectrographs, and space-born observatories with spectro- 
scopic capabilities in the ultraviolet (HST, FUSE, GALEX), infrared (SPITZER), 
X-ray (xmm-Newton, Chandra), and in the y-ray (HETE-2, swift) domains; 

(ii) the development of simulations of complex physical processes in the interior 
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and outer envelopes of massive stars, such as rotation and magnetic fields; and (iii) 
the advancement of model atmospheres and radiative transfer techniques. 

Modeling the atmospheres of hot stars is a tremendous challenge due to severe 
departures from Local Thermodynamic Equilibrium (LTE) because of the intense 
radiation, low densities, and the presence of supersonic velocity fields initiated 
by the transfer of momentum from the stellar radiation field to the atmospheric 
plasma. Radiation-driven winds are fundamentally important, in providing energy 
(kinetic and ionizing radiation) and momentum input into the ISM, in creating 
wind-blown bubbles and circumstellar shells, and in triggering star formation. 
They can affect stellar evolution by modifying evolutionary timescales, surface 
abundances, and stellar luminosities. The mass loss does not only influence the 
stellar evolution but also the atmospheric structure, and winds need to be properly 
modeled to derive the correct stellar parameters by means of quantitative spec- 
troscopy. 

Th e last review on radiation-driven winds from hot stars was given bv lKudritzki and Pulsl 
(120001) who concentrated on the derivation and calibration of the so-called wind- 
momentum luminosity relation (WLR) and its potential application as an extra- 
galactic distance indicator. Thereafter, efficient methods have been developed that 
allow for theoretical predictions o f mass-loss rates as a fu nction of stellar param- 
eters in the complete upper HRD dVink et alll2000L l200lh . and these can be used 
to check our standard picture of radiation-driven winds through comparison with 
empirical results. Owing to the enormous improvements in observational capabil- 
ities and diagnostic methods, a number of puzzling phenomena have meanwhile 
emerged, and this has resulted in exciting speculations and theoretical develop- 
ments which are by no means settled. To name two of the most prominent, there 
is the so-called weak-wind problem which indicates that the mass-loss rates from 
late 0-/early B-type dwarfs might be a factor of 10 to 100 lower than theoreti- 
cally expected, challenging our understanding of radiation-driven winds. Second, 
there is the issue of wind-clumping which refers to small-scale density inhomo- 
geneities distributed across the wind. Because clumping "contaminates" almost 
all our mass-loss sensitive diagnostics, it might result in severe down-sizing of 
previously derived mass-loss rates, and consequently impact significantly our un- 
derstanding of stellar and galactic evolution. 

It is the goal of this review to provide an overview of these developments, 
to outline their implications, and to indicate possible directions for future work. 
Because of the breadth of the topic and its complexity, we cannot cover all as- 
pects. Instead, we primarily focus our review on the winds from single, "normal" 
hot massive stars in well-established evolutionary stages such as OB-dwarfs, gi- 
ants, and supergiants, and on some of the most relevant aspects of the winds from 
Luminous Blue Variables (LBVs) and Wolf-Rayet (WR) stars (detailed in a re- 
cent review by Crowther 2007). We also consider the winds from Central Stars 
of Planetary Nebulae (CSPN) which are thought to obey similar physical rules as 
their massive 0-star counterparts. For a review on physical p r ocess es related to 
wind- wind collisions in binary systems, we refer to lDe Becker! (120071) . 

This review is organized in such a way that most chapters can be read inde- 
pendently from each other. Sect. 12.1 [ provides the basic line-driven wind theory in 
terms of a 1-D, stationary model with statistically distributed lines. In Sect. 12.21 
this "standard model" is extended by accounting for the interaction with rotation 
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and magnetic fields, which are not only important when the field strengths are 
large, as it is the case for few "normal'Q OB stars, but also when the wind density 
is quite low and relatively modest fields are present. 

In Sect. |3] we compare and contrast the various stationary wind models pro- 
viding theoretical predictions for the mass-loss rates as a function of spectral type, 
metallicity, and the proximity to the Eddington limit. We focus on predictions ob- 
tained by Monte Carlo methods for stars in various evolutionary phases (OBA, 
LBV, WR) and for various chemical mixtures. These prescriptions are used in 
most modern evolutionary calculations and are generally referred to when obser- 
vational findings are compared to theory. An issue of particular relevance is the 
so-called bi-stability jump. Objects below a specific effective temperature (around 
25 kK, kK = kilo Kelvin) should have systematically larger mass-loss rates and 
lower terminal wind velocities than objects above this jump temperature, due to 
an abrupt change in the ionization equilibrium of iron, the most important element 
with respect to line-driving in the inner wind. 

After a brief discussion of various wind-diagnostics and their pros and cons, 
recent observational results and their comparisons with theoretical predictions are 
summarized in Sect. HI We first demonstrate how the continuing development of 
NLTE model atmospheres has led to severe parameter changes due to the inclu- 
sion of metal-line and wind-blanketing. Most affected is the effective tempera- 
ture scale of Galactic OB-stars, with the new temperatures significantly below the 
older ones. We subsequently present observationally derived wind-parameters for 
OBA-stars, LBVs, and WRs, and we discuss their implications. Special emphasis 
is given to the VLT- FLAMES survey of massive stars in the Galaxy and the Mag- 
ellanic Clouds (lEvans et alllioOSL 120061) which allowed for a derivation of stellar 
and wind parameters of a large sample of objects, enabling a detailed empirical 
study of the mass loss-metallicity dependence. We also cover topics such as the 
winds from BA-supergiants beyond the Magellanic Clouds, the winds from CSPN, 
and we end with an excursion to super-Eddington, continuum-driven winds which 
might be relevant for the giant eruptions of LBVs. 

At first glance. Sect. |4] might give the impression that theory and observations 
largely agree and that only a few items remain to be clarified before hot-star winds 
can be regarded as "understood". Indeed this was a widely held belief at the end 
of the 1990's. In the meantime, however, the issues of wind clumping and weak 
winds (see above) have now emerged as urgent. 

In Sect. [5] we specify the "weak wind problem" and discuss various physical 
processes which might be responsible for the apparent dilemma. We demonstrate 
the challenges for conventional mass-loss diagnostics that become rather insen- 
sitive for very low mass-loss rates. As an alternative diagnostic, we discuss the 
potential of the NIR Bra-line as a tool to enable more precise measurements and 
that might allow for a quantification of the mass-loss "deficit" with respect to the- 
oretical predictions. 

Sect. |6] is devoted to another central topic of current activity in stellar wind 
research, comprising various aspects of wind-structure and time-dependence - 
neglected in the standard model. After a summary of observational findings, we 
provide theoretical arguments which show that line-driven winds are subject to 
a strong instability inherent to the line acceleration itself. This line-driven (or 
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de-shadowing) instability has been proposed to be the origin of wind-clumping, 
triggering the formation of small-scale structure. In Sect. 16.31 we discuss this in 
detail. Clumping has severe implications for the interpretation of observed line- 
profiles, particularly with respect to the derived mass-loss rates. Conventional 
recombination-based p^-processes such as H^, IR, and radio-excess become very 
sensitive to clumping, resulting in overestimated mass-loss rates when clumping 
is present but has been neglected in the analyses. We discuss various diagnostic 
methods to investigate the clumping properties of stellar winds and to derive "true" 
mass-loss rates. Corresponding results disagree with respect to the estimated de- 
gree of mass-loss reduction, ranging from factors in between 2 and 10, and even 
more. Potential problems with these diagnostics have been identified by different 
research groups and are likely related to the porosity of a clumped wind medium 
and the corresponding velocity space, or vorosity. We finish the section with some 
independent insights from stellar evolution, which suggest that the reduction of 
the mass-loss rates is more likely to be moderate. 

Turning to larger-scale structure, Sect. l6.4l covers the discrete absorption com- 
ponents (DACs) and modulation features, as well as corresponding time-dependent 
models of co-rotating interaction zones (CIRs) and co-rotating, azimuthally ex- 
tended structures which might be responsible for these other structural phenom- 
ena. 

Since the first X-ray measurements of massive stars became available, it was 
clear that they are strong X-ray emitters. In Sect. 16.51 we review the status of 
the X-ray wind line emission, as observed by XMM-Newton and CHANDRA. 
This emission has been interpreted to either result from magnetically confined 
winds or to result from shocks embedded in the wind, likely related to the line- 
driven instability or clump-clump collisions. Recent diagnostics has enabled to 
"measure" the onset of this emission, and current analyses utilize the line emission 
as a tool to deduce the degree of wind-clumping, its radial distribution, and the 
shape of the clumps. We report on corresponding results and implications. 

In Sect. |7] we sunmiarize the most important achievements obtained during 
the last decade of wind research and provide some future prospects. 



2 Line driven winds of liot stars - tlieoretical background 

As each ph oton carries m omentum, hv /c, it was speculated as early as in the 
1920s (e.g., iMilnd [l92^ that radiative acceleration by spectral lines might be 
capable of ejecting metal ions from stellar atmospheres. However, it was not un- 
til UV balloon flights in the late 1960s that the theory of radiation line driving 
became the leading theor y for the steady outflows from hot luminous OB stars. 
ILucv and Solomon! ( Il970l) and lCastor et all ( Il975h realized that if the acceleration 
exerted on the metal ions could be shared with the more abundant hydrogen and 
helium species in the plasma, it could result in a significant mass-loss rate M, of 
the order of lO'^M^yr"^, which might affect the evolution of massive stars in a 
substantial manner. 
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2.1 The standard model (1-D, stationary, smooth) 

Momentum transfer by line scattering. In the following, we will concentrate 
on those lines that contribute to most of the total acceleration, i.e., resonance = 
scattering lines from ionic ground states (C,N,0,. . . ) or low-lying meta-stable lev- 
els (iron group elements). The principal concept of momentum transfer by line- 
scattering is that the initial photon originates from a particular direction in the 
stellar photosphere, whereas the subsequent re-emission is more or less isotropic, 
or, to be more precise, at least fore-aft symmetric with respect to the radial di- 
rection. This change in direction angle 6 leads to a radia l transfer o f momentum, 
AP = /;/c(VinCos0in — VoutCosQout) (for details, see, e.g..l Pulslll987f) . and this an- 
gle change (together with the Doppler-effect, see below) is the key of the momen- 
tum transfer and associated radiative line acceleration. 

When integrating this momentum transfer over all scatterings, the losses or 
gains due to re-emission processes cancel because of fore-aft symmetry, resulting 
in a radially directed, total line acceleration, ^J.™ • "^^^^ acceleration acts on the to- 
tal wind plasma, as long as the momentum gained by the metal ions can be shared 
with the more abundant ions of H and/or He which are (usually) less accelerated 
due to far fewer lines close to the stellar flux maximum. 



Coulomb coupling. For denser winds, the momentum transfer from the active 
metal ions to the passive hydrogen and helium ions is easily fulfilled via Coulomb 
collisions, but this may no longer hold for tenuous winds. The one-fluid approxi- 
mation is applicable as long as the active ions are slowed down efficiently through 
electric charge interactions with the passive particles, i.e., if the relevant time scale 
for momentum transfer due to collisions is small in comparison to the time the ac- 
tive ions (accelerated by line acceleration) need to drift apart from the passive 
ones. For a one-component medium, this drift velocity, vj, scales with (with 
density p; for details, see Springmann 1991). Thus, decoupling (vd > Vth, with Vth 
the thermal velocity, here for the dominant passive ion) and subsequent ionic run- 
away can be expected to become most relevant for stars with small m ass-lo s s rates 
and large terminal wind velocities V oc ( Sprinsmann and Pauldrach 1992 , Babell 



ll995Lll996LrKrticka and Kubatll2000LlOwocki and Pulsll2002LlKrti^a et ai2003h . 



Basic assumptions, equation of motion. In order to predict the physical structure 
of a line-driven wind, let us first employ the so-called standard model. This model 
assumes the wind to be stationary, homogeneous (e.g., no shocks or clumps) and 
spherically symmetric, which is justified as long as rotation and magnetic fields 
(covered in Sect. 12.21) can be neglected, since all external forces are purely radial. 

In this case, the total mass flux (or mass-loss rate) through a spherical shell 
with radius r surrounding the star is conserved, as may be seen from the equation 
of continuity 

M = 471 r^p (r)v(/') = constant. (1) 
The corresponding (stationary) momentum equation reads 



dv _ GM 1 dp 
dr r^ p dr 



= 2 T— + ^rad, (2) 
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with inwards directed gravitational acceleration and outwards directed pressure 
{p) term and radiative acceleration. The necessary condition for initiating a stellar 
wind is that the total radiative acceleration, g^^d = ^''"d + ^rad'' exceeds gravity 
beyond a certain point in the outer photospheric layers. With the equation of state, 
p = p, where a is the isothermal speed of sound, and using Eq.[Tl the equation 
of motion is given by 

, fl2 2a^ da^ GM 

vi"^)^:r = ^5 r + ^>-ad- (3) 

dr r ar 

The radiative acceleration. The prime challenge lies in accurately predicting the 
grad term in the equation of motion. For free electrons (the Thomson acceleration), 
this is simply proportional to the Thomson opacity, = Sf.p {s^ proportional to 
cross section) times the flux: 

^rad = ^^^^ ^ ^grav-T, (4) 

with Eddington's F. However in reality line scattering is the more dominant con- 
tributor to the overall radiative acceleration. The reason is that (i) line scattering is 
intrinsically much stronger than electron scattering, due to the resonant nature of 
bound-bound transitions (Gavlev 1995), and (ii) that although photons and matter 
are only allowed to interact at specific frequencies, they can be made to resonate 
over a wide range of stellar wind radii via the Doppler effect (see reviews by 
lAbbottlll984IOwockilll994^ . 

For a single line at frequency v and line optical depth T illuminated by radially 
streaming photons, the line acceleration can be ap proximated in t erms of purely 
local quantities within the so-called Sobolev theory dSobolevI 19601) . Under typical 
conditions, this approximation applies for the supersonic part of the wind, and, to 
a lesser extent, also for the transonic region as long as opacity, source function 
and velocity gradient do not change significantly over a velocity interval Av = Vth, 
corresponding to a geometrically region of width ArKi Vth/(dv/dr) =: the 
so-called Sobolev length. Insofar, the line optical depth in expanding atmospheres 
has an entirely different character than in hydrostatic photospheres: in winds, the 
interaction between photon and absorbing atom is restricted to a region of width 
L^°^, which is usually very narrow, and the line optical depth becomes an (almost) 
local quantity. 

Using this approximation, the line acceleration can be expressed as 

Jine ^vV dv 1 

= ), (5) 

where Ly is the luminosity at line frequency, and the (radial) line optical depth, in 
the Sobolev approximation, is 

T«TS°'' = ^A/(dv/dr), (6) 

with X the frequency integrated line-opacity and X the transition wavelength. Note 
that for optically thin lines, T < 1, the acceleration becomes indendent of the ve- 
locity field and recovers the same functional form as Eq. |4l with a^L replaced by 
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XLy, whereas for optically thick lines, T > 1, the line acceleration depends on the 
spatial velocity gradient dv/dr. This dependence (unique in physics and being a 
consequence of the Doppler effect) is the origin of a number of peculiar features 
related to radiative line driving, as we will see in Sect.|6] 

The final exe rcise is now to sum over all contributing lines, which, follow- 
ing the ideas by ICastor et"ail d 19751 "CAK"), is conventionally done by using 
so-called line-strength distribution functions, which describe the statistical depen- 
dence of the number of lines on frequency position and line-strength, k]^, 

kL = ■ (7) 

P ^eVth 

measures the line opacity in units of the Thomson opacity and, for resonance 
lines and frozen in ionization, remains constant throughout the wind. From in- 
vestigating an ensemble of C III lines (in LTE), CAK modeled the corresponding 
distribution function by a power law, 

dN{v,ki^) = -Nofy{v)ki^"-^dvdkL (8) 

with < a <1, N„a normalization constant and /v(v) the frequency distribution, 
which is rather independent from the distribution in line-strength. More elaborate 
calculations including ions from all participating elements have confirmed this 
approach (see Fig.[Tl left panel). 

Combining the radiative line acceleration (Eq. |5]) and the distribution of lines, 
(Eq. |8]l, the total line acceleration can be calculated by integration, and, following 
CAK, expressed in terms of the Thomson acceleration (Eq. |4| multiplied by a 
factor M, the so-called force-multiplier, 

g=M(0 = ^r«o.(^)« (9) 



(for details see, e. g.. iPuls et al and for an alternative, sometimes advanta- 



geous, formulation lGavlevlll995h r The optical depth parameter, t. 



dv/dr k]^{x=\) 

is the optical depth of a line with an opacity equal to Thomson-scattering, or, alter- 
natively, corresponds to the inverse of the line-strength of a line with unit optical 
depth. Because of the power-law distribution, the line-acceleration is dominated 
by lines distributed around a line-strength ki^{T = 1) = t~^{r) (± 1 to 2 dex above 
and below this value, depending on a). 

The force-multiplier parameter k is proportional to the effective, i.e., flux- 
weighted number of lines stronger than Thomson-scattering, A'effH and can be 
interpreted as the fraction of flux that would be blocked already in the photo- 
sphere if all lines were optically thick, a does not only correspond to the slope of 
the line-strength distribution, but, more generally, quantifies also the ratio of the 

^ For the above power-law distribution, A'etf = N„/{1 — a) x /(J°(LvV/L)/v(v)dv, k = 
(vth/c)r(a) X A'eff, with Vth the thermal speed of a representative ion (chosen as hydrogen by 
CAK) and Gamma-function F(a) fa 1/a for the values considered here. 
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line acceleration from optically thick lines to the total one. In case of deviations 
from a power-law distribution, the radi ative acceleratio n at location r is controlled 
by its local slope a around ki^{x — 1) (iPuls et alll2000h . 

This line force description can easily be extended with an additional depen- 
dence via («eii/W)^ to account for ionization effects dAbbottll 19821) . in the real- 
istic case that the ionization is not frozen in (Wg 1 1 is the electron density in units 
of 10^' cm~^, and W the dilution factor of the radiation field). Under typical O- 
star conditions, 5 is small (of the order of 0.1, see below), but can reach larger 
(positive or even negative) values if either the ionization equilibrium changes sig- 
nificantly throughout the wind (which might happen in dense winds) or the wind 
is mostly driven by hydrogen and helium lines, in case of extremely metal-poor 
stars (see pagell lb. Note that in any case 5 is not a free parameter, but has to be 
determined - in parallel with the other two force-multiplier parameters k and a 
- in a self-consi stent way from the (NLTE-) line opacities and the radiation field 
(see, e.g., Pauldrach et al 1994). 

Finally, the so-called finite cone angle correction factor has to be included to 
accou nt for non-radially streaming photons (jPauldrach et aL198^lFriend and AbbottI 
^Mi- 



Solving the equation of motion. Also for the complete line sample, the radia- 
tive acceleration depends on dv/dr (now through a power of a, Eq. |9|, and it is 
conventionally assumed that this term has a similar meaning as the velocity gra- 
die nt entering t he inertial term on the left hand side (hereafter Ihs) of Eq. [3] (but 
see lLucvlll998h . In this case, the equation of motion becomes non-linear in the 
variable J = r^vdv/dr, and can be approximated by 

y^-GM{\-r)+fkLM'"-y" (11) 

where / is either constant (in the CAK case) or depends on {r^v,Av /dr^M) if 
one includes all subtleties. As the enthalpy term in the energy equation is much 
smaller than the potential and kinetic energy, the above equation has been sim- 
plified by neglecting the gas pressure. For constant /, this equation can be easily 
solved, and also the complete equation provides no real difficulty (see CAK and 
iPauldrach et al|[l98^ . Note that these solutions involve the presence of a critical 
point which controls the mass-loss rate. For models including non-radial photons, 
this critical point is located close to the stellar photosphere and connects a some- 
what shallow (with respect to y), sub-critical flow smoothly with a steeper, super- 
critical one. For further discussion regarding solution topologies and non-Sobolev 

line-fo rces, w e refer the reader to Poe et al ( 1990), Owocki and Puis ( 1999), Feld meier and Shlosmanl 
(I200dl2002h and lMaduraetaii(l2007h . 

As a final result, one obtains the following velocity law and scaling relations 
for Voo and M: 

( \- (]_j>i\P / /3 = 0.5 (CAK case) 
v[r)-v^{v K/r) , 1 j3 = 0.8 (0-stars) . . . 2 (BA-supergiants) ^^^^ 



_ 2GM{l-r) 1 _ r = ij^r CAK case 

''--'--^ ) -^~^esc, |coe«2.25i^ incl. all details 



(13) 



10 



J. Puis, J.S. Vink & F. Najarro 



Moc -r))>->/"' (14) 

with photospheric escape velocity Vesc (already corrected for the Thomson accel- 
eration!) and a' = a — 5. For 0-stars, 5 is usually small (order of 0.1), a' w 0.6 
and k of the order of 0.1. From Eq. [121 we can also derive the behaviour of the 
optical depth parameter, t. It is proportional to (v(r)/v„o)'/^~^, i.e., remains con- 
stant throughout the wind for the CAK case, and increases from outside to inside 

for j3 > 0.5. 

These relations (e.g., Kudritzki et alll 19891 : for first order corrections to Eq. [T3] 



see iKudritzki and Pulsll26ool . their Eq. 15) compare very well with numerical re- 
sults from integrating the complete equation, but are only val id for spatially con- 
stant force multiplier parameters, k, a and 5 . iKudritzkil ( l2002h provides elaborate 
additions to consider depth dependent parameters as well. 

Further assumptions involve the adoption of a core-halo structure that neglects 
continuum formation in the wind (important for dense winds), and the neglect of 
multi-line effects: in the approach as described above, each line is allowed to in- 
teract with unattenuated photospheric continuum radiation, irrespective of its fre- 
quency separation fr om other lines. This neglect leads to overestimated radiation 
forces in most cases ( Friend and Castor 1983 , Abbott and Lucy 1985 , Puis 19871 



iLucy and Ab"bot3ll993U Springmannl ll994lGayley et aliri99llVink et afcOOOh .ln 



Sect. |3] most of these assumptions will be relaxed and more quantitative predic- 
tions will be provided. Note, however, that the basic scaling relations remain rather 
unaffected by these improvements. 

The wind-momentum luminosity relation. Using the above scaling relations, one 
can construct a modified wind momentum rate, Omom 

= Mv„ (/?*//?0)i/2. Given 

that Voo scales with the escape velocity (Eq. [13]) , Z)mom scales with luminosity and 
effective line number alone, 

as long as a' — 2/3 (remember that a' w 0.6 for 0-star). In this case, the effec- 
tive stellar mass M(l — F) cancels in the product Mvoo when applying the scal- 
ing relations, which is rather convenient since masses are difficult to determine 
from observations and, moreover, s eem to disagree when compared with evolu- 
tionary predictions (e.g.. |Pulsll2008l and references therein). Taking the logarithm 
of Eq. ([TSl l and assuming a' — 2/3, we find 

logZJmom w x\og{LlLQ)+D (16) 

(with slope X and offset D, depending on the flux-weighted number of driving 
lines), which is known a s the "W ind momentum luminosity relationship (WLR)" 
([Kudritzki et al 1995, Pu is et al| [l996) and can be considered to be one of the big 
triumphs of the radiation-driven wind theory. This relationship can, e. g., be uti- 



lized to determine extragalactic distances by purely spectroscopic means (IKudritzki et all 
[1999, Kudritzki and Puis 2000, but see Sect. I4.2T I) , and it also played an instru- 
mental role in determining the empirical mass-loss metallicity dependence for O 
stars in the Local Group, which will be discussed in Sect. [4| Equally important. 
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observed and theoretically predicted WLRs can be compared to test the predic- 
tions and validity of the radiation driven wind theory for various spectral types, 
and to identify present shortcomings (e.g., Sect.|5]|. 

Predictions from line-statistics. What determines the slope ? One may now ask 
what controls the slope of the line-strength distribution function (Eq. Ull, and in 
particular, why is a w 2/3? Since the line-strength k]^ is proportional to the prod- 
uct of absorber density per unit mass (depending on abundance, ionization frac- 
tion and excitation) and oscillator-strength, the answer to this question is rather 
difficult, in particular since a large fraction of the line acceleration, at least in the 
lower wind, is due to lines from iron group elements (see below) with a complex 
electronic configuration. In the following, we will highlight certain aspects (for 
details, see Puis et al 2000). 

Most interestingly, it is rather easy to show that a would be almost exactly 
2/3 if all lines would originate from hydrogenic ions and would be resonance 
lines, i.e., excitation effects could be neglected. Then the line-strength distribution 
function depends solely on the distribution of the oscillator strengths /, and from 
the well known Kramers-formula for hydrogenic ions with transitions 1 ^ « (with 
principal quantum number n) and corresponding selection rules, one finds that 
/(l,n) w const/n^ and d/V/d/ /"^Z^. Since, under the discussed conditions, 
this slope corresponds to a — 2, we immediately obtain a = 2/3. 

Similarly important is the result that for trace ions one stage below the major 
ionization stage the equality a + 5 w 1 is valid throughout the wind. Thus, for hot, 
metal-poor stars driven predominantly by Hi and He II, the 5-term can become 
much more influential on the wind solution than in "conventional" winds driven 
by ions o f the domi r iant st age, where 5 w 0.1. This analytical result has been 
proven bv lKudritzkil(l2002h by means of consistent wind calculations, and is of 
relevance in the context of radiation driven winds from the First Stars. Note that in 
these extreme cases the scaling relations for M and Voo as provided above become 
invalid, since they are based on small values for 8. 

For non-hydrogenic ions (with many and low lying levels), the situation is 
more complicated. Per ion, the typical line-strength distribution consists of two 
parts, a steeper one dominated by excitation effects and a flatter one following 
the oscillator strength distribution, which, after summing up over all participat- 
ing ions, becomes almost irrelevant. Most important for the shape of the total 
distribution is the difference in line-statistics between iron group and light ions 
as well as their different (mean) abundances. Since the former group comprises 
a large number of meta-stable levels, the line number from iron group elements 
is much higher, especially at intermediate and weak line-strengths. Additionally, 
this number increases significantly with decreasing temperature (more lines from 
lower ionization stages), which becomes important with respect to the so-called 
bi-stability mechanism (see below). 

For solar abundances, iron group elements dominate the distribution at low and 
intermediate values of line-strengths, i.e., they dominate the acceleration of the 
inner wind (large t) and thus determine the mass-loss rate. Lines from light ions 
(including hydrogen under A-star conditions), on the other hand, populate the high 

end and accelerate the outer wind (low t), i.e., determine v«,. Typically, this part 
of the distribution is steeper than the rest (e.g., Fig.[T]), due to excitation effects. 
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log k^ metallioity 7 (log) 

Fig. 1 Left: Logarithmic plot of the line-strength distribution function for an O-type wind of 
40 kK (solar abundance), and corresponding power-law fit. 

Right: Predictions from line-statistics. Dependence of mass loss (via , cf. Eq.[T7} on metal- 
licity, for reff= 4 kK (dotted) an d reff= 10 kK (dashed). The slopes are 0.56 and 1.35, respec- 
tively. Data from lPuls et aiH2000l Table 3). 



and only a small number of lines (lKrtickall2006 l: "a few dozens") is responsible 
for this acceleration! Insofar, the rather large observed variance in Voo/vgsc even 
for stars of similar spectral type (Sect. l4~2!T] i is not surprisi ng, since subtle effects 
(small variations in temperature, density, composition - cf. Pauldrach et al|[l990l 
their Fig. 8 - and also X-rays, see Sect. 16.1b can have a significant effect on the 
population and line-strengths of light ions. 

Note finally that the above results, which are based on the investigation of line- 
strength distribution functi o ns alone, have be en confirme d by direct win d calcula- 
tions, see lPauldrachl (Il987l) . lvink et all d 199911200 IH and lKrtickal (l2006h . 



Predictions from line-statistics. Dependence on temperature, wind-density and 
metallicity. Summarizing the above paragraphs, deviations from a perfect power 
law are to be expected, mainly due to the different distributions of "light" and 
"heavy" ions, their specific dependence on temperature and the fact that the dis- 
tribution at the high end is controlled by excitation effects for the predominant 
ion species, inducing a steeper slope. Translated to the effective value of a (i.e., 
evaluated at = f~^) which controls the line-force (and the slope of the WLR), 
this quantity should decrease with decreasing Teg, decreasing wind density and de- 
creasing metallicity Z. In all these cases, the WLR is expected to become steeper 
than in the "standard" situation encountered for 0-stars where a is of the order 
of 0.6. . . 0.7. For A-supergiants, e.g., line-statistics predict a w 0.45 such that in 
addition to a steeper WLR also the scatter introduced by a non-vanishing mass 
term should become significant. 

Besides the variation of a with Z, the major consequence of a different metal- 
licity is a variation of k, i.e., a different effective line-number. Roughly speaking, 
the higher the metallicity, the more effective lines are to be expected, and vice 
versa. Since varies proportionally to Z, as long as the ionization equilibrium re- 
mains rather unaffected from the different metallicity, the assumption of a power- 
law line-strength distribution function (Eq.[8| immediately implies No (Z) oc Z^~", 
with similar scalings for k and A^eff - By means of Eq. [141 the mass-loss rate should 
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then scale according to 

M{Z)ocZ^. (17) 

From Fig.[ll right panel, we see that the predicted dependence of k{Z) oc z'~" is 
actually present, and the derived slopes (5 = (1 — a) /(a — 5)) are consistent with 
prototypical values, namely (a, a')=(0.68, 0.58) for the 0-type and (a, a')=(0.48, 
0.38) for the A-type models, respectively (lower a because of lower Teff: more Fe 
lines from lower stages). Moreover, these numbers are also consistent with results 
derived from the analysis of a large sample of 0-stars in the Galaxy and the Mag- 
ellanic Clouds (Sect. l4~2TI) . Detailed measurements of the metallicity dependence 
of A-supergiant winds are still missing though (but see page lSOl) . 

77ie bi-stability mechanism. Most of the previous paragraphs may have given the 
reader the impression that - save for their intrinsic luminosity dependence - the 
mass-loss properties from the winds of massive stars are more or less identical, but 
nothing is further from the truth: the wind properties of massive stars are intrigu- 
ingly diverse ! This diversity is, of course, related to the variation of effective line 
number and a' as a function of spectral type and metallicity and to the different 
efficiencies of multi-line effects as a function of wind density. 

Of particular interest is the so-called bi-stability mechanism, which was first 
encountered in the model calculations of t he famous Luminous Blue Variable 
(LBV) P Cygni b y iPauldrach and Pulsl (ll99Q) . At a critical effective temperature 
around 21 000 K dLamers et al 19951). the over all properties of Voo and M are ex- 



pected to change drastically. IVink et all d 19991) showed that when a massive star 
evolves from hot temperatures to lower ones, the dominant wind driving element 
of iron (Fe) recombines from FelV to Fe III, which is predicted to result in an 

increase of M by a factor 5 and a drop of Vo o by a facto r 2. 

As was pointed out already by Pauldrac h and Pulsl (il990h . such a recombina- 
tion can be triggered by the optical depth of the hydrogen Lyman continuum]! 
As long as this is optically thin (e.g., for hotter temperatures), the EUV contin- 
uum is rather strong, such that most iron resides in Fe IV, which has fewer lines 
than Fe III, resulting in a lower density, faster wind. When the Lyman continuum 
becomes optically thick, there is almost no EUV flux left, and iron recombines 
to Fe III with more lines, inducing a dense, low velocity wind. Since the optical 
depth of the Lyman continuum is closely coupled to wind density via recombina- 
tion effects, 

oc p2 oc (^)' oc (M(i -r)) « (M(i -^))-'•■-^ (18) 

the abrupt (see below) change in wind conditions does not only depend on Tfff, 
but also, and severely, on effective mass (e.g.. Figs. 4 and 5 in jPauldrach and Pulsl 
[1990). As we will see in Sect. l2.2Tl this can play an important role in understand- 
ing the B[e]-star phenomenon. 

^ At least for dense winds such as from P Cygni considered in the original explanation. For 
most "normal" massive stars as studied by Vi nk et al (J999, 2000), the Lyman continuum re- 
mains optically thin between photosphere and critical point, and the abrupt change of the Fe 
ionization equilibrium is rather due to conventional temperature effects and the runaway in 
density as described below. For these objects we consider TLy as an alias for the Fe lll/Fe iv 
ionization equilibrium. 
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The reason that the wind should react to the different accelerating conditions 
by an actual "jump" in M and Voo is that the physics of the line driving implies that 
when the density increases due to an increasing M, the inner wind is anticipated to 
recombine even more, leading to a run-away mass-loss increase until all iron has 
accumulated in Fe III. The line-driving in the outer wind is expected to react to 
an increased inner density by lowering its terminal wind velocity (remember that 
•?rad °^ P Eq.|9]l. In Sect. l4.2.Tl we will discuss the empirical constraints on the 
bi-stability mechanism. 

Also here, the spectral-type dependence of the mass-loss rate and terminal 
velocity occurs in conjunction with changes in A'gff and a', which could subse- 
quently be reflected in the slope and offset of the WLR. Let us emphasize that the 
spectral-type de pendence of the WLR represents a most act ive aspect of hot-star 
wind research (feenaglia et ail l200llMarkova and Pulsll2008h . 



2.2 The impact of rotation and magnetic fields 

All previous results have been derived under the assumption of spherical sym- 
metry, with gravity and radiative acceleration as the only external forces. In the 
following sections, we will summarize alterations due to rotation and magnetic 
fields, if present. 

2.2.1 Rotation and stellar winds 

(Almost) all massive stars start their evolution as rapid rotators, and remain rapidly 
rotating during the largest part of their life-time, at least with respect to the most 
decisive quantity, namel y the so-called "critical ra tio" between surface rotation 
Vrot and critical velocity dMeynet and Maedel2000l) . Vrot itself decreases (almost) 
continuously during the main sequence and evolved phases, where the rate of this 
decrease is larger for the more massive stars due to the larger mass loss and ac- 
companying larger removal of angular momentum. For supergiants, the reduction 
is due to an increasing stellar radius as well. 

The role of rota tion as a fundamen tal stellar parameter of massive stars has 
been highlighted bv lLanger et all (Il997l) . and effects on stellar structure and evolu- 
tion were detailed by the Geneva group (A. Maeder, G. Meynet and co-workers) in 

a series of publ ications. For corresponding revie ws, we refer to [Maeder and Mevne^ 

(I2000bll2003bl) (see also material presented in lWolf et all 19991 and lMaeder and Eeiiensl 
120041) . 

Rotation affects the physics and diagnostics of line-driven winds from mas- 
sive stars at least in three ways: (i) the dynamics is modified, (ii) the diagnostics, 
particularly for M, needs to consider rotation due to deviations from spherical 
symmetry, and (iii) line-profile variability is induced due to an interplay between 
disturbances in the photosphere (stellar spots, non-radial pulsations), rotation and 
wind-dynamics (co-rotating interaction regions, CIRs, see Sect. 16.41 ). 

1-D models, scaling relations. First attempts to include rotation into the stellar 
wind dynamics were based on a 1-D solution in the equatorial plane, and have 
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been performed by iFriend and Abbott (Il986h and lPauldrach et all(ll986h . with ap- 
plications regarding outflowing disk s of Be-stars b y de Arauid (Il995h. and iri ves- 



tigations of the sol ution topology bv lCure and Riall (l2004l see also ICurell2004l and 
iMadura et afcoOTt) . 

Under the assumption of purely radial line-forces, the equation of motion is 
governed by central forces, i.e., the angular momentum is conserved, L oc rv^ = 
const, and the azimuthal velocity in the equatorial plane is given by v,|)(r) = 
Vy:oi{Rt)Rt./r, with polar velocity ve(r, 9 = 90) = and co-latitude 9. Thus, the 
equation of motion remains almost identical to Eq. [3] or the simplified version 
Eq. [m except that the effective gravity (here defined as the difference between 
gravitational and Thomson acceleration) has to be corrected for the centrifugal 
acceleration v^^/r, 

r-^ ^ r Vait 

with Vciit the 1-D critical velocitjQ, Vcrit = Vesc/\/2- In the following, we will 
mostly consider models not too close to critical rotati on {D. ^ 0.75 ) . For a discus- 
sion of more extreme situations in one dimension, see ICurel (12004 . ICure and Riail 
(12004) and lMadura et all ( l2007l) . Since in the former case the critical point (page|9]l 
of the non-rotatin g solution remains almost unaffected by the centrifugal term due 
to r„i r w R* (^e.g.,'Pauld rach et allll98dlOwocki and ud-Doulall200llMadura et all 



the effective mass can be replaced by M(l — F) (1 — Q ), and M scales, to 



I2007h . 

a very good approximation, with 

M(i2) wM(0)(l-X22)i-V«'. (20) 

Similarly, since Voo scales with Vesc, a rather good approximation for the depen- 
dence of Voa on fl is 

v<„(r2)wv„e(0)(l- 122)1/2. ^21) 

For independent arguments leading to this equation, see lGavlevI (l2000l) . Compar- 
isons of exact vs. approximate solutions for M and Voo as a function of Q have 
been given, e.g., by .Puis et al (.1999,) . 



Wind c ompressed disks and zones. In a seminal paper, 'Biork man and Cassinellil 
(Il993l) included the latitude dependence of the centrifugal forces into the previous 
1-D description in order to generalize the approach. Since (at that time) it was 
still undebated to assume line-forces which are purely radially directed, under 
this assumption the angular momentum is conserved for each particle starting at 
a given co-latitude < 9q < TZ, and its trajectory is restricted to an orbital plane 
tilted by 9o with respect to the equatorial one. Thus, a 1-D solution is still possible 
per orbital plane, if one accounts for the centrifugal acceleration as a function of 
co-latitude, 

R 2 

^cent = ( Vrot ) sin Qq ) ^ — ^ , (22) 



for the radiatively reduced gravity as defined above and neglecting gravity darkening. 
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and the scaling relations for the equatorial flow (Eqs. |20l f2]\ can be generalized 
with respect to do, 

M(r2,0o)°= (1-^22 sin^Qo)'"'/"', v„{n,do) {I - sin^ do)^^^ . (23) 

In this description, M increases and Voo decreases towards the equator, respectively. 
Since the particle is restricted to its orbital plane, the corresponding azimuthal 
angle, 'i>'{r), increases for increasing r and can be calculated from integrating 
d4''(r) /dr = (r, 9q) / {rVr{r, do)), using a j3-velocity law (Eq. ll2l) in combination 
with Voo(0o) for the radial velocity component v^. 

<P'(r) may become large if either Vrot is large, Voo(0o) is small, the velocity 
field is flat (j3 large) or the particles start from equatorial regions. In any case, 
however, the particles move towards the equator and, if 'P'{r) became n/2, would 
collide in the equatorial plane with particles from the opposite hemisphere which 
have started at similar (absolute) latitudes. Since the equatorwards directed motion 
corresponds to a non-vanishing polar velocity component in the stellar frame, and 
this component is small, but supersonic for large <P\ a shock develops, and a 
"wind-compressed disk (WCD)" will be created, confined by the ram pressure of 
the wind. The WCD model h as been confirmed i n all aspects by numerical models 
using identical assumptions dOwocki et alll 19941 and Fig.|2l left panel). 

Even if the deflection in <P' is too small to allow for the creation of a disk, 
i.e., the equatorial plane cannot be reached by certain trajectories (because of too 
low rotation or too steep a velocity field) the above scaling relations (Eq. |23]) 
remain valid, and a non-spherical, oblate wind is predicted, in this case denoted 
as a " wind-compressed zone (WCZ)" model. For such models, Petrenz and Pul^ 
(Il996) have calculated profiles as a function of Vrot and v sin /, and concluded 
that standard 1-D diagnostics might lead to severe overestimates of M (50 to 70%) 
in those cases were v sin i is small and the star is observed pole-on. 



Non-radial line forces: the inhibition effect. After publication, the WCD model 
was regarded as a major step towards an explanation of disk phenomena in hot 
massive stars (e.g., Be-stars) until it be came c lear that it was based on an insuf- 
ficient assumption ( Owocki et al 1 996 . 1 1 998ah . Since, in a more-D situation, the 
optical depth depends on the directional derivative of the local velocity field (the 
radial gradient dv/dr in Eq. |6]has to be replaced by n • V(n • v)), and the radia- 
tive acceleration of strong lines is proportional to the integral over solid angle of 
this derivative times direction n, the line force does not only depend on the radial 
velocity gradient, but also on other velocity gradients present in the flow. Most 
importantly, the dependence of v,- on polar angle (increasing towards the pole, see 
Eq. |23] | leads to a small but substantial polar acceleration, which is directed to- 
wards the pole as long as dvr/dd < OEI i.e., as long as the equatorial outflow is 
slower than the polar one, which is (almost) inevitable for a rotating line-driven 
wind (due to the decrease of the effective, centrifugally corrected gravity towards 
the equator and Voo °= v^sc)- This polar component of the line force, 

g]!"^ ~ 4 dvr/dr \ r ) 



5 more precisely, if (r|:^ -I- < 



Mass loss from hot massive stars 



17 
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- -■ o - I ■ . - 1 ^ 

log density/ {g/'cm**3,) 



log dcnsity/(g/'cm+*3,) 

Fig. 2 Density and polar velocity of radiation driven wind models for a rapidly rotating B-dwarf, 
using different assumptions. The rotation axis is directed along the ordinate. 
Left: radial line force wind compressed disk. Middle: non-radial line force, uniform bright- 
ness temperature. Right: non-radial line force + gravity darkening. 

Model parameters (polar values): T^ff^ 20 kK, logg= 4.11, 4 M = 7.5 Mr^. Vi-ot= 
350 k ms~', V e sc(pol e) = 843 kms^', Vesc(equator) = 585 kms^', Vci-it= 487 km s^' (see Ea.l26t. 
From iPetrenzin99^ . 



is sufficient to stop the small, equatorwards directed vq component of the WCD 
model a nd to induce a polewards directed velocity field, as shown by Owocki etall 
( Il998ah . This is the so-called "inhibition" effect (lOwocki et allll996), which can 
be see n in all numerical models including non-radial line forces (e.g.. lOwocki et ad 
119961 and Fig. El middle panel). 



Gravity darkening: The QF limit and prolate winds. Because of the centrifugal 
forces, not only the wind is disturbed, but also the stellar surface distorted. The 

corresponding shape can be calculated from a Roche potential (see Collins 1963', 

Collins and Harrington 1966 , and for applications iCranmer and OwockL1995, . ,Petrenz and Pulsl 
1996h . [n this case, and with angular velocity O), 



pmax 

^2 = ^ = ^^, (25) 

Oicdt Vcrit ^eq('") 



^ First o bservations of such a surface distortion have been obtained by 
iDomiciano^ e Souza et al (2003) for the rapidly rotating Be-star a Eri with the VLTl, re- 
sulting in a ratio between equatorial and polar radius of 1.56±0.05. The shape of the observed 
distortion seems to indicate that the conventional assumption of a uniform rotation might not 
apply for this star. 
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Q becomes a function of equatorial radius /?eq (which itself is a function of Q). 
Q = I (i.e., "total gravity" = at the equator, see below) occurs when the equa- 
torial radius has reached its maximum extent, Z?™" =1.5 Rp, with polar radius 

The critical velocity, Vcrit, requires special attention, due to the problem of 
gravity darkening. Af ter some controversial discussions arising from the sugges- 
tionof an Q. limit by iLanged (I1997L Il998h which h as been criticized by Glatzej 
([1998) (because of disregarding gravity darkening), ' Maeder and Mevnetl (i2000ah 
performed a detailed study on the issue. According to this study, for not too large 
r the critical speed has to be calculated from the actual mass (independent from 
F) and from the maximum equatorial radius, consistent with the arguments pro- 
vided by Glatzel, 



vcnu = (^)^ = (^)^ r< 0.639. (26) 

Note that this value is lower by a factor of ^J\2/'i) than in the corresponding 1-D 
case, due to the difference of polar and equatorial radius. 

For such a distorted surface with latitude dependent effective gravity (in the 
following defined as gravity corrected for centrifugal terms), one has also to con- 
sider gravity darkening as a consequence of the von Zeipel-theorem dvon Zeipell 
119241) which states that the surface flux of a rotating star is proportional to the 
local effective gravity^Thus, the effective temperature depends on the latitude (in 
the example of Fig.|2] T^^ (pole) = 21.3 kK and Tgff (equator) = 16.9 kK), which 
has to be accounted for in the calculation of the occupation numbers (through the 
radiation field). 

The principle pr esence of gravity darkening has been recently confirmed by 
lMonnieretailJ2007[) . who combined the NIR radiation of the rapidly rotating A7V 
star Altair collected from four telescopes of the CHARA interferometric array 
to obtain an intensity image of the surface of this star. Interestingly, also their 
measurements (remember the problem for a Eri) showed a somewhat different 
temperature distribution compared to the standard approach assuming a uniformly 
rotating star and the von Zeipel theorem. 

In any case, since the flux is modified by gravity darkening , also the F-term 
needs to be modified (lMaedeilll999LlMaeder and Meynedl2000ah . Within the stan- 
dard approach and to a good approximation, it turns out that 



in 



(27) 



if one expresses the "total" gravitational acceleration (centrifugally corrected grav- 
ity, geg, minus radiative acceleration) as gtot = ^eff(l — Tq), to preserve the for- 
mulations obtained so far. The function f{fi) varies from 1 at low v^ot to 0.813 w 

' Though the original version of this theorem was restricted to the case of rotational laws that 
can be derived from a potential (e.g., cylindrical rotation), more recent investigations by Maede^ 
j 19991) have proven that this theorem is also applicable (within a 10% error) to the case of 
"shellular" rotation (co = (o(r) for not too fast rotation) w hich has bee n proposed as the relevant 
rotational law at the surface levels of non-magnetic stars l lZahn|[T993) . 
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1^(2/3) at critical rotation, and can be neglected in most cases. For large F then, 
the critical speed needs t o be redefined, and can be approximated (for exact ex- 
pressions and details, see lMaeder and Mevnetll2000ah by 



^2' 



vcrit,2«vcrita(^) (l-r)2 r> 0.639. (28) 



which results from the cond ition gtot = when fh = 1, the so-called ^2F limit 
dMaeder and Meyne3l2000ah : For large F, the critical velocity is significantly re- 
duced by the proximity to t he Eddington-limit. Thus, there is indeed a lower crit- 
ical speed (as suggested by iLangeJl 1 997l [r998l) . at least when F > 0.64, but it is 
important to emphasize that the i2F limit occurs over all latitudes, not just near 
the equator, as assumed in the original idea of the D, limit. 

Note that the redefined I}2-term needs to be considered in fast rotating models 
close to the Eddington-limit, a problem which has been neglected in all hydrody- 
namical wind calculations performed so far. In the present context, however, the 
most important effect of gravity darkening concerns the line force, due to the mod- 
ified illumination of the lines. In a rotating wind illuminated by a uniformly bright 
disk, the mass flux increases from pole to equator according to Eq. 123 Fl The in- 
clusion of gravity darkening into the luminosity dependence of the mass-loss rate 
(Eq.IHl, 

thus leads to as\ additional dependence on (1 — I2^sin^ 0)]'/" , such that in total 

M(0) oc V«'Meff(0)'"'/"' oc (1 - 12^ sin^ 0) + ' (gravity darkening included) 

(29) 

the situation becomes reversed and the mass flux is predicted to increase from 
equator to pole, in paralle l with the correspon di ng terminal velocity, inducing a 
prolate wind structure (cf. lOwocki eta l' 1998a'9,|PeiuEessy et al 2000 and Fig.|2l 
right panel ). The present-day stellar wind of r\ Car has indeed been found to be 
elongated dSmith et alll2003al Ivan Boekel et alll2b03h with a position angle equal 
to that of its ejecta. 

Though the predicted angular dependence of the mass-loss rate is not affected 
by the modified r^-term (except for its influence on Q. via Vcrit), this modification 
plays an important role when comparing the total, polar-angle integrated mass- 
loss rate with mass-loss rates from non-r otating models. To a simila r degree of 
precision as valid for Eq. |27] one obtains dMaeder and Mevneil2000ah 

^(Vrot) ^ (1-n'/"'"' .30. 

M(0) ('i_4_vW^_j.xl/a'-l' 
^ ^crit, 1 ~ 

Thus, no dramatic influence of rotation on the total mass loss is to be expected in 
most cases. 

^ Note tiiat tlie inhibition effect does not change this basic result, since it affects "only" the 
small polar velocity component and thus the disk formation. 

' The formal divergence at the CIF limit is artificial, since in this case the material can no 
longer be lifted from the gravitational potential. 
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Fig. 3 Predicted WLR for rapidly rotating B-type stars (vi-ot= 0.85 v^nt) with Tf.ff= 20 kK and dif- 
ferent radii. Compared are the modified wind-momenta D,^om calculated from total (diamonds), 
equatorial (squares) and polar (triangles) mass-loss rates and average terminal velocities result- 
ing fr om 2-D, NLTE simulati ons, and from corresponding 1-D models with Vi-ot = (crosses). 
From lPetrenz and Pulsl ( 120001) . 



2-D NLTE models and the WLR of rapidly rotating stars. Strictly speaking, the 
scaling relations derived in the previous sections are only valid if the line-statistics 
remains unaffected from the modified wind structure. In order to investigate the in- 
fluence of density, velocity an d radiation field on the oc cupation numbers and line- 
forces as a function of (r, 6). IPetrenz and Pulsl (l2000l) calculated self-consistent, 
2-D, NLTE models for rapidly rotating B-star wind£3. to investigate a spectral 
regime where the ionization structure is most sensitive to local conditions and 
variations of the radiation field. For all models, a prolate wind structure has been 
found, in accordance with the principal predictions from above. 

As obvious from Fig. [51 the total wind momenta are hardly affected by rota- 
tion, since the increase in M, Eq.[30l is compensated by a decrease in the average 
Voo. Of course, the wind-momenta differ significantly when observed either pole 
or equator on. This problem can be mitigated when M diagnostics is used which 
scans mostly the lower wind region, e.g., H^, since in these regions the density 
contrast between polar and equatorial regions is lower. Consequently, WLRs de- 
rived from 1-D analyses should remain rather unaffected from deviations from 
spherical symmetry, when samples of significant size are used, and objects with 
low V sin / (likely to be observed pole-on) are avoided. The scatter in the WLR 
increases due to rotational effects, of course. 



B[e]-supergiants an d the bi-stability. B[e] supergiants (for an improved classi- 
fication scheme, see iLamers etallll998h show a hybrid spectrum (narrow B aimer 
lines and low ionization metal lines in emission, but also UV and Balmer P Cygni 
lines and a strong IR excess), which has been explained by an outflowing disk 



computational feasible at that time because they applied a Sobolev line transfer and assumed 
an optically thin continuum. 
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Fig. 4 Onset of bi-stability (co-latitude 6 where TLy = 1 is reached) as a function of TLy(polar) 
and three rotational speeds, Q = 0.5,0.7,0.9 (solid, dotted, dashed), f or a B -supergiant with 
7^2=20 kK. Optical depth of Lyman continuum calculated according to lLame rs and PauldracQ 
( fT99l ). Left: original version. Right: mass-loss rate calculated including gravity darkening. The 
formation of a disk requires much larger polar mass loss and rotation than in the original version, 
see text. 



with high density, low veloc ity and low ionizatiqn and a fast polar wind with low 
density and high ionization dZickgraf et alll 19861 Il989h . 

Th e origin o f this phe nomenon has been suggested by iLamers and PauldrachI 
(ll99lL see also lLamers" et al 1999) to rely on a combination of rotation and bi- 
stability. Since, for a rotating wind, M was thought to increase and v„ to decrease 
towards the equator (the impact of gravity darkening on M was yet unknown), the 
optical depth of the Lyman continuum 

TLy=/[(^)',7;,d] (31) 

becomes a strongly increasing function of 6, particularly if including the von 
Zeipel theorem when calculating Ty^d{(^) which establishes the ionization equilib- 
rium03] Lamers and Pauldrach (1991) were able to parameterize the optical depth 
as a function of latitude-dependent gravity, TLy °= (^poie/^(0))^^- Because of the 
dependence of the bi-stability jump on this optical depth, the latitude where TLy 
reaches unity separates two regions, a polar one with low density and high velocity 
(predominantly driven by Fe IV) and an equatorial one with high density and low 
velocity, driven by Fe III (see Sect. |2.1| ). The contrast between both regions results 
in an equatorial outflowing disk. From Fig. ID left panel, one can see that even a 
thin polar wind with TLy = 0.01 (corresponding to M= 0.35 •lO~^M0yr~') can 
switch almost abruptly to TLy > 1 at w 40° for Q = 0.9 in the outlined scenario. 

After the impact of gravity darkening on the mass-loss rate had been realized, 
lOwocki et all(ll998a[) pointed out that the creation of B[e] star disks is much more 
difficult than previously thought since M decreases equatorwards (Eq.|29ll, in con- 
trast to the original approach. Thus, the increase of TLy towards equatorial regions 
is much weaker, TLy °= (gpoie/giS))^'^ , and even for highest rotational speeds a 
disk formation becomes impossible as long as the polar wind is not considerably 
dense (Fig.|4l right panel). 

" In a consistent calculation, the von Zeipel theorem has be employed both for the NLTE- 
aspect (occupation numbers) and for the flux-dependent line-acceleration. 
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Based on their previous investigations. lCure et all (l2005h suggest an alternative 
possibility to explain the B [e] phenomenon. They propose that for near-critical ro- 
tation the wind can switch to a slower, shallow-acceleration solution which can 
lead to a density enhancement and, again in combination with the bi-stability 
mechanism, to the formation of a dense equatori al disk. From analyti c consider- 
ations and (1-D) time-dependent hydrodynamics. iMadura et al (l2007h confirmed 
most of these results, in particular the slow acceleration for 12 > 0.85 and a cor- 
responding increase of density. However, they argue that an inclusion of gravity 
darkening and 2-D flow effects still renders the formation of dense disks as un- 
likely. 

Let us also note thatlZsargo et all (l2008b l) have identified an additional prob- 
lem. From 2-D, NLTE calculations by means of their axial-symmetric, NLTE 
code ASTROROTH they concluded that hydrogen recombination in the disks of 
B[e]-supergiants is almost impossible, due to strong ionization from excited lev- 
els. Thus, the chance that TLy becomes much larger than unity and enables line 
driving by Fe III is rather low, even if the problem of gravity darkening was not 
present. 

On the other hand, IPelupessy et all (120001) solved the momentum equation for 
sectorial line-driven winds from rotating stars including both the stellar oblateness 
and gravity darkening, as well as force-multiplier parameters specifically com- 
puted for the bi-stability regime. In their computations for stellar parameters far 
from the bi-stability limit, they found the mass-loss rate and the terminal wind 
velocity to increase from the equator to the pole (thus confirming the reliability 
of their approach), however when the counteracting effect of wind bistability was 
included, they were able to produce a density contrast between the equator and 
the pole of a factor 10 for a relatively modest rotational velocity of 12 = 0.6. One 
of the attractions of the scenario remains that it may naturally account for the 
expectation to find outflowing disks only around B-type stars. 

From the above discussion it should be evident that more work (in 2-D) is 
needed to prove whether or not the bistability mechanism is the root cause for the 
existence of disks around B[e]-supergiants. 

2.2.2 Winds and magnetic fields 

In contrast to the Sun and other cool stars where magnetic fields are thought to 
be generated through a dynamo mechanism related to convective motions in the 
H/He recombination zones, hot stars have been considered to be free of strong, 
dynamo-created surface magnetic fields until recently, because of the absence 
of strong outer convection zones (H remains completely ionized even through- 
out the photosphere). There are, however, several other possibilities to generate 
surface magnetic fields, (i) either related to a thin convection zone due to re- 
combination of He III; (ii) due to the strongly convective cores of massive stars 
which might g ive rise to d ynam o-generated flux tub es diffusin g to the surface 
dCassinelli and MacGregoii l2000i but see iMacDonald and MullanI 12004 (iii) by 
a dynamo that is operating i n shea r-unstable gas in the radiative stellar enve- 
lope (MacDonald and Mullanl l2004l) : (iv) formed already in an early convective 
phase during stellar formation; (v) or even through compression of int e rstella r 
magnetic flux during the initial collapse. It has been shown by IMossI (l200lh . 
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iBraithwaite and Spruij (l2004l) and iBraithwaite and Nordlundl (l2006h that a dy- 
namically stable configuration of such fossil fields on long time-scales is actually 
possible. 

Another promisin g mechanism is the Tayler-S pruit dy namo (ISpruilll999L 2002 
and g eneralizations bvM aeder and Meynetl2004| ; see also lMullan and MacDonalc 
l2005h which does not need a convection zone but is based on a strong instability 
for generating magnetic fields in radiative layers of differentially rotating stars. 
This dynamo has been incorporated into recent stellar structure/evolution calcula- 
tions ( Maeder and Mevnet 2003a, .2004, .2005„ Heger et al 2005, Brott et ali2008i) 
to understand the internal angular mo mentum transport (though the induced mix- 
ing processes might be too strong, e.g.. lBrott et alll2008h . Note, however, that this 
dynamo is able to create magnetic fields (in the azimuthal direction, of the order of 
a few 10^ Gaus s (G)) only outside the stel lar core and not in the outermost layers 
of the star (e.g., lMaeder and Meynelll2005h . 

In recent years, considerable effort has been spent to measure surface mag- 
netic fields, by means of spectropolarimetry (exploiting the circular polar ization 
of the Zeeman compon ents) and a so-called least square deconvolution d Semell 
ll989LlDonatietallll997h which allows to measure the longitudinal component of 
the S-field averaged over the stellar disk. Though large samples of OB-stars (> 
45) have been an alyzed (see the m aterial presented by Sch nerr 2007, and refer- 
ences therein and IPetit et alll2008h . significant field strengths above the present 
detection threshold (1-a: 40-100 G) were detected for only ten stars (excluding 
chemically peculiar Ap/Bp stars), where in most cases the observations could be 
fitted by an oblique dipole (but see iDonati et all 12006 for the complex magnetic 
configuration in T Sco). Among these ten stars, there are three OB-stars with well 
known peculiarities in their spectra (0^ Ori C (04-6V), HD 192612 (06.5f?pe- 
08fp) and t Sco (B0.2V)), with polar field-strengths from 500-1500 G (assuming 
a dipole), and thre e 15 Cep stars which in cidentally were found to be nitrogen rich 
dMorel et alll2006h . rVerdugo etail ( l2003h rep ort the non-detectio n of surface mag- 
netic fields in a sample of 12 A-supergiants. iDonati et all ( 1200^ ) conclude that in 
non-peculiar hot stars, B is weak and/or acts on small scales (spots?). 

As should be clear from these introductory remarks, magnetic fields might 
be large in the stellar interior, but only few "normal" OB-stars show indications 
of strong surface fields, at least on larger scales. Even weak fields, however, can 
affect the wind, as we will see in the following. We will report on recent progress 
on magnetic winds, at first for a simplified case (no rotation) and then for the 
problem of field aligned rotation and the limiting case of an oblique rotator with 
very high field-strength. 



Magnetically channeled line-driven winds - The confinement par ameter. In order 
to expla in the observed X-ray emission from the AOp star IQ Aur. lBabel and Montmerld 
( Il997a f) suggested that the confinement of the wind by the strong magnetic field 
leads to a collision from the wind components of the two hemispheres in the 
closed magnetosphere and thus to the formation of a strong shock. Subsequently, 
iBabel and Montmerld ( fl997b) used this "magnetically confined wind shock (MCWS) 
model to explain the periodic variations of the X-ray emission from the 0-star 
d ^ Ori C (see above), requiring a polar field strength at the stellar surface of w 
300 G (present knowledge: w 1100 G. IDonati e"taiir2002) . where the variability 
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should be caused by a circumstellar disk predicted by the MCWS model. In these 
investigations, a prescribed, ^ec/ magnetic fie ld configuration had be e n used . 

This approximation has been relaxed by -Doula and Owockil (12001 see 
also lOwocki and ud-Doula 2004 for a comprehensive analytic description), who 
investigated the interaction between S-fields and winds in detail, based on 2-D 
time-dependent magneto-hydrodynamics for the simplified case of no rotation and 
a purely radial line-force. Thus, the external forces in the momentum equation are 
(effective) gravity, line-acceleration and Lorentz-acceleration, (47rp)^' (V x B) x 
B, and one has to solve for the S-fiel43 in parallel with the wind-dynamics . In the 
following, we will summarize the major findings bv lud-Doula and Owockil ( |2002|) 
and subsequent work. 

As it turned out, the most important parameter controlling the structure of the 
modified wind, T]* (see below), is related to the ratio of magnetic to wind energy, 

' ■£wind'^pvV2 Mv„ {l-R./r)P' 

assuming a j3-velocity field for the wind and B{r) °^ (7?*/r)^ (c-g-, ^ = 3 for a 
dipole field). The scaling factor in front of the last expression, evaluated at the 
magnetic equator under the assumption of a dipole configuration (5dipoi(^*j 6) = 
Sp(cos^ + 1 /4sin^ 0)^/^) has been denoted as the "confinement parameter", 

_ B^e = 90°)/?.^ _ (gp/2)27?,2 ^ B\,,R\, ^ ^^^^ 

In the latter expression, the polar field strength is measured in units of 100 G, the 
stellar radius in units of 10 Rq, the mass-loss rate in units of lO^^Moyr"' and the 
terminal velocity in units of 1000 kms^'. 

As we will see below, a value of 77* = 1 relates to the onset of larger struc- 
tural modifications due to B. To reach this, e.g., for the prototypical 04(f) super- 
giant C, Pup (RiQ w 2,M_6 w 4, vg = 2), a field of Bp w 320G would be required, 
whereas for the sun (where the magnetic field has a dominating influence on the 
wind) one finds (with M_6 w 10~^, vg = 0.5,6p w 1 G) a considerable value for 
the confinement parameter, T]* (0) w 40. 

For so-called weak winds (Sect. |5]), on the other hand, a magnetic field of 
only Sp < 30G (well below present detection limits) is needed to reach t]* = 
1 (assuming M < 10~^ M0yr~' and Voo= 2000 kms"'). We note that T Sco (a 
prototy pical weak-w inded star) reaches 77* w 80 for the measured field strength of 
500 G (IDonati et alll2006l) . i.e., extreme effects are to be expected. 

The reason to denote r|^, as a "confinement" parameter relates to the Alfven 
radius of the wind. Since MHD waves propagate with the Alfven speed, 

B V 1 

^A = 7:; — 777? ^ Ma = — = — (34) 
(47rp)V2 VA ^/r| 

the Alfven radius can be calculated from the corresponding Mach number Ma (/?a) = 
1, i.e., from the condition T](r, 6) = 1. For a dipole field, one obtains 

Ra^R* (t?*<1), RA^r]*^^R* (t?*>1). (35) 

V • B = and dB/dt = 'V x (v x B) in the MHD approximation of infinite conductivity. 
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Since the Alfven radius c orresponds roughly to the rn aximum radius of closed 
loops (see the discussion in lud-Doula and Owockill2b02h . the confinement param- 
eter controls this size, i.e., controls whether the wind is confined in such loops. 

Wind- structure as a function 0/7}^,. In their calculations. lud-Doula and Owockil 
(l2002h derived the wind-structure from evolving a model with initial conditions 
of a dipole field superposed upon a spherically symmetric wind. In dependence of 
the confinement parameter, three exemplary topologies have been found: 

(i) For moderately small confinement, T]* = 0.1, the surface magnetic fields are 
extended by the wind into an open, nearly radial configuration. There is still a no- 
ticeable global influence from B on the wind, enhancing the density and decreasing 
the flow speed near the magnetic equator. 

(ii) At intermediate confinement, 77* = 1, the field lines are still opened by the 
wind, but retain near the surface a significant tilt which channels the flow towards 
the magnetic equator. The polar velocity components can become significant, of 
the order of a few 100 kms~', and a thin disk begins to develop. 

(iii) For strong confinement, rj^ = 10, the field remains closed in loops near the 
equatorial surface, and the wind is accelerated and channeled upwards from foot- 
points of opposite polarity. These flows collide near the loop tops, with shock ve- 

locityjumps of up to 1000 kms^', leading to hard X-ray emission in t he keV band. 

These fi ndings are in qualitative agreement with the MCWS model bv lBabel and Montmerld 
(Il997ah . a nd also in more q uantitative agreement with X-ray observations for 

Ori C dGagne etalll2005l and Sect. 16.5.11 ). Since the shocked material at the 
loop tops is very dense, its support by magnetic tension becomes unstable, lead- 
ing to complex infall patterns along the loop lines down to the star. Even for large 
77* , the faster radial decline of the magnetic energy (compared to the kinetic en- 
ergy of the wind) leads to a final dominance of flow effects, and the field lines are 
extended into an open configuration at larger radii. 

Globally, the mass flux of the outer wind increases towards the equator due to 
the tendency of the field to divert the flow into this direction, whilst the base mass 
flux in equatorial regions (and bey ond for large t]*) is lower than the field-free case 
(see lOwocki and ud-Doulall2004l for further details). The terminal velocity, on the 
other hand, is larger than for a non-magnetic wind at all latitudes except for the 
equatorial plane, where it becomes significantly reduced due to the higher mass 
flux. Similar to the case of rotating winds, the to tal mass flux is barely affected as 
long a t]h, ^ 10. For further scaling relations, see lud-Doula et al (120081) . 

In summary, non-rotating magnetic winds develop a geometrically thin, slow 
and dense disk, superimposed by a fast, thin polar wind, quite different from ro- 
tating non-magnetic winds. Non-radial line forces should be irrelevant since the 
polar velocities are much larger than in the WCD case. X-rays are to be expected 
from the channeled flows colliding at loop tops and from shocks neighboring the 
compressed disk. Finally, it might be speculated whether oblique magnetic rota- 
tors can explain part of the observed UV variability (Sect. 16.1b and induce CIRs 
(Sect. 16.41) . because of the large density/velocity contrast which would occur be- 
tween the magnetic and rotational equator in such stars. 



Field alig ned rotation and the rig idly rotati ng Magnetosphere rn odel. In a follow- 
up study, lud-Doula et all (l2008l . see also lOwocki et all 12005b included rotation 
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(aligned with the magnetic field to avoid 3-D calculations) and investigated the 
question whether magnetic fields could spin up the stellar wind o utflow into a 
"magnetically torqued disk", as suggested by 'Cassi nelli et all (I2002D . The closed 
loops present for larger T]* tend to keep the outflow in rigid body rotation (i.e., 
Vij>{r) °= r), and it might be possible to propel material into a Keplerian disk if 
the Alfven radius is somewhat larger than the Keplerian co-rotation radius, /?k ~ 
(vrot/vorb)"^^"'^*; where rigid-body rotation would yield a balance between cen- 
trifugal and gravitational acceleration. Vorb = s/GMJrI is the orbital speed near 
the equatorial surface Note that strong confinement and rapid rotation are re- 
quired to obtain /?a > ^k- Although various combinations of rotational speeds and 
magnetic field strengths have been considered, and the analysis has been extended 
to 77* = 100, no stable Keplerian disk has been found to be formed in any of these 
cases, simply because most of the material does not have the appropriate velocity 
for a stationary orbit. However, there were clear indications of a "quasi-steady" 
rigid-body disk, as discussed further below. Note also that in the opposite case of 
slow rotation (more precisely, for /?a < Rk)^ the basic results as obtained for the 
non-rotating case remain valid. 

For the limiting case of very high field strength, t]* °o. lTownsend and Owockil 
([2005) developed the alternative "rigidly rotating magnetosphere" model, origi- 
nally designed to investigate the prototypical Bp star a Ori E with B 10^ G. In 
this model, the field lines behave like rigid tubes, and the outflowing material is 
constrained along corresponding trajectories fixed by t he prescribed field geome- 
try (see also the rigid field hydrodynamics approach bv lTownsend et al.2007,) . 

For sufficient rotation then, minima in the effective potential along each field 
line should gradually fill with plasma, forming an almost steady, co-rotating mag- 
netosphere. For the oblique rotator a Ori E, this should occur in the form of 
two co-rotating clouds. Time-series of synt hetic H^^t profiles for s uch a geometry 
are consistent with the observed behaviour dTownsen d et al"2005'). As mentioned 
above, the time-dependent calculations for large 77* by ud-Doula et al (2008) have 
confirmed the basic picture of a rigid-body disk (for aligned rotation), but with 
a complex dynamic behaviour with infall and outflow limiting the growth of the 
disk. E.g., the eventual centrifugal breakout of such material might suggest a new 
heating mechanism (via magnetic reconnection) to explain the har d X-ray flares, 
as obs erved in Bp-stars CGroote and Schmitti2004) and simulated bv iud-Doula et all 
(120061) . 



3 Stationary models of radiation driven winds 

3 . 1 Predictions from the improved CAK approach 

As was noted earlier, the origina l CAK a pproach to compute the wind struc- 
ture has been massively improved ([Friend and Abbott 1986, Pauldrach et al 198^ 
lPauldrach|[l987[ iPauldrach et alll 1994 l20ol ! Whilst the CAK models were com- 
puted in local thermodynamical equilibrium (LTE), departures from LTE have 
been included because of the strong influence of the radiation field and the low 



This approximation for /?k assumes the star to be far away from the QV hmit. 
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densities in the atmospheres of hot massive stars. Another significant improve- 
ment involved the inclusion of millions of lines (predominantly from iron-group 
elements) from some 150 different ions - in contrast to the single Cm list used 
in the original work of CAK. Within the current state-of-the-art modified CAK 
modeling, the rate equations, radiative transfer equations, energy equation, and 
approximate hydrodynamic equations, using force-multiplier parameters, are all 
accounted for. More simplified models relying on an approximate solution for 
the rate-equations and an analytic approach to solve the hydrodynamic equa- 
tions (with consistent, depth dependent force-multiplier parameters) have been 
presented by lKudritzkil(l2002h . 

A complementary approach currently in use concerns the hydrodynamical 
method of Krticka (2006). As in the models by Pauldrach and co-workers ("WM- 
basic", see Sect. 14.2.11 . the various equations are consistently solved, but the con- 
cept of force multipliers is dropped, though the Sobolev approximation is still ap- 
plied. The major advantage of this method is that it allows for a multi-component 
description of the hydrodynamical equations, which enables an investig ation of 
potent ial metal-ion decoupling from the H/He plasma. A drawback of the lKrtickal 
models is that line-blocking and multi-line effects are not accounted for - 
limiting their applicability to weaker winds. 



3.2 Predictions using a Monte Carlo radiative transfer approach 

A diffe rent approach involves the "Monte Carlo" method developed by I Abbott and Lucvl 
(Il985h . where photons are emitted at an inner boundary, and their scattering his- 
tory is tracked on their journey outwards. At every interaction some amount of 
momentum and energy is transferred from the photon to the ion under consider- 
ation. From a macroscopic point of view, the overall integrated mass-loss rate is 
obtained from global energy conservation: 



where A L is the total amount of energy that has been transferred per second from 
the radiation field to the scattering ions in the wind. The main drawback of this ap- 
proach is that one first needs to establish an appropriate velocity law, v(r), gener- 
ally using an empirically pre-determined Voo, although in pri nciple it is possible to 
iterate and obtain dyna mically consistent mass-loss rates (see lSpringmann and Puis 
ll998LIVinket"aill 19991) . 

One of the major advantages of the Monte Carlo method is that it easily al- 
lows for multi-lin e scattering, which becomes important in denser winds, as al- 
ready outlined by iFriend and CastoJ (Il983l but see also rail987h . Before the 
year 2000, theoretical mass-loss rates fell short of the observed rates for dense O 
star and WR winds, whilst for weak winds the oft-used single line approach over- 
estimated the mass-loss rates con siderably. Although early Monte Carlo simula- 
tions (e.g. lAbbott and Lucy||l985h only treated line scatt ering, this approximation 
has been relaxed in the mass-loss rate determinations bv lVink et al (2000, 2001), 
although further improvements to the method, e.g., in terms of line branching (iSirnI 
12004') and wind clumping need to be considered in future model generations. 




Vesc^) = 



(36) 
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Here we wish to mention that the high mass-loss rates derived for dense winds 
when accounting for muhi-line scattering (as in the Monte Carlo method) do not 
violate the principle of momentum conservation. An appealing picture is that of 
"hemispheric scattering" where a single photon can be launched in one direction, 
transfer its momentum to an ion, and travels towards the opposite hemisphere - 
repeating the same process over and over again, thereby pumping up the momen- 
tum transfer in the wind. One may note that all what really happens is that at each 
iteration the photon is slightly red-shifted, until it eventually runs out of energy. 
In reality, w e know that hemispheric scattering w ill not occur. It has been demon- 
strated (e.g.. lSpringmannll994lGayley et alll995h that when one employs realistic 
line lists, photons follow paths that are less extreme, but the essential point is that 
multiply scattered photons are always adding radially outward momentum to the 
wind, and the wind momentum may easily exceed the so-called single-scattering 
limit, i.e., Mv„/{L/c) can become larger than unity. 

I Vink etall(l200(]L 1200 Ih used the Monte Carlo method over a wide range of 
stellar parameters and derived a "mass-loss recipe" using multiple linear regres- 
sions to their model results. On the hot side of the predicted bi-stability jump at 
^ 25 000 K (but see Sect. l4.2.Tl) . the rates of OB supergiants should roughly scale 
as 

M oc l2.2 ^-1.3 (v„/Vesc)"^ (37) 

when L is measured in 10^ Lq, M in 30 Mq and T^ff in 40 kK. This relation 
predicts that the mass-loss rates should strongly depend on the luminosity (L?'^), 
steeper than M oc i'-^ which is often quoted and relies on the simplified scaling 
relation Eq.[T4]with a' w 0.6. The reason is that the efficiency of multi-line effects 
accounted for in Monte Carlo simulations increases with increasing wind density, 
i.e., for more luminous stars, whereas in (modified) CAK models these effects 
are normally neglected, particularly in the scaling relations presented in Sect. 12. II 
Because the mass-loss rate also scales with stellar mass as M^'^, and applying 
a typical M — L ratio of L M*^ with x = 2 - 3 for massive stars, one obtains an 
overall M scaling as /^i 55-1.77^ agreement with observational findings. 

The basic succe ss of the Mo nte Carlo method is exemplified when comparing 
Figs. 1 and 3 from lVinfl (1200 6*), which display the degree of agreement between 
modified CAK models and observations on the one side and Monte Carlo models 
and observations on the other: by properly including multiple scatterings in the 
predictions, the results are equally successful for relatively weak (with M ^ 10~' 
M0yr~') as well as denser winds (with M ^ 10~^ M0yr~'). The predictions can 
be conveniently expressed via the WLR. For 0-stars hotter than 27 500 K, the 
relation is shown in Fig. |5] and given by Eq. (116b with x = 1.83, corresponding 
to a' = 0.55. For objects on the cool side of the bi-stability jump, the mass-loss 
rate increases and the slope of the WLR becomes x = 1 .9 1 , corresponding to a' 
= 0.52. 

M onte Carlo predictions for hot stars of lower mass were computed bv lVink and Cassisil 
(|2002h , whilst modified CAK predictions for Central Stars of Planetary Nebulae 
were computed by Pauldrach et al. (1988, 2004), see pagelSTl 
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Fig. 5 Predicted WLR for O stars hotter than 27 kK for a range of (L,M) -combinations in the 
upper HRD. From lVink etail ( l2000h . 



3.3 Predictions for models close to the Eddington limit 

The predictions of Eq. (|37] | are valid for models that are located at a sufficient dis- 
tance from the Eddington limit, with F < 0.5. There are two regimes where this 
is no longer the case: (i) stars that have formed with large initial masses and lumi- 
nosities, i.e. very massive stars (VMS) with M ^ 100 Mq, and (ii) less extremely 
luminous "normal" stars that approach the Eddington limit when they have lost a 
substantial fraction of their initial mass. Examples of such objects are the LBVs 
and the more common evolved WR stars. 

IVink and de KoteJ (|2002|) and ISmith et all ( |2004) showed with their Monte 
Carlo computations that when lowering the mass for constant luminosity, the 
mass-loss rate increases more rapidly than Eq. (137b indicates. For these LBV mod- 
els, M M"'^ instead of M~'^ as for normal OB supergiants. In other words, 
when stars start to lose mass, not only does the mass-loss rate increase due to a 
lower stellar mass, but as the mass-loss rate increases more rapidly this leads to a 
strong positive feedback. 

An illustration o f this feedba ck-effect for VMS in the mass range 100 — 300 
Mq was provided by IVinkI (l2006h where it was noted that VM S mass loss drasti- 
cally increases when the objects approach the Eddington limit. IVink and de Koted 
(|2005) also simulated mass loss close to the Eddington limit when performing a 
pilot study for evolved late-type WR (WN and WC) winds, as a function of metal- 
lici ty. They showed that these mass-loss rates were a factor 6-7 higher than using 
the lVink et ai (l2000l) recipe for OB supergiants - for objects with similar stellar 
parameters, as a result of the closer proximity to the Eddington limit. 

A different appr oach to the problem of winds close to the Eddington limit was 
employed by Grafe neT and Ham^itil (|2008|) . One should note that their WR wind- 
models were not applied to evolved WR stars, but to WR stars in which the broad 
emission lines form due to a large intrinsic luminosity. For WR wind-models, the 
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Fig. 6 Mass-loss predictions as a function of Z, for objects with Z/Zq > 1/100. The solid 
line indicates the M(Z) dependence from M onte Carlo simula tions by I Vink et all i2001[) and the 
dotted line corresponds to the predictions bv lKudritzkH i2002D . See text. 

CAK formalism is no longer applicable, as a myriad of weak iron lines provide 
massive line overlap, and the wind has become optically thick. 

Th e fra mework of optically thick winds has been employed in studies of lKatol 
(1985*) and lNugis and LamersI (l2002h where the basic assumption is that the criti- 
cal poin t of the problem is the sonic point, similar to recent suggestions by iLucvl 
(l2007allbh regarding "normal" 0-star windQ, and that the corresponding condi- 
tions can be described in the diffusion limit. In the transonic region, the Rosseland 
mean opacity shows an outwards increasing behaviour and, through the use of the 
equation of continuity, M follows from the specific condition s at th e sonic point. 
In the more advanced studies of iGrafener and HamannI 120051 1200^") the momen- 
tum equation of the entire wind is accounted for and a strong dependence of M on 
the proximity to the Eddington limit is also predicted for these models (see also 
Sect. 14331) . 



3.4 Predictions for very low Z and Pop III stars 
3.4.1 The observable massive stars at low Z 

IVink et ail ( l200lh applied their Monte Carlo models to OB-stars within a Z-range 
representative for the observable Universe, Z /Z,^ > 1/100. The corresponding re- 
sults are compared to those from lKudritzkil (l2002h in Fig.|6] Over the last decades 
themodified CAK models predicted M Z'" with m in the range 0.47-0.94 dAbbottl 
119821 iKudritzki et al 1987. Kudritzki.2002) . The main reason for this scaling is 
that the CAK k parameter is strongly Z-dependent, whilst second order effects 
that relate to how a varies with Z have a larger effect on the terminal wind veloc- 
ity instead (see pages[TTJf.). 



A discussion of this assumption is beyond the scope of this review. 
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As can be noted from the fi gure, the Vink et aj (l200lh predictions exhibit a 
steeper slope {m = 0.85) than the lKudritzkiT i 20021) predictions, with m = 0.5 - 0.6. 
Again, the reasons are thought to be related to the neglect of multi-scattering in 
CAK-type models, i.e., the corresponding single-line approach could result in both 
an M — L dependence and an M-Z dependence that are too shallow. On the other 
hand, most results computed using the Monte-Carlo method have not accounted 
for a Voo-Z dependence, which is anticipated to result in an overestimate of m. 
Therefore, a value of in = 0.7 has been rec ommended as th e appropriate value 
for the mass-loss scaling with Z for OB-stars dVink et alll20 0l L[Krticka 20061. The 
observed metallicity dependence of OB-stars will be discussed in Sect. 14.2.1] 



3.4.2 Objects at very low Z and Pop III stars 

Extending the predictions to very low Z/Zq < 10~^, where stars can no longer be 
observed individually, it has been shown that M keeps dropping until the winds 
reach a regime where they become susceptible to ion-decoupling and multi-com- 
ponent effects dKrticka et alll2003h . The only way to maintain a one-fluid wind 
model is through an increase of the Eddington factor by pumping up the stellar 
mass and luminosity. Such an increase in stellar mass is by no means artificial, as 
the earliest generations of massive stars are anticipated to be more massive than 
today's population since stars formed with fewer metals to provide the cooling - 
at the very least resulting in a larger Jeans mass. 

In case of Pop III stars with truly "zero" metallicity, i.e. only H and He present 
in the atmosphere, it seems unlikely that these stars would develop line-driven 
winds of significant strength (lKudritzkill200llKrticka and Kubatll2006l) . Nonethe- 
less, other effects may contribute to the driving. Interesting possibilities include 
stellar rotation and pulsations, although pure vibr ation models for P op III stars 
also indicate little mass loss via pulsations alone dBaraffe et alll200lh . Perhaps a 
combination of several effects would result in large mass loss close to the Edding- 
ton limit. Moreover, we know that even in the present-day Universe a significant 
amount of mass is lost in LBV type eruptions (potentially driven by continuum 
radiation pressure, see page l56b which migh t be also relevant for the First Stars 
dVink and de Koteij|2005llSmith and Owockill2006h . 



3.4.3 Winds enhanced by self-enrichment 



Given that the First Stars are potentially massive, luminous, and rapidly rotating, 
it is not inconceivable that, despite the fact that the first generation(s) of massive 
stars start their evolutionary clock with few metals, the objects may enrich their 
outer atmospheres with nitrogen and carb on due to rotational mixing tMevnet et all 
l2006h . inducing a strong line-driven wind dVinkl2006LrGrafener and Hamannl2 008). 

In a first attempt to investigate the effects of self-enrichment on the total wind 
strength, Vink and de Koter ( 2005) performed a pilot study of WR mass loss ver- 
sus Z. The prime interest in WR stars here is that these objects, especially those 
of WC subtype, show the products of core burning in their outer atmospheres. 

Although the l ast few decades provided increasing evidence that WR winds are 
radia t ively-driven dL ucv and Abbott 1993 , Gavlev et al 1995. Sprin gmann and Puis! 
ll998LlNugis and Lamers.2002 . ,Grafener and Hamann.2005i) . the question of a WR 



32 



J. Puis, J.S. Vink & F. Najarro 



-4r 




log (Z/Ze) 



Fig. 7 Monte Carlo WR mass-loss predictions as a function of Z. The dark line represents the 
late-type WN stars, whilst the lighter dashed line shows the results for late-type WC stars. The 
slope for the WN models is similar to the predictions for OB-supergiants, whilst the slope is 
shallower for WC stars. At low Z, the slope becomes smaller, whilst it flattens off entirely at 
Z/Zq = 10^"*. The computations are from lVink and de Koterl i2005D . 

Z-dependence remained controversial, with some stellar modelers extrapolating Z 
dependencies from O stars, whilst others assumed WR winds to be Z-independent. 

The reasoning behind the assertion that WR winds may not be Z-dependent 
was that WR stars enrich themselves by burning He into C, and it could be the 
large C-abundance that is the most relevant ion for the WC wind driving, rather 
than the sheer number of Fe lines. Figure |7] shows that despite the fact that the C 
ions overwhelm the amount of Fe, both late-type WN (dark line) and WC (light 
line) show a strong M-Z dependence, basically because Fe has such a complex 
electronic structure. 

The implications of Fig.Qare two-fold. The prediction that mass loss no longer 
decreases when Z/Z© drops below ^10^^ (due to the final dominance of driving 
by carbon lines) suggests that once massive stars enrich their outer atmospheres, 
radiation-driven winds might still exist, even if the object started its life with negli- 
gible amounts of metals. Whether the rate of mass loss is high enough to seriously 
alter the evolutionary tracks of the First Stars also depends on other physical fac- 
tors, such as the proximity to the Eddington- or the I2F-limit. 

The second point to address regarding Fig.|7]is the new result that the mass- 
loss rates of WR stars drop steeply with Z at subsolar Z. This is anticipated to have 
important consequences for black hole formation and th e progeni t or ev olution of 
long duration gamma-ray bursts. The collapsar model of lWooslevI ( I1993D requires 
a rapidly rotating core before collapse, but at solar Z stellar winds are anticipated 
to remove significant angular momentum (Zahn 1992). The WRM-Z dependence 
from Fig. |7] provides the intriguing possibility to maintain rapid rotation because 
WR winds are predicted to be weaker at lower Z. 

Interestingly, the WR M — Z depen dence also appears to be a ble to reproduce 
the drop in the WCAVN ratio at low Z dEldridge and Vinkll2005) . Although there 
are various reasons to assume that additional effects such as rotational mixing and 
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binarity play an important role for the formation of WR stars, the results suggest 
that the mass-loss-Z dependence is a first order effect as far as the WCAVN ratio 
is concerned. 



4 Observed wind parameters 

As we have seen in the preceeding sections, radiation driven wind theory provides 
precise predictions for the two global parameters of hot stellar winds, namely the 
terminal velocity, Voo, and the mass-loss rate, M, as a function of stellar parameters 
(including abundances, rotation rates, magnetic field-strengths). The most promi- 
nent prediction concerns the presence of a wind-momentum luminosity relation (if 
"contaminating" effects do not dominate), which should allow for an easy check 
of the theory when the global wind parameters are known from observations. 

Later on, we will summarize recent results on such comparisons. In order to 
do so, we need to know how and to which precision these wind parameters can be 
"measured". 

A few caveats. Most of the following investigations (still) rely on the assump- 
tions of the standard model (Sect. lZTI) . partly extended to allow for the presence of 
inhomogeneities (clumping) and X-rays (see below). Within this model, the mass- 
loss rate follows from the continuity equation (Eq.[Tl, and practically all diagnos- 
tics base on the approach that certain observational features (spectra or SEDs) are 
fitted by corresponding radiative transfer calculations on top of a wind model with 
a density expressed in terms of the mass-loss rate, p(r) = M/{Anr^v{r)), and a 
/3 -velocity field. 



Remember that this functional form relies on the theoretical predictions (Eq.ll2b. 
The constant b fixes the onset of the wind, usually obtained from requiring a 
smooth transition between photospheric density-Zvelocity-stratification and wind 
conditions, or by imposing a certain minimum velocity v{R^) = Vmin, of the order 
of the isothermal speed of sound. 

By varying M,Voo and j3, an optimum fit is aimed at and the so-called "ob- 
served wind parameters" are just the outcome of such fitting procedures. Insofar, 
already their determination relies on a consistent description of the wind and the 
photosphere, since not only the wind properties but also additional atmospheric 
conditions such as Teff, logg and abundances affect the radiative processes which 
are modeled to fit the observations. Whenever "observed wind parameters" are 
presented, one should be aware of the fact that already these are the result of 
diagnostic techniques based on a substantial amount of theoretical modeling. De- 
pending on the degree of sophistication of the underlying model but also on the 
part of the spectrum and the specific process which is used for the determination, 
the reliability of the results can vary drastically. 

Even worse, by investigating the scaling properties of particular solutions, it 
turns out that such fits are not unique, but depend, subject to the considered pro- 
cess, on certain combinations of wind parameters and stellar radius. Important 
quantities are the "optical depth invariants" Qres for resonance lines with a line 




(38) 



34 



J. Puis, J.S. Vink & F. Najarro 



opacity o= p and Q for recombination based line processes (such as Ha) with opac- 
ities o= p^, 



M _ M 



Gres = 77-^. Q = TTT-TTJ (39) 



(details are given in the following). Individual mass-loss rates can be obtained only 
when the terminal velocity and the stellar radius (depending on distance which is 
problematic for Galactic objects) are known, or, at least in principle, if resonance 
and recombination processes are fitted simultaneously by means of NLTE atmo- 
sphere codes. Otherwise, the result of such fits is not the mass-loss rate but one of 
the above scaling invariants. Thus, even if the fit looks perfectly, the same needs 
not to be true for the derived mass-loss rate. 

The standard model is certainly an (over?)-simplification. Significant evidence 
for deviations from spherical symmetry, non-stationarity, dumpiness and shocks 
has been found from the beginning of wind observations on. Aspects related to 
rotation and magnetic fields have already been discussed in Sect. 12.21 and we 
will turn to problems related to time-dependence and inhomogeneity in Sect. |6l 
Though these phenomena had been widely ignored with respect to diagnostical 
methods until the late 1990's, increasing observational and computational feasi- 
bilities have lead to the conclusion that certain aspects must be included to avoid 
obvious inconsistencies. In due course, simplified approaches have been devel- 
oped to account for these aspects, where the inclusion of wind-clumping and the 
X-ray/EUV emission from shocked material into the present generation of NLTE 
codes are the most prominent examples of such updates. Note, e.g., that the optical 
depth invariant Q as defined in Eq.[39]needs to be modified with respect to the so- 
called clumping factor if the winds are considered as inhomogeneous (Sect. 16.31 
Eq.m. 

Obviously, a detailed inclusion of deviations from the standard model requires 
significant effort regarding both diagnostics and the development of new models 
and methods, a task which is just at the beginning to date. Examples are the devel- 
opme nts of multi-dimensional NLTE codes such as ASTROROTH ( Georgiev et all 
l2006h mentioned already in the context of B[e] stars (Sect. l2!2Tt and a 2-D de- 
scription of rotational effects within the spectrum synthesis based on NLTE model 
atmospheres, as reported by ,Bouret et aL (.2008.') . Further progress can be expected 
in the near future. 



4. 1 Diagnostic methods 

Diagnostic methods can be roughly divided into two different classes. Approxi- 
mate methods concentrate on the analysis of one specific process and try to de- 
duce specific (wind-)parameters by a simplified description, usually by consid- 
ering the wind physics alone and requiring additional information from external 
photospheric diagnostics. 

Alternatively, modern diagnostics rely on the use of "unified" NLTE model at- 
mospheres. Stellar and wind parameters are derived simultaneously by optimizing 
synthetic spectra and SEDs for a large wavelength range, comprising diagnos- 
tics of a multitude of processes with different dependencies and scaling relations. 
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In order to draw correct conclusions, these processes have still to be understood 
though. 



4.1.1 Approximate methods 



A detailed discussion and comparison of the v arious approxim ate methods to de- 
rive global wind parameters has been given by lKudritzki and Puis, (.2000.) . and we 
will summarize only important aspects. 

The two most prominent line types formed in a stellar wind are P Cygni pro- 
files with a blue absorption trough and a red emission peak, and pure emission 
profiles (or absorption lines refilled by wind-emission). Their different line shapes 
are caused by different population mechanisms of the upper (emitting) energy 
level of the transition. 

In a P Cygni line, the upper level is populated by the interplay between ab- 
sorption from and spontaneous decay to the lower level, usually the ground state 
of an ion. This process is called line-scattering. If, on the other hand, the upper 
level is populated via independent processes, e.g., by recombination into this level 
or decay from levels above, a pure emission line (or an absorption line weaker 
than formed by purely photospheric processes) will result. The former line type 
as seen in the (F)UV is majorly used to determine the velocity field in the wind, 
particularly Voo, and the latter (H^, but also He li 4686, Br^) to derive M, or, more 
precisely, Q. 



Analysis of UV P Cygni profiles. (F)UV P-Cygni lines from hot stars (e.g., C IV, 
N V, Si IV, O VI, P V) are usually an alyzed by means of the so-called "SEI" method 
(Sobolev plus exact integration, see lLamers et allfl987 h. based on a suggestion by 
iHamann (1981a). The central quantity which controls the line formation process 
is the (radial) Sobolev optical depth of the specified line (cf. page|7]i 

TsobW = §^^, X{r) = —Mir). (40) 

dv/ar Voo m^c 

f is the oscillator-strength and «i the lower occupation number of the transi- 
tion, neglecting stimulated emission. Here and in the following, radii r and ve- 
locities V are measured in units of the stellar radius and the terminal velocity, 
respectively (remember that the line optical depth in a stellar wind is an almost 
purely local quantity, cf. page|7]). Relating the occupation number, ni, to the local 
density, this quantity can be expressed by 

^ /\ n\ 171 \ t \ ^ {ne'^)/{m^c) Ak 

Tsob y = 2 A lA ^ Kr)=E{r)q{r)- j — fX (41) 

r^vav/ar R-tVoo^ 4nmu l+4yHe 

where E is the excitation factor of the lower level (= 1 for resonance lines coupled 
to the ground-state), q the ionization fraction, the abundance of the element 
with respect to hydrogen and Fhc the helium abundance. Obviously, this expres- 
sion is scaling invariant with respect to the quantity Qres = M/ (7?* Voo^) (cf. Eg. 1391) 
as long as the ground-state population of the considered ion is proportional to p 
throughout the wind, e.g., if the ion is a major one. From the latter equation and 
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Fig. 8 SEI fits to the P Cygni profiles of O vi, P v, N v and C iv of ^ Pup (041(f)) as observed by 
COPERNICUS <Morton an d Underbill 1977). An almost unique solution with = 2350 kms^' 
and j3 = 0.8 is found for all four profiles, with a "micro-turbulence" of 10% of v^. The upper 
panels show th e derived run of the quantity k{v) <x Mq (Eq. 1411 . Note the black troughs in N v 
and C IV. From lHasej l fT99l . 



the assumed j3 -velocity law, it should be clear that the maximum set of deducible 
parameters from a fit to the observations consists of (v„, j3, k{r)), see Fig. [8] 

Thus, information on the mass-loss rate is hidden in the quantity k{r), and 
mass-loss rates can be derived from the P Cygni profiles of resonance lines (if at 
all) only when external knowledge about ionization, abundance and stellar radius 
is available. Exemplarily, Howarth andPrinja (1989) and Haser (1995) inferred 
the product M{q) for large 0-star samples from unsaturated profiles , where {q) 
is a suitable (spatial) average of the corresponding ionization fraction. 

Most of the "strategic" P Cygni lines, however, are saturated under typical con- 
ditions (C IV, N V in early 0-stars and Si IV in late 0-supergiants). In these cases, 
the derivation of mass-loss rates (or more precisely, of the product M{q)) becomes 
impossible, and only lower limits are accessible. The reason that these lines are of- 
ten saturated has to do with the particular dependence of the line optical depth on 
the product [r'^vAv / Ar)~^ , which, as already stated on page[lOl varies according to 
y(^^) 1/^-2 Pqj. jjjjgg a rather constant ionization fraction throughout the wind 
and typical 0-star parameters, this variation is only mild, and a line most likely 
remains optically thick throughout the wind (i.e., the P Cygni profile is saturated) 
when it has a large optical depth already in the sonic region. In this case, "only" 
terminal velocities (from the frequency position of the blue absorption edge) and 
veloc ity field parameters B (from the shape of the em ission pe ak) can be measured 
(e.g., iHowarth and Prinjall 19891 iGroenewegen etalll 19891 and lHasedll995h . Typi- 
cal values j3 = 0.7 ... 1 .5 have been found for OB-stars, where supergiants show a 
clear tendency towards higher i3. For A-supergiants, values a s high as j3 = 3 . . . 4 
have been derived (e.g.. lStahl et allll99lLlKudritzki et allll999h . 

The most severe problem of P Cygni line diagnostics is related to the presence 
of so-called black troughs, which are extended regions in the absorption part of 
saturated profiles with almost zero residual inte nsity (see Fig. [Hi. These troughs 
cannot be explained without further assumptions. iHamannl ( Il981at ib) introduced a 
highly supersonic "microturbulence" (of the order of 10% of Voo) throughout the 
entire wind to overcome this and related problems. Since in this case the intrinsic 
absorption profile is extremely broad, the observed troughs can be simulated, in 



Mass loss from hot massive stars 



37 



paral lel with the red-shifted emission peaks (see also iGroenewegen and LamersI 
ll989land |Haser 1995| for the inclusion of a depth-dependent velocity dispersion). 

Lucvl (fl982. ,1983h suggested an alternative, more physical explanation for 
the generation of black troughs due to enhanced back-scattering in multiply non- 
monotonic velocity fields, which most likely are present according to observations 
and time-dependent hydrodynamical simulations (see Sect. |6ll. Detailed UV line 
formation calculations have shown that black troughs can b e actually created in 
the profiles of such models ( Puis et allll993allOwockilll994bl) . 

Note finally that until very recently UV resonance lines have been considered 
as being (almost) uncontaminated by clumping effects. As we will see in Sect. l6.3l 
this might not be true, due to the effects of porosity or "vorosity" (= velocity 
porosity). 



Ha diagnostics. To avoid the problems inherent to the M-determination from 
UV P Cygni lines (which, moreover, require space-bound observations), the com- 
monly accepted standard diagnostics to derive M relies on the wind emission in 
Hcf, where the influence of ionization (which is almost complete) and abundance 
becom es almost negligible. Th e ide a to use Hg as a mass - loss indicator goes 



Klein and Casto^ (19781) andlLeithered (ll988h . lDrew] (Il990l) . iGabler et all 



back to ^ -. . , 

(Il99(lh . IScuderi et all(ll992h and lLamers and Leit herer^ C199J) applied this method 



to derive first results for a variety of hot stars. Pulsetal (1996) extended the 
method to obtain Z-fast analysis tool and to eliminate some systematic errors in- 
herent to previous approaches. 

The basic difference between Hg (and similarly behaving lines such as He II 
4686 and Brg) and P-Cygni lines is related to the source function, which in the 
former case is rather constant throughout the wind, as long as the involved levels 
are predominantly fed by recombinations or decays from upper levels. This con- 
dition is met in O- and early B-stars. For cooler stars, the population by pumping 
due to Lyman lines becomes equally important, until, in the A-supergiant domain, 
these lines and the Lyman continuum become optically thick and the correspond- 
ing transitions reach detailed balance. Consequently, the second level of hydrogen 
{n2, the lower level of H^) becomes the effective ground state of Hi and the line 
behaves as a scattering line, thus displaying a typical P-Cygni signature when 
the winds are sufficiently strong. In such cases, the profile is usually visible until a 
blue edge frequency corresponding to Voc, which allows to derive al so this quantity 
and thus enables a complete wind analysis in the optical alone (see lKudritzki et all 
[1992). 

Since neutral hydrogen is almost always a trace ion in the winds of OBA stars, 
the opacity of Hg scales with p^, 

, r U ^ \^ r\ X{r)lR* b2{r) ^^^^ 



i.e., the scaling invariant is (Eq. |39l ), and b2 is the NLTE departure coefficient 
of n2 which can be rather easily calculated^ Although similar arguments hold for 

When n2 becomes the effective ground state, the opacity scales with /W with dilution 
factor W, see lKudritzkT and Puls,.200a . 
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the formation of He II 4686, the corresponding departure coefficients have a more 
complicated (and more uncertain) stratification which favors the use of Ha- 

Scali ng relations for the equivalent width of in 0-stars have been pro- 
vided bv lPuls et all (Il996l) . and typical errors for Q are of the order of 30%, where 
these errors can reach values below 10% for strong wind emission. Also here, 
the velocity field parameter j3 can (and has to) be determined from the shape of 
the em ission peak or w ing, resulting in typical values j3 w 1 for 0-type super- 
giants (iPuls etall ll99^. For weaker winds, on the other hand, the wind emission 
is low and "hidden" within the (rotationally broadened) photospheric profile, and 
accurate determinations of M become more difficult. Detailed NLTE calculations 
based on unified model atmospheres (Sect. l4~L2] l are required to derive meaning- 
ful results in such cases. But even if doing so, uncertainties up to a factor of two 
in M can arise if the photospheric absorption profile is only marginally refilled. 
Such uncertainties are rooted in the strong dependence of the mass-loss rate on j3 
and the fact that this parameter can no longer be derived (because of the hidden 
wind emission ) but has to be a ssumed instead (see, e.g., the corresponding error 
analysis in ,Repolust et alll2004h . 

At even lower wind-densities (roughly below a few times 10"** MQyr~^ in the 
case of 0-/early B-stars), completely looses its sensitivity to mass loss, due to 
a vanishing wind emission. For such low M (which is present, e.g, in many SMC 
0-stars and so-called weak- winded stars, see Sect. [51, other diagnostics must be 
used, either UV resonance lines or Bra (see below). 

The major problem of Ha, however, is its dependence. Any significant 
degree of clumping will lead to an overestimate of the mass-loss rates if this is 
neglected in the analysis. Somewhat advantageous, however, is the fact that H^ 
remains optically thin in the largest part of the emitting wind-volume, such that 
porosity/vorosity effects can be neglected, contrasted to the analysis of (strong) 
UV resonance lines. Further details and consequences are discussed in Sect. 16.31 



Thermal radio and FIR continuum emission. A conceptually different approach 
to "measure" mass-loss rates is to exploit the information contained in the (F)IR, 
(sub)millimeter and radio continua. Actually, this approach provides a remarkably 
reliable and easy-to-use tool, since it is based on rather simple processes. 

The basic idea is to measure the excess relative to the flux predicted by photo- 
spheric models, which is emitted by free-free ("Bremsstrahlung") and bound-free 
processes in the wind. This excess flux becomes significant at longer wavelengths, 
due to the dependence of the corresponding opacities and the corresponding in- 
crease of the effective IR/radio-photosphere, and can be translated into a mass-loss 
rate. In the following, we will concentrate on purely thermal emission, whereas 
non-thermal effects are discussed further below. 

The approach to exploit the radio emission (here, the "excess" corresponds 
to the totally emitted flux, due to negligib le photospher ic emission) was i ndepen- 
dently developed by I Wright and Barlo w (197513 and Panagia and Felli (1975),_ 



and applied to larger 0-star samples bv lAbbott et aU198(lil981a) and .Lamers and Leithered 



(0993). A generaliz ation towards shorter IR-w avelen gths has been provided in a 
series of papers by Lamers a nd WatersI (Il984al lbh and lWaters and LamersI (Il984h . 
who expressed the flux excess in terms of a curve of growth, from which the 



these authors considered IR emission as well. 
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velocity law and the mass-loss rate can be derived simultaneously. A first appli- 
cation of t his method was th e analysis of IRAS observations at 12, 25 and 60 jim 
of ^ Pup ( Lamers et allll984) . Further references (regarding radio, submm and IR 
observations and analyses) can be found in [Kudritzki and Puis (2000). 

Following iLamers and WatersI (Il984ah . the combined free-free and bound free 
opacity (in units of cm~^) at frequency v can be written as 

= 3.692 • 10« (1 - exp(-^)) ? {g{v, T)+b{v, T))^, (43) 

with electron temperature T, the mean value of the squared atomic charge, 
g{v,T) and b{v,T) the Gaunt factors for free-free and bound-free emission, y 
the ratio between electron and ion density, and n,- the ion density, which relates to 
the mass density via p = nijimn with atomic weight jUmn- For long wavelengths 
and (almost) completely ionized H and He, the major dependencies are 

which increases strongly with A and p. Inserting typical values, it is easy to show 
that from the FIR on the continuum becomes optically thick already in the wind, 
and the emitting wind volume increases as a function of A. E.g., the radio photo- 
sphere of typical 0-supergiants is located at 100 stellar radii or even further out. 
Exploiting this knowledge, one can approximate v(r) w Voo at radio wavelengths, 
and an analytical solution of the radiative transfer problem becomes possible for 
an isothermal wind: 



F..^23.2(g)''^ (|)'^\ (45, 



when Fy is the observed radio flux measured in Jansky (10~^^ Wm~^ Hz~^), M 
in units of M0yr~\ Voo in kms^^, distance d to the star in kpc and frequency V in 
Hz. Thus, the spectral index of thermal wind emission is 0.6, since 

Fv- ((v^(v,r))''"ocvO^ (46) 

The scaling invariant of both radio and IR emission is M/R^ rather similar to 
the quantity Q, if the object's angular diameter d/R^ remains fixed during the 
analysis. For typical 0-stars located at 1 kpc, Fy is of the order of 0.1 mJy in the 
radio range (2 - 20 cm). Thus, OB-stars are generally weak radio sources (which 
prohibits the use of this method for extragalactic work!) and useful radio obser- 
vations can be obtained only for nearby stars with rather dense winds. Moreover, 
since the spectral index of a non-thermal component (see below) is negative (e.g., 
= —0.7 for synchrotron emission), the radio domain can be more easily contami- 
nated than shorter wavelengths such as the (sub)millimeter regime. 

Again, due to the scali ng with p^, the majo r pro blem is the impact of clump - 
ing, already dis cussed bv lAbbott et all (Il981al) and iLamers and WatersI (Il984bl) . 
IPuIs et all (120061) have used this additional dependence to derive the radial stratifi- 
cation of the clumping factor, see Sect.| 



40 



J. Puis, J.S. Vink & F. Najarro 



Non-thermal (NT) emission. So far, we have concentrated on the thermal emis- 
sion. However, there is ample evidence of non-therm al emission in hot stellar 
win ds. First observation al findings have been reported by' White and Beckeil (Il983h 
and lAbbottetall ( Il984l) , andfBieging et al ( 1989) concluded that roughly 30% of 
all massive stars were non-thermal em itters. To date, there are 16 0-type stars con- 
firmed as non-thermal radio sources dvan Loo et alll20d^ and references therein). 
A recent review on non-thermal emi ssion from massive stars (with emphasis on 
massive binaries) has been given by 'De Becker' (|2007|), which includes the most 
up-to-date census of non-thermal radio emitters (the aforementioned 16 0-stars 
plus all known WR sources) including their multiplicity status. 

Although it seems to be clear that the non-thermal emission is due to syn- 
chrotron radiation from relativistic electrons, different agents responsible for the 
acceleration are di scussed in the literature: wind accretion onto compac t objects 
dAbbott et al"l984'), magnetic reconnection in single or colliding winds (Pollocld 
[1989, see also Litvinenko 2003) and first-order Fermi acceleration (Fermi 1949), 
which, in the presence of strong hydrodynamic sh ocks, is also referred to as the 
Diffusive Shock Acceleration (DSA) mechanism dBlandford and Ostrikerlll978l 
Bell 1978a. b). With respect to massive stars, the latter process has been invoked by 
Whitd (119851 . for important modifications to the original idea see lChen and Whita 
19941) to explain the NT emission in suspected single stars, where the required 
shocks should be th ose embedded in t he win d and created by the line-driven in- 
stability (Sect. 16.2b . lEichler and Usovl dl993h . on the other side, showed that also 
the shocks created by stellar wind collisions in early type binaries are potential 
sources of the observed N T emission, in line with even earlier suggestions by 
IWilliams et al dl98lll990l) regarding the non-thermal emission from the WC7-I-05 
binary WR140. More recently, jPittard et al (200l) and iReimer et all d2006h have 
developed models applying the DSA mechanism to the corresponding physical 
conditions of hydrodynmical shocks in colliding wind binaries. 

The presently favorized hypothesis is that non-thermal emission from massive 
stars is generally procuded by the latter mechanism. Two elements have simulta- 
neously led to this idea. On the observational side, a number of systematic ob- 
servational studies of NT emitters to investigate their multiplicity have revealed 
that the vast ma jority of them are confirmed (or at least suspected) binaries (see 
iDe Beckeiil2007) . Note, e.g., that already Dougherty and Williams, (.2000,) found 
strong indications that at least for WR-stars a massive compa nion is a prerequisite 
for the appearance of NT emission. On the theoretical side, v an Loo et a il (I2OO6L 
based on Van Loo 2005) showed that the observed spectral shape of non- thermal 
emission cannot be reproduced by current (1-D) hydro-simulations applying the 
wind-instability scenario in presumed single 0-stars. They suspected that also all 
0-stars with non-thermal radio emission should be members of binary or multi- 
ple systems, where the non-thermal component is created in the shocks associated 
with colliding stellar winds. 

In the present context, NT emission is particularly important because of its im- 
pact on the radio spectrum and flux. If present, it can strongly perturb mass-loss 
rate determinations based on the thermal radio flux, as, e.g., pointed out bv lSteveni 
([1995). In this context and assuming the above hypothesis of NT emission cre- 
ated in colliding winds, high angular resolution observations allowing for a direct 
disentangling of purely thermal (single star winds) and non-thermal (colliding- 
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winds) emission in wide binaries will prove as highly valuable. However, let us 
also note that an additional thermal radio cont ribution may be produced by th e 
shocked gas in the colliding-wind region itself ( Stevenslll99i iPittard et alll200^ . 
therefore contributing to some additional confusion. 

4.1.2 Methods based on NLTE atmospheres 

As outlined above, more recent wind diagnostics is mostly based on NLTE model 
atmospheres. At least three different methods can be envisaged: 

(i) Synthetic spectra from self-consistent atmospheric models (comprising a NLTE 
+ hydro-description, where the wind-structure relies on a line force resulting from 
self-consistent NLTE occupation numbers and the corresponding radiation field, 
see, e.g.. Fig. 1 in Pauldrach et al 1994) are fitted to observations which cover a 
significant spectral range and are sensitive to both photospheric and wind condi- 
tions. The only input (= fit) parameters are Tgff, \ogg, /?*(or different combinations, 
e.g., Tgff, L*, M*), the individual elemental abundances, and, if necessary, (strati- 
fied) clumping factors and a description of the X-ray/EUV emission. Besides the 
fact that such an approach is very time-consuming, it relies on the assumption that 
the complete underlying wind-physics is correctly treated. This method does not 
allow for an independent test of the theory itself. All inaccuracies/problems of the 
theory will, of course, show up in deviations between synthetic spectra and ob- 
servations. Unfortunately, however, these will be interpreted in terms of peculiar 
abundances, clumping factors, radii etc., since for a given set of stellar parameters 
only one specific wind-structure (i.e., a specific combination of M, Voo and velocity 
field) is possible, which may not be correct. The better the theory, the better the 
results from such an approach, of course. An interesting application of method (i) 
will be presented in the context of determining stellar parameters of Central Stars 
of Planetary Nebulae via wind diagnostics. 

(ii) Synthetic spectra from consistent atmospheric models (NLTE + hydro-descrip- 
tion) are fitted to observations. In contrast to method (i), M and Voo are adapted to 
fit the wind-lines, but not directly (as in method (iii) discussed below), but by 
varyin g th e force-multiplier parameters k, a, 5 (or equivalent quantities) in a rea- 
sonablqi_| way. The fit parameters are as above, augmented by the force-multiplier 
parameters. The advantage of this method relies in the possibility to obtain a phys- 
ically justified stratification of the velocity field (and the indirect possibility to 
check the wind-theory, by checking whether the observations allow for reason- 
able force-multiplier parameters), but this advantage can also become a disad- 
vantage when the actual velocity field cannot be matched by the hydrodynamic 
approach, due to shortcomings in the description. Method (ii) is usually applied in 
investigations using the atmosphere/wind-code WM-Basic (see Sect. 14.2.11) . and 
impressive fits of the complete UV region have been obtained f or a (small) num- 
ber of s tars, e.g., C P up, a Cam (09.5Ia) dPauldrach etalll200ll) and HD93129A 
r02If*. rTaresch et a L ,1997) . in the latter case including the optical region. Simi- 
lar analyses o f a larger sample of Galac tic O supergiants and dwarfs h ave been 
performed by ' Bianchi and Garcial (l2002h and lGarcia and Bianchil (|2004|) (see be- 
low). Note that in almost all these cases the inclusion of the observed X-ray/EUV 



within the range predicted by line-statistics, see pagesll Iff. 
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emission (interpreted as due to the emission from cooling zones of shocks em- 
bedded in the wind, see Sect. 16.11) has been cr ucial to explain the ionizatio n of 
highly ionized species such as O VI and S VI, see lPauldrach et all(ll994ll200ll) and 
[Macfarlane et al ( 1993). These ions which show strong P Cygni lines in the FUV 
(e.g.. Fig. [Ill cannot be created in "standard" NLTE calculations based on a wind 
with an electron temperature of the order of Tgff, a problem realized already at 
the beginrung of ho t star UV spe ctrosc opy and d enoted as superionization (e.g., 
iLamers and Mortonlll976llCastoiill979llHamann|[l980h . 

(iii) The last method enumerated here is the most frequently used. Synthetic spec- 
tra from a NLTE model are fitted to the observations, where the wind-structure is 
analytically described via a j3 velocity law (Eq. [38l) , and a smooth transition be- 
tween the analytical wind structure and a quasi-hydrostatic photosphere is adopted. 
Input-parameters (in addition to those from method (i)) are M, and j3. Such 
models (and those calcula ted by method (ii)) are called Unified Model Atmo- 
spheres dGabler et aJ 1989 1). and are standard nowadays. Corresponding codes will 
be discussed in the next section. 

Let us finally note that not only different methods can be (and are!) combined 
with (different) atmospheric codes, but that also the fit-optimization method can 
be different. Even in recent investigations, a "fit-by-eye" method has been fre- 
quently used, but more objective and automatic methods such as the optimization 
by genetic algorithms (Mokiem et al.2005i) or usi ng predefined grids in combina- 
tion with a y^/maxin ium likelihood optini ization (lLefevej2007h or via a principal 
component analysis dUrbaneja et alll2008l) have also entered the spectral analysis 
of hot stars with winds. 



NIR spectroscopy. The dramatic progress of IR astronomy in the last two decades 
has opened a completely new window for the systematic investigation of hot star 
winds, particularly for the analysis of dust-enshrouded, very young objects and 
samples in highly reddened clusters. Note also that the current generation of 10m 
class ground-based telescopes has been optimized to observe in this wavelength 
range. 

Also the IR spectra from hot stars show pure emission lines, wind-contami- 
nated absorption lines and sometimes even P Cygni profiles, and these lines have 
to be analyzed by means of NLTE models due to strong NLTE effects in certain 
transitions. Since the IR continuum of objects with strong winds is formed already 
in the wind, however, IR lines may sample different depths inside the wind and 
provide additional information about the shape of the velocity field and particu- 
larly on the clumping properties. 

The winds of hot stars of extreme luminosity and with strong IR emission 
characteristics in the Galactic Center have been investigated by different groups 
(e.g., i Naiarro et al 1997b, Najarro 2006, Martins etal 2007), and Lenorzer et aj 
ilOO^ . iRepolust et all ((20051) and lGeballe et all (l2006h have demonstrated in how 
far stellar and wind parameters of "normal" OB-stars may be constrained by IR- 
spectroscopy alone. One of the most important applications of IR line diagnostics, 
however, will be the measurement of mass-loss rates from stars with very weak 
winds by means of Bra, as detailed in Sect.jSj 
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4.2 Recent results 
4.2.1 OBA-stars 

Most of the relevant results in the field of quantitative spectroscopy of massive 
hot stars that have emerged with the new millennium are the outcome of more 
than ten years of considerable efforts made by different research groups to incor- 
porate, in a realistic way, the effects of line-blanketing in NLTE-modeling of the 
atmospheres of massive stars. Indeed, the possibility to include the role of several 
thousands (up to millions) of metal lines was intimately c oupled to the revolu- 
tion of computer facilities. First approaches were made by 'Abbott and Hummed 
(fl 985) who modeled the effect of "wind blanketing" by introducing a wavelength 
dependent albedo representing the radiation reflected back onto the photosphere 
by electron and line scattering from the wind. Abbott and Hummer showed that 
the main result of "wind blanketing" was to heat the photosphere from the sur- 
face to the depth of continuum formation ("back-warming"), so that the resulting 
spectrum would resemble that of a hotter star. Thus, depending on the strength of 
the stellar wind (M), the albedo could lead to a net reduction of the derived stellar 
Tgff by up to 20%. The main caveat of this approach was the neglect of photo- 
spheric "line blanketing", so that a unified picture could not be provided. As well, 
effort has b een undertaken to incorporate the effe cts of photospheric "flux block- 
ing" alone (lHerreroetalll2000l . IVillamarizll200lh . However, since only blocking 
was considered but no back-warming, revised Teff values were significantly higher 

than models without blocking. 

One of t he major steps in the blanketing affair was presented by lAndersoiil 
(I1985L 119891) who introduced the so-called "superlevels". Instead of dealing with 
all individual energy levels of metal species, several states with close enough ener- 
gies were packed into one superlevel. Thus, the number of statistical equilibrium 
equations could b e drastically reduced saving an enormous amount of computing 
time (but see, e.g., Lucv 2001 on potential problems inherent to this approach). Af- 
ter some pioneering work (Hauschildt 1992), by the end of the nineties, fou r major 
codes, PHOENIX iH auschildt 1992), CMFGEN (Hillier and Milleilll998h . WM- 



Basic (iPauldrach etal.19 98, 2001) and TLUSTY (,Hubenv.l998i plane-parallel, no 
wind) could address line blanket ing in a rea listic way. Just a c ouple of years later 
two more codes PoW R dOrafener et aJ2002l ') and FASTWIND dSantolaya-Rey et all 
ll997L lPuis et alll2005l) becam e available. For a brief comparison of the features of 
these codes, see lPulsl ( l2008l) . 



Effects of line blanketing. To illustrate the two main direct effects of line-blan- 
keting, we have made use of CMFGEN and computed a model for a late B- 
hypergiant with a strong wind. Fig.|9]shows the comparison of the blanketed/non- 
blanketed cases for the emergent flux distribution and temperature structure of 
the atmosphere. The models confirm the effects anticipated from early on and 
discussed above. We clearly see how the blocked flux in the UV emerges at opti- 
cal and longer wavelengths. As the flux is blocked and a certain amount of pho- 
tons are backscattered, flux (luminosity) conservation demands a more efficient 
photon diffusion in the inner photosphere. Since this effect is controlled by the 
temperature gradient, the temperature is increased in this region (back-warming) 



44 



J. Puis, J.S. Vink & F. Najarro 




3.0 3.2 3.4 3.6 3.8 4.0 4.2 4.4 10 11 12 13 

Log \ A Log N_ (cm-") 



Fig. 9 Effects of line blanketing (solid) vs. unblanketed models (dashed) on the flux distribution 
(logFv (Jansky) vs. logA (A), left panel) and temperature structure (7(10'* K) vs. log iie, right 
panel) in the atmosphere of a late B-hypergiant. Blanketing blocks flux in the UV, redistributes 
it towards longer wavelengths and causes back-warming. 
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Fig. 10 Effects of line blanketing on the hydrogen ionization structure (log H I /H ll vs. log^e, 
left panel) and line profiles (H„, right panel) of a mid B-supergiant. The enhanced ionization 
caused by blanketing (solid vs. dashed) weakens considerably the Ha -profile (dashed-dotted) 
with respect to the unblanketed model (dashed) and requires an increase of M by up to a factor 
of three in B-supergiants to balance the effects of blanketing and to fit the observations (solid). 



as shown in Fig. |9] Back-warming, in combination with the enhanced radiation 
field due to backscattering, will lead to enhanced ionization in the inner parts 
of the atmosphere. The question arising now is how the spectral lines will be 
affected and whether the effect of blanketing could be considered as a general 
shift in the ionizing condi tions throughout the whole atmosphere (as found by 
I Abbott and Hummenll985h or whether it will affect differently the inner photo- 
sphere and the outer wind. Thus, we want to know whether this enhanced ioniza- 
tion will be kept throughout the atmosphere or whether the stellar wind will block 
enough ionizing radiation and decrease the ionization degree in the outer wind 
parts. The outcome of this competition will depend on the parameter domain in 
which the star is located, and will be very sensitive to stellar properties such as 
effective temperature, gravity and wind strength. For an illustration, we will com- 
pare the situation for a mid B-supergiant and a very hot early 0-supergiant. 

Fig. 1 101 shows the effects of line blanketing on the hydrogen ionization struc- 
ture of the mid B-supergiant. Despite the strong wind, the enhanced ionization 



Mass loss from hot massive stars 



45 




ior , ; • 

2.0 2,5 3,0 3,5 4,0 

Log X (A) 



Fig. 11 Effects of line blanketing on the helium (left) and hydrogen (middle) ionization struc- 
ture of an 03 If supergiant. Solid: Tgff= 45.5 kK, blanketed; dashed: Tgff= 45.5 kK, not blan- 
keted; dashed-dotted: reff= 50 KK, not blanketed. See text. Right: Flux distribution of models 
with (dashed) and without (solid) blanketing matching the stellar line spectra and photometry. 

caused in the inner photosphere is held throughout the wind. Further, the increased 
depletion of H I in the atmosphere through blanketing requires an enhancement 
of M by up to a factor of three to compensate the extra ionization and match 
the observed profile (see Fig. [TOll. Thus, for mid B-supergiants, blanketing af- 
fects majorly the derived mass-loss rate, while milder changes appear on the de - 
rived effective tempe rature. Similar trends have been found by lUrbanejal (|2004|) . 
ICrowther et ail (l2006h and lMarkova and Pulsl ( l2008l) . 

Fig. [TT] displays the effects of line blanketing on the He and H ionization 
structure of an 03 If star. For the case presented in Fig. (TT] we have selected 
three different models. The first one corresponds to a blanketed model with Tf.ff 
= 45.5 kK, logg = 3.7 and R = 15.6 Rq. The second model has the same pa- 
rameters but no blanketing, while the third one, again not blanketed, is hotter, 
Tgff = 50 kK (log^ = 3.8 and R =14.6 Rq), and reproduces equally well the ob- 
served spectra of the O supergiant and the observed V magnitude. All models 
have the same wind-strength parameter Q. From Fig. [TTJleft we see that blanket- 
ing (solid) severely enhances the He ionization in the line forming region with 
respect to the cool unblanketed model (dashed), so that a degree of He ionization 
similar to the hot unblanketed model (dashed-dotted) is obtained by decreasing Tf-g 
by almost 5000 K. In other words, the basic diagnostic for temperature determi- 
nations in O stars, namely the He ionization balance, reacts strongly to blanketing 
and causes the observed reduct ion on the effective temperature scale (for further 
details, see lRepolust et alll2004l) . 

Inspection of the impact of blanketing on the H ionization structure of O stars 
(Fig. [TTl- middle) tells us that only minor changes are produced in the ionization 
degree of H and, therefore, only small variations should be expected in the in- 
ferred mass loss rates (assuming they are derived from H^, as Hell 4686 would 
be severely altered). Finally, it should be noted that the relatively strong wind of 
the star blocks the He II ionizing radiation and forces recombination of He III to 
He II, with important implications on the number of He II ionizing photons (see 
also FiglTllright). Thus, this is one case where the strong wind controls the ion- 
ization and outweighs the enhancement produced through backwarming in the 
inner regions. 

First impact: The new temperature scale for massive OB stars. IVacca etall(ll996t) 
presented a compilation of the spectroscopic determinations of effective temper- 
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Fig. 12 The tem perature sca l e for G alactic 0-dwarfs (left) and supergiants (middle). The solid 
lin es are for the IVacca et al' (1996^ scale, the dashed li ne for the "th eoretical" scale defined 
by iMartins et all 12005a), filled symbols are data from iRepolust et all ^2004) and open sym- 
bols are data from Herrero et al (2002). Right: The new temperature scale for B-supergiants 
from Markova and Puis ( 2008), combining results from Urbaneja ( 2004), Crowther et al ( 2006), 
l^zvbilla et all d 2006) , iLef ever et all I I2007I) and their own data. See text. 



atures of massive OB stars. They gave preference to the most recent calculations 
which at that time were mostly based on plane-paralle l, hydrostatic, unb lanketed 
NLTE model atmospheres. Interestingly, in their paper. [Wcca et all d 1996 ) warned 
about the significantly lower temperatures obtained by the few available wind- 
blanketed analyses when compared to unblanketed results. 

Once the codes cited above became available, a series of papers were published 
mostly devoted to the issue of temperature scales. The fi rst "modern" calcu lations 
pointing to a cooler temperature scale were those from lMartins et all (120021) using 
CMFGEN. These authors limited their calculations to OB-dwarfs, so that the in- 
fluence of m ass-loss effects were negligible. Therefore, the main differences with 
IVacca etai! were clearly due to line-blanketing. These differences could reach up 
to 4000 K for early types, and decreased towards 09 and BO types, as can be 
ap preciated in Fig. [121 

iHerrero et all (l2002h provided a temperature scale for supergiants in Cyg 0B2 
using FASTWIND. They found differences of up to 8 000 K. In this case both mass 
loss and line -blanketing played a role. These authors also showed that two stars 
with the same spectral type and luminosity class may have different effective tem- 
peratures if their wind densities are different. While the results from Herrero et aj 
were based on the analyses of only seven Cyg 0B2 supergiants, iRepolust et al 
(12004 ) presented an analysis of 24 stars (17 supergiants and giants and 7 dwarfs), 
based on a slightly improved version of FASTWIND that confirmed the same 
trends. 

(Almost) identical trends have been found by iMartins et all (12005 al) who pro- 
videj^ by means of CMFGEN models, new temperature scales for massive OB 
stars of differe nt luminosity classes. These scales agree quite well with those from 
IRepolust et all and confirm that new models result in effective temperatures that 
are several thousands Kelvin cooler for early and intermediate spectral types, de- 
creasing towards late spectral types. A comparison of all these different tempera- 
ture scales is presented in Fig. [12| 

Likewise detailed studies have been performed for B-supergiants dUrbanejai 

l2004 |Crowther et al 2006, Lefever et al 2007, Mark ova and Pulsl2008l) which show 



among other scaling relations for atmospheric parameters and ionizing photon numbers. 
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that also the temperature scale of Galactic B-supergiants needs to b e revised down- 
wards . Compared to the unblanketed and "wind-free" scale by iMcErlean et all 
(Il999h . this amounts from 10 up to 20%, the latter value being appropriate for 
stronger winds. The resulting new Teg-scale (described by a third order polynomial 
fit, with gray-shaded l-a uncertainties) is displayed i n Fig. | 12f right , and agrees 
perfectly with the results from an independent study bv lTrundle et all (l2007h . 

Similar investigations have been performed for OB-stars in the Magellaruc 
Clouds (MCs), in particular to study metallicity effects on the effective tempera- 
tures and mass-loss rates. Massev et al (2i)04, 2005) investigated a large sample of 
MC 0-stars by means of FASTWIND, and provided a spectral-type-Teff calibra- 
tion for the SMC. For the LMC, the situation remained unclear, since their sample 
was concentrated towards the hottest objects, 02-04. Overall, it turned out that 
for a given spectral sub-type, TeffCSMC dw) > 7;ff(MW dw) w TeffCSMC sg) > 
7eff(MW sg ), where the Teff-s cale for SMC 0-stars differs much less from the 
unblanketed I Vacca et all d 19961) calibration than the scale for their Galactic coun- 
terparts ('dw' = dwarfs, 'sg' = supergiants). This finding has been attributed to less 
blanketing and weaker winds because of the lower SMC metallicity (see below). 
For dwarfs, most of the results derived by Massey et al. are in go od agreement 
with investigations of different samples performed with CMFGEN dMartins et all 
[2004) and TLUSTM3(Bouret et al 2003) , whereas a large number of MC super- 
giants analyzed by means of CMFGEN dCrowther et ajl200l iHillier et allbool 
lEvans et all2004al) turned out to be significantly cooler, even cooler than implied 
by the Galactic scale. In certain cases, this might be explained by wind-blanketing 
effects, since some of the discrep ant objects are r ather extreme, but other cases 
certainly need further inspection. , Heap et all d2006h analyzed a sample of SMC O- 
dwarfs and giants, and compared their results (using TLUSTY) with other stud- 
ies. Though they stress the large scatter of Tgff within individual sub-types, they 
found fair agreement with othe r investigations, except for Galactic data derived by 
iBianchi and Garcial d2002h and lGarcia and Bianchil d2004h by means of a pure UV 
analysis (using WM-Basic), which seem to suggest systematically cooler temper- 
atures than all other studies. Also this discrepancy has to be investigated in the 
near future. 

The Vir-FLAMES survey of massive stars. One of the most important recent 
projects on OB-stars was the VLT-FLAMES survey of massive stars (FLAMES = 
Fibre Large Array Multi-Element Spectrograph). By means of this campaign, the 
massive stellar content of 8 young and old Galactic/MC clusters has been spectro- 
scopically investigated, in order to answer urgent questions regarding (i) rotation 
and abundances (rotational mixing), (ii) stellar mass loss as a function of metal- 
licity and (iii) fraction and impact of binarity. In total, 86 0-stars and 615 B-stars 
have been observed at high resolution. An overview of objects and objectives has 
been given by Evans et al (2005, 2006), and a summary of important results can 
be found iniEva ns et al (2008). 

iMokiem et'al (.2006 . ,2007 a'FI studied the 0-/early B-star targets of the FLAMES 
survey in the SMC and LM C, respectively. They confirmed the basic results from 
iMassev et all d2004l |2005|) . but refined the spectral-type- rgff scale, particularly 



augmented by a CMFGEN analysis of the wind parameters. 

by means of a genetic algorithm combined with FASTWIND, see Sect. l4.1.2l 
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with respect to the LMC objects. They showed that, at least for 0-dwarfs (which 
are not "contaminated" by additional wind-effects), the effective temperatures for 
a given spectral sub-type decrease with increasing metallicity, i.e., rejf(SMC) > 
TeffCLMC) > TeffCMW). A similar result was derived by Trun dle et al (2007) for 
the FLAMES B-type dwarfs, in this case based on TLUSTY model atmospheres. 



Wind properties of OB stars at different metallicities. Wind-properties of Galac- 
tic/MC OB-stars (primarily mass-loss rates, M, and velocity field parameters, j3) 
have been determined by numerous investigations (overlapping with those men- 
tioned above and detailed in Puis 2008), from Hq (partly combined with He II 
4686) and the UV. In most cases, terminal velocities, v^ c, have been adopted from 
UV-measurements (see Sect. 14.1. it and/or calibrations dKudritzki and Pulsll200(]|) . 
Wind-momentum luminosity relations ( WLR) have been inferre d and compared 
with theoretical predictions, mostly from I Vink et all (l2000l . l200lh . see Sect. [3] In 
summary, the following results have been obtained: 

(i) For most O-Zearly B-stars, the theoretical WLR is met. Notable exceptions are 
0-supergiants with rat her dense winds, where the "observed" w ind-momenta ap- 
pear as "too large" (seel Markova et alll2004l . lRepolust et alll2004) . which has been 
attributed to wind-clumping, see Sect. 16.31 and low luminosity 0-dwarfs, where 
the "observed" wind-momenta are considerably lower than predicted, which has 
been denoted as the "weak wind problem" and will be covered in Sect.|5] 

(ii) The bi-stability jump in the B-s tar domain is repre s ented by a gradual de crease 



in Voo over the bi-stability region dEvans et alll2004bl ICrowther et alll2006l). bv a 



factor of roughly 2.5 with respect to the ratio of Voo/vesc- According to lMarkova and Pula 
(|2008), the limits of this region are l ocated at somewh at lower Teff (betweeen 
18 to 23 kK) than those discussed by Vink et al (2000). Most important, how- 
ever, is the finding that M changes over the bi-stability region by a factor (in be- 
tween 0.4 and 2.5, iMarkova and Pulsll2008h which is smaller than the predicted 
factor of 5. Thus it seems that the decrease in Voo over the bi-stability region is 
not over-compensated by an increase of This finding seems to apply not 
only to the bi-stability region itself, but to the complete "low temperature region" 
(i.e., all mid/late B-type stars), where the predicted mass-loss rat es are larger than 



those found from NLTE analyses based on unified atmospheres dVink et all200(]L 



Trundle and Lennon 2005, C rowther et alll200^lLefever et alll2007llBenaglia et al 



2007, Markova and Puis 2008h . This problem is reminiscent of the findings by 
Kudritzki et al ( 19991) who derived wind-momentum rates of mid B-supergiants 
which were much smaller (by a factor of 10) than those of early B-supergiants. 
Bearing in mind that Kudritzki et al (1999) used unblanketed atmospheres, even 
if applying an upwards "correction" of their mid B-supergiant mass-loss rates by 
a factor of three (accounting for the arguments from pagel45]), their results would 
still be in contrast to the predictions. 

Nevertheless, a certain effect dividing hotter from cooler winds is present in- 
deed, evident from the decrease in the ratio Voo/vgsc and likely related to t he princi- 
pal bi-stability mechanis m. Moreover, as shown bv lBenaglia et all d2007l see also 
iMarkova and Pulsll2008l) . a local increase of the ratio Mlv^ seems to be present 

^' at least not if the winds from hotter objects are not substantially stronger clumped than those 
from the cooler ones. 
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Fig. 13 Comparison of observed (solid, with grey shaded l-a confidence intervals) and theo- 
retically predicted (dashed, Vink et al 2000, 2001) WLRs for Galactic, LMC and SMC 0-/early 
B-stars (from top to bottom). Adapted from .Mokiem et al.2007h . see text. 

at least inside the transition zone, which would partly support the predictions by 
Vink et al., though not regarding the mass-loss properties of objects below the 
jump. Further effort to clarify the situation is certainly needed. 

(iii) Combining the results of the various investiga tions of Galactic OB-stars and 
A-supergiants (the latter from lKudritzki et aHl 19991) . the corresponding WLRs (as 
a function of spectral type) extend over significant ranges in Tgff, assuming that 
clumping is not a strongly varying function of spectral type. Excluding weak 
winds, three rather well-defined relations for O and early B-stars (roughly above 
23 kK), for mid and late B-supergiants (between roughly 18 kK and 10 kK) and 
for A-supergiants seem to exist, where the slope of the former two relations is 
rather similar, and the slope of the latter is steeper , in accordance with theoret- 
ical predictions. In particular. iKudritzki et all d 19991) derived a slope correspond- 
ing to a' = 0.38 ± 0.07 for their sample of A-supergiants, which is in perfect 
agreement with corresponding results from line-statistics (page [13]) . Within the 
(re-defined) bi-stability region (18 to 23 kK), on the other hand, the scatter seems 
to be mi ich larger, possibly due t o ongoing changes in the global ionization equi- 
librium dMarkova and Pulsll2008h . 

Translated into the number of driving lines (or force-multiplier parameter k, 
cf. Eg. list, the above results imply that this number remains rather co nstant over 
large ranges in T^ff (again in accordance with theoretical predictions, IVink et all 
I2OOOI and Sect. 12.11) . which are divided by the transition from Fe IV to Fe III and 
from Fe III to Fe II as the major driving agents. The "only" difference between 
these principal predictions and the observations is a quantitative one, regarding 
the location of the (first) bi-stability region and the degree to which extent A/eff (or 
k) changes from one region to the other. 

(iv) For a given luminosity, the mass-loss rates of SMC-stars are lower than for 
their Galactic counterparts, consistent with theory dMassev et all 12004. 2005|). A 
more precise quantification of the metallicity dependence of the winds from O- 
/early B-stars was possible within the FLAMES survey. From the analysis of the 
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Fig. 14 WLR for Galactic and extragalactic A supergiants. Dashed: Linear regression for Galac- 
tic and M3 1 objects. Dotted: Galactic relati on scaled to the mean abundan ce found for NGC 300 
and NGC 3621, Z/Z; = 0.4. Adapted from lBresoUn and KudritzkHllOOl . 



SMC/LMC objects by iMokiem et all (l2006l) and iMokiem et al d^OTa), respec- 
tively, and in combination with data from previous investigations, Mokiem et a| 
(l2007b ) derived the WLRs for Galactic, LMC and SMC objects, with rather nar- 
row l-<7 confidence intervals, and showed that the wind-momenta strictly increase 
with metallicity Z (i.e., the WLR of the LMC lies in between the corresponding re- 
lations for the MW and the SMC, see Fig.Qill. Using Z(LMC) w O.SZg, Z(SMC) 
w 0.2Z©, and allowing fo r a modest "clumping correcti on" for dense winds (fol- 
lowing the arguments by iRepolust et all 12004 see also iMarkova et alll2004] and 
Sect. 16.31 1, they obtained M oc 20.72±o.i5 -pj^j^ result is in very good agreement 
with theoretical predicti ons, both from (statio nary) models (Sect. 13.4.1b and from 
line statistics (page fTTl) . Mokiem et all (l2007bl) pointed out that the derived metal- 
licity dependence should remain unaffected from potential future changes in the 
absolute values of mass-loss rates due to global clumping corrections, as long as 
the clumping properties were similar in winds of different metallicity. 

Beyond the Magellanic Clouds. As already noted in our discussion of the (theo- 
retical) concept of the WLR (Eq.[T6l), one of the early aims to use this relationship 
was to derive extragalactic distances by purely spectroscopic means. Particularly 
advantageous for this purpose are early A-supergiants, because they are the visu- 
ally brightest "normal" stellar objects. Also, problems of multiplicity and crowd- 
ing are (almost) negligible because of the ir brightness, their sho rt life-times, and 
their relatively old evolutionary ages (see lKudritzki et alll2008al) . Moreover, their 
Hfy-li i ie can be used to measure M and v«, simultaneously (e.g., ' Kudritzki et all 
1 19991 iKudritzki and Puis 2000). Subsequently, a number of these objects have 
been analyzed in Local Group galaxies and beyond to compare with Galactic 
counterparts and to check the theoretical concept. A summary of the results has 
been given by Bresolin and Kudritzki (2004), highlighted in Fig. [14] which dis- 
plays the corresponding WLRs. The relation for "solar" metallicity is defined by 
Galactic objects plus two M31 stars, and additional wind-momentum rates for 6 
A-supergiants in NGC 300 and one object in NGC 3621 are displayed. Assuming 
that the latter stars have an abundance similar to the mean found for their host 
galaxies, Z w O.4Z0, the dotted line provides the theoretical expectation when the 
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Galactic WLR is scaled to the lower metallicity. It is evident that the WLR concept 
seems to hold. Details regarding extragalactic samples, observations and analyses 
can be found in lBresolin and Kudritz kil (120041) and references therein. 

After these promising investigations had been finished, however, no further 
atte mpt was made to de rive distances via the WLR. The reason for that was 
that lKudritzki et all (12003 1 had developed an alternative possibility to use late B- 
/early A-supergiants as standard candles, by means of what they called the Mux- 
weighted Gravity Luminosity /Relationship. The authors showed that the flux- 
weighted gravity, g/Tf.f{^, and the absolute bolometric magnitude are strongly 
correlated for these objects, since they pass their short evolutionary phase with 
nearly constant luminosity and mass, and because luminosity and mass are re- 
lated via L oc with jc w 3. Thus, the underlying difference in both methods is 
that the WLR method exploits the information on radius contained in the stellar 
wind, whereas the FGLR method relies on evolutionary facts. 

As shown in the original paper and a follow-up study ( Ku dritzki et alll2008bl) . 
this new method to derive distances seems to have advantages over the WLR 
method, because of the high precision that is possible even for low -re solution ob- 
servations (0.1 mag in the distance module when analyzing 10 target stars) and its 
simplicity. E.g., the FGLR method requires much less (and only photospheric0) 
parameters to be determined (Teg, log^, and metallicity), and is also much less 
affected by additional dependencies. Remember, e.g., that the mass term in Eg. [141 
only partly vanishes for A-supergiants, since a' ^ 0.4 and not close to 2/3 as in 
the OB-star case. Moreover, it is still unclear whether and in how far clumping 
affects the WLR method. 

A fir st impressive demonst ration of the power of the FGLR method has been 
given by lUrbaneja et all (l2008l) who determined the distance to the Local Group 
Galaxy WLM with a distance modulus of 24.99±0.10 mag (which is only 0.07 mag 
larger than the most recent Cepheid distance). Further applications within the 
Arauc aria project ("Measuring improved distances to nearby galaxies". iGieren et all 
l2005h are to be expected soon. 

Central Stars of Planetary Nebulae. Although Central Stars of Planetary Neb- 
ulae (CSPN) are objects of significantly lower luminosity and of completely dif- 
ferent evolutionary status, they have winds with a mean density similar to massive 
stars, and their wind-momenta seem to correspond to an extrapolation of the WLR 
of 0-stars dKudritzki et al|[l997l I2006L iHultzsch et af2007h . though with a larger 
scatter. This finding must be regarded as an encouraging success of the interpreta- 
tion of winds in terms of radiative driving and of the concept of the WLR. Let us 
point out that recent analyses by means of FASTWIND have suggested that some 
(but not all) of these winds are strongly clumped, which could be quantified be- 
cause of the different reaction of the He II 4686 and the H„ line on wind-clumping 
at temperatures around Teg w 30 kK ("Kudritzki et alll2006h . 

In an interesting study, Pauldrach et al (2004i) analyzed the winds of 9 CSPN 
by means of WM-Basic and method (i) as outlined in Sect. 14. 1.2] i.e., they derived 
the complete set of stellar parameters by comparing self-consistent wind models 
with UV-observations. The information provided by the wind features (Voo, M) 



at least when discarding extreme objects with wind-affected SEDs. 
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permits to derive the stellar radius via such models, since the corresponding de- 
pendencies Voo = /(vesc) and M = f{Lit,Mif,r) allow for a unique solution for /?* 
giv en the observat ional constraints (for a discussion of the underlying principles, 
see lKudritzki et a]| [l992). By this method, all quantities such as luminosity, mass, 
and distance can be constrained without relying on the (core)mass-luminosity re- 
lation tak en from stellar structu re theory and post-AGB evolutional Notably, the 
results bv lPauldrach et all (|2004) showed severe departures from this generally ac- 
cepted relation, and masses from 0.4 to 1 .4 Mq have been derived, with five out 
of the nine sta rs located close to (but not above) the Chandrasekhar mass limit for 
white dwarfs. IPauldrachI (l2005h argues that these objects might belong "to a not 
yet understood subgroup of CSPNs that evolve to white dwarfs which can end up 
as supernovae of type la" (see also Maoz 2008). For further progress on this issue, 
it needs to be clarified in how far the wind and spectrum synthesis calculations 
are not biased by contaminating effects. An important step towards this is the de- 
termination of precise distances to CSPN (to enable a measurement of the stellar 
radius) for which fairly resolved spectra can be obtained as well. 

4.2.2 Luminous Blue Variables 

What is an LBV? Luminous Blue Variables represent a short-lived (^ 10^ — 10^ 
years) phase of massive star evolution in which the stars are subject to significant 
effective temperature changes. They come in two flavors. The largest population 
of ^30 LBVs in the Galaxy and the Magellanic Clouds is that of the S Doradus 
variables with magnitude changes of 1-2 magnitudes on timescales of years to 
decades. These are the characteristic S Dor variations, which are represented by 
the dotted horizontal lines in Fig. [l5j The general understanding is that the S Dor 
cycles occur at approximately constant bolometric luminosity (which has yet to 
be proven) - principally representing temperature variations. The second type of 
LBV instability involves objects that show truly giant eruptions with magnitude 
changes of order 3 — 5 during which the bolometric luminosity most certainly 
increases. In the Milky Way it is only the cases of P Cygni and Eta Carina which 
have been witnessed to experience such extreme behaviour. 

Whether these two types of variability occur in similar or distinct objects is 
not yet clear, however in view of the "unifying" properties of the LBV P Cygni it 
is highly probable that the S Dor variables and giant emptors are related, that they 
are in a similar evolutionary state, and th at they are subject to the same type of 
instabilities near the Eddington limit (see IVinl3l2009l for a review). We note that 
the true nature of the LBV instability has yet to be revealed. 

LBV parameters and abundances. Reliable determinations of stellar properties in 
LBVs became possible only after the advent of unified model atmospheres. Even 
then, however, in some extreme o bjects like r\ Car, the e normous complexity of 
the observed ground based spectra (iHillier and AlleiJ 19921) . contaminated by mul- 
tiple emission regions, hinders a proper quantitative spectroscopic analysis. In this 
case, stellar parameters need to be determine d either from "c leaner" wavelength 
regions (e.g., from the millimeter continuum, ICox et alll 19951) or from "cleaner". 



which is usually done when purely spectroscopic analyses such as outlined above are per- 
formed. 
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Fig. 15 The position of the LBVs in the upper Hertzsprung-Russell diagram. Note that the 
general population of blue supergiants that occupy the same part of the HRD are not plotted. The 
slanted band running from 30 kK at high L/L,.;, to 15 kK at lower luminosity is often referred 
to as the S Dor instability strip. The vertical band at a temperature of ~ 8000 K represents the 
position of the LBVs in outburst. The vertical line at 21 000 K is the positio n of the observed b i- 
stability jump (the center of the observed bi-stability region). Adapted from lSmith et all i2004D . 



resolved observations (e.g., using the HST. iHillier et alll200lL or the recent, re- 
ally promising res ults from long baseline interferometers such as amber-vlti, 
IWeigelt et afcoOTl) . 

The outstanding wind density (with lO~^M0yr~^as estimated for this ob- 
ject, |CoxeLalJ^95|,|Hnnereta||200^ places /?(tross = 2/3) at 80% of the ter- 
minal velocity, impeding any derivation of the hydrostatic radius. Nevertheless, 
CNO abundances could be determined and are consistent with those found for the 
surrounding nebula (for a recent review, see Naiarro and Hillier 2 008) . 

For less extreme LBVs such as P Cygni (M^ 10~^Mf7^vr~^ iNajarrolllOOlh . 
more detailed investigations could be carried out. A major breakthrough was ob- 
tained by Langer et al (1994) who presented a quantitative study o f this object, 
combining evolutionary and hydrodynamical aspects. Subsequentlv. lNajarro et all 
fl997 a) were able to accurately constrain the main stellar properties, including 
Helium abundance, by means of homogeneous, unblanketed atmospheric models 
and UV to radio observations. These result s were only m inorly altered by utiliz- 
ing blanketed models including clumping (lNaiarroll200ll) . which confirmed the 
presence of CNO processing and provided, for the first time, direct metallicity 
measurements for an LBV. 

Apart from the problem of determining the (pseudo-hydrostatic) core radius 
(see also below), two other major difficulties arise when analyzing the spectra of 
LBVs. The first one is the determination of r^ff. Onl y if both He I and He II are 
available (as for, e.g., He3-519 and AG Car, ISmit het al 1994) a robust estimate 
may be obtained (although in some cases the He I spectra may be controlled by the 
presence of a binary companion, e.g., in 77 Car). Otherwise, the effective tempera- 
ture has to be obtained from simultaneous fits to lines from species with different 
ionization potentials (e.g. Si II, Mg II, Fe II). This is the case for the tw o LBVs in 
the Quintuplet cluster at the Galactic Center (the Pistol Star and #362, iFiger et al 
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Fig. 16 The possible formation of a "pseudo-photosphere" for a low mass LBV. The difference 
in inner (dashed) and apparent temperature (representative for the size of the computed pseudo- 
photosphere) is plotted against the stellar mass. These computations have been performed for 
a constant luminosity of log L/L.-^. = 5.7. The mass is gradually decreased until the object 
approaches the Eddington limit, and the apparent temperature drops as a result o f low effec- 
tive gr avity and very high mass loss, forming a pseudo-photosphere. Adapted from lSmith et all 
( l2004h . 

I1998L iGeballe et alll2000l lNaiarroll2008h . The second problem arises for low ex- 
citation LBVs and is related to the Helium abundance. Given the low excitation, 
He II recombines in the inner regions and the weak He I lines form very close to 
^(fRoss =2/3). In consequence, a degeneracy of the helium abundance appears 
(as for, e.g., HDE 316285. .Hillier et al 1998 ) as basically for each H/He ratio there 
is a M-value that reproduces the observed H and He I lines equally well. 

Do LBVs form pseudo-photospheres? It is evident that extreme objects such as 
77 Car with present-day mass-loss rates as high as 10^^ M^yr^' form pseudo- 
photospheres, but it is not generally accepted whether the mass-loss rates of the 
more common S Dor variables are large enough to form false photospheres. The 
issue is of particular interest for the following reason: whilst it appears that the 
photospheric temperatures of the objects in Fig. [15] change during HRD transits, 
there is the alternative possibility that the underlying object does not change its 
temperature but changes its mass-loss rate i nstead - resulting in the formation of 
an optically thick wind (e.g. lDavidson|[l987h . As a result of increased mass loss it 
is hypothetic ally possible to form a pseudo-photosphere. Until the late 1980s this 
was the leading idea to explain the color changes of S Dor vari able s. Using more 
advanced NLTE model atmosphere codes. lUeitherer et all d 19891) and lde Koter et all 
(|l996) showed that the color changes consistent with measured LBV mass-loss 
rates are not large enough to make an LBV appear cooler than the temperature 
of its underlying surface. But is the issue of pseudo-photospheres in LBVs com- 
pletely settled? 

In NLTE codes used to investigate pseudo-photosphere formation in LBVs, 
the surface layers are based on the physical relation between L, R, and T, with 
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L = 4nRf^aBT^'^- As the inner radius, Rin, is chosen to lie in the photosphere, the 
input temperature Tj^ does not necessarily equal the output apparent temperature, 
and one therefore defines the "effective" temperature T^-ff at a position where the 
therm alization optical depth equals 1 / a/3 in the optical V band (see de Koter et all 
1 19961 and references therein). For the optically thin winds of OB stars, Tgff equals 
Tin, but for some LBVs that lose mass at rates as high as M ^ lO^'^MQyr^' cases 
might exist where there are significant differences between Ti„ and Teff. If the 
wind is dense enough for the optical photons to emerge from a layer of rapid wind 
acceleration, the LBV can be considered to be forming a pseudo-photosphere. 
This may be favoured by lower LBV masses, as proximity to the Eddington limit 
provides larger mass loss and an increase of the photospheric scale-height. 

Despite the proximity of LBVs to the Eddington limit, current c onsensus is 
that op tically thick winds are generally not present in LBVs. However. lSmith etall 
(120041) studied a special group of "missing" LBVs where the luminosity is rel- 
atively large (with log L/Lq = 5.7) but the stellar masses are possibly much 
lower than those of the classical LBVs (at log L/Lq ^ 6). In other words, the 
lower-luminosity LBVs may well be in closer proximity to their Eddington limit. 
Figure [16] shows the potential formation of an optically thick wind for a rela- 
tively low-mass LBV. The size of the temperature difference (dashed vs. solid) 
is a proxy for the extent of the pseudo-photosphere. The figure demonstrates that 
for masses in the range 15-25 Mq, the winds remain optically thin, but when the 
stellar mass approaches 10 Mq and the star enters the regime near the Edding- 
ton limit, the photospheric scale-height blows up - resulting in the formation of 
a pseudo-photosphere. It should be noted that the reality of this scenario hinges 
critically on the input stellar mass, and LBV masses are poorly constrained from 
observations. 

TTie line-driven winds during S Doradus variations. Whilst most Galactic and 
Magellanic Cloud LBVs have been subject to photometric monitoring, only few 
have been analyzed spectroscopically in sufficient detail to understand the driv- 
ing mechanism of their winds. As reported above, mass-loss rates are of the or- 
der of 10^^ - 10^^ Mgyr^', whilst terminal wind velocities are in the range 
^ 100 — 500kms~'. Of course, these values vary with L and M, but there are 
also indications that the mass loss varies as a function of Tgff when the S Dor vari- 
ables transit the upper HRD on timescales of years. It is this aspect that provides 
us an ideal laboratory for testing the theory of radiation-driven winds. 

The Galacti c LBV AGCar i s one of the more comprehensively monitored 
S Dor variables. IStahl etall (1200 Ih investigated its mass-loss behaviour d uring the 
1990s - modeling the line profile in detail. Ivink and de Koted ([2002) showed 
that the Stahl et al. mass-loss rates rise, drop, and rise, in line with radiation-driven 
wind models in which the variations in mass-loss rate and wind velocity are at- 
tributable to ionization shifts of Fe - the dominant element driving the wind. 

It is relevant to mention here that this variable wind concept (basically wind 
bi-stability) has been put forward as a potential explanation for the circumstel- 
lar density variations inferred from quas i-periodic modul ations in the radio light- 
curves of some transitional supernovae dKotak and Vink| [2006). Additionally, the 
same mechanism is able to account for wind-velocity v ariations in the P C ygni 
absorption line spectra of supernovae such as SN2005gj (iTrundle et alll2008[) . As 
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current stellar models predict massive stars with M > 30 Mq to explode at the end 
of the WR phase, rather than during or after the LBV phase, the ultimate impli- 
cations could be far-reaching - impacting even our most basic understanding of 
massive-star evolution. 

Continuum driven winds during giant eruptions. Whilst, during "quiet" phases, 
LBVs lose mass most likely via ordinary line-driving, they also appear to be 
subject to one or more phases of much stronger mass loss. For instance, the 
giant eru ption of Tj Car with a cumulative loss of ^10 Mp between 1840 and 
1860 (Smiihetal||2003b) which resulted in the Homunculus nebula corresponds 
to M w 0. 1 -0.5 Mpyr" \ which is a factor of 1000 larger than that expected from a 
line-driven wind at that luminosity. Such a strong mass loss has been frequently at- 
tributed to the star approaching or even exceeding the Eddington limiO but line- 
driven mass loss is difficult to invoke, since the formal divergence of M for F ^ 1 
(Eq.O is accompanied by a vanishing terminal speed, Voo °= Vesc °= (1 — F)^-^ — > 0. 
(Remember that the wind-momentum rate is almost i i idependent of F). 

Building upon pioneering work by lShavivl (fl998l . l200lh . lOwocki et"ail (l2004) 
investigated this problem in some detail and developed a simple theory of "porosity- 
moderated" continuum driving in stars that exceed the Eddington limit. At first, 
they showed that the divergence of line-driven mass loss for F w 1 is actually lim- 
ited by the photon tiring associated with the work needed to lift material out of 
the star's gravitational potential. The authors argue that continuum-driven winds 
in super-Eddington stars with a mass loss close to this tiring limit, Mtir = LJ 
{GM^/R^), will result in a stagnating flow at a finite radius, which should lead 
to extensive variability and spatial structure. This should also be true for other 
instabilities (such as a convective one), which are expected to be present in the 
envelopes and atmo spheres of stars close to or beyond the Eddington limit. 

As first noted by lShavivl (Il998h . the porosity of such a structured medium can 
reduce the effective coupling between the matter and radiation by lowering the 
effective opacity in deeper layers (i.e., f^ff < 1 for F > 1), thus enabling a quasi- 
hydrostatic photosphere, but allowing for a transition to a supersonic outflow when 
the structures become optically thin, i.e., f^ff — > F for decreasing optical depth. As 
detailed in Sect. 16.3.21 the effective opacity can be approximated by 



for the optically thick and thin case, respectively, where k" is the mean opacity 
allowing for micro-clumping, Tc is the continuum op tical depth and h = /fi 
the porosity length introduced by lOwocki et all ( |2004|) with / the size and L the 
separation of individual clumps. In other words, the porosity length is the photon's 
mean free path for a medium consisting of optically thick clumps. 

Based on this porosity length and the Ansatz that it should scale with the grav- 
itational scale height H (in analogy to the mixing length theory), the authors de- 
rived scaling laws for a porosity-moderated, continuum-driven mass-loss rate from 
stars with F > 1 . For a super-Eddington model with a single porosity length hKiH 






(47) 



here and in the following, we will consider a generalized Eddington limit, accounting for 
Thomson scattering and bf/ff continuum processes. 
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and r » 1, M w L^/{ac), consistent with the earlier findings bv lShavivl ( l200lh . 
This is much higher than typical for line-driven winds but still only a few percent 
of the tiring limit. Even higher M as implied from the ejecta of 77 Car could be 
obtained by a power-law porosity model (in analogy to the line-strength distribu- 
tion) which should mimic a broad range of clump scales. For power law exponents 
Up < I, M becomes enhanced over the single-scale model by a factor increasing 
with F'/"''~^ approaching the tiring limit under certain circumstances. Together 
with quite fast outflow speeds, Voo w i^(vesc0 and a velocity law corresponding 
to j3 = 1, the derived wind structure can actually explain the observational con- 
straints of giant outbursts in T] Car and other LB Vs. Moreover, the porosity model 
retains the essential scalings with gravity and radiative flux (the von Zeipel the- 
orem, cf. Sect. 12.2.1b that would give a rapidly rotating, gravity-darkened star an 
enhanced polar ma ss loss and flow speed , similar to the bipolar Homunculus neb- 
ula. As stressed by ' van Marie et a il(l2008h. continuum driving does not require the 
presence of metals in the stellar atmosphere. Thus, it is well-suited as a driving 
agent in the winds of low-metallicity and First Stars and may play a crucial role 
in their evolution. 

4.2.3 WR-stars 

The subject of WR stars was recently reviewed bv lCrowttie^ ( l2007t) to which we 
refer for details. Here we briefly discuss some of the recent progress in WR wind 
research. 

WR stars can broadly be divided into nitrogen-rich WN stars and carbon (and 
oxygen) rich WC (and WO) stars. The principal difference between the two sub- 
types is believed to be that the N-enrichment in WN stars is merely a by-product 
of H-burning, whilst the C in WCAVO stars is a sign of the fact that He-burning 
products have reached the stellar surface. As a result, WC stars are thought to be 
more evolved than WN stars. 

We re-iterate that a classification as a WR star reflects purely spectroscopic 
terminology - signaling the presence of strong and broad emission lines in the 
spectra which are produced by strong winds. Such spectra can originate in evolved 
stars that have lost a significant amount of their initial mass, or alternatively from 
an object that has formed with a large initial stellar mass and luminosity. This 
latter group of WR stars may thus include objects still in their core H-burning 
phase of evolution. 

During the last decade, we have witnessed at least two important develop- 
ments in the study of WR stars. One advancement concerned the inclusion of 
line-blanketing in NLTE models which, contrary to the situation for the O stars 
(see Sect. 14.2.11). led to an increase in T^ff and L for WR stars ( ISchmutzlll997l 
iHillier and Milleii[l999l) : The strong dependence of the ionization equilibrium of 
the W R wind on the interaction of the He 11 La line with metal lines close in wave- 
length dSchmut z 1997) a nd the correct trea tment of the He II n^l series near the 
He II Lyman limit ( Hillier and Milledll998h lead to a lower ionization in the wind. 
Thus, a higher temperature is required to reproduce the observed WR lines. 

A second development was the result of the realization that WR winds are 
clumped (e.g. Moffat et al. 1988) leading to a decrease in the empirical WR mass- 



calculated here without the correction for F. 
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Fig. 17 Comparison of mass-loss rates from WR and Galactic OB-supergiants. Solid and dotted 
lines represent mean relations for H-poor WN (solid) and WC stars (dotted) as provided by 
INugis and Lamers (2000), with \Y = 0.9?,. Z = 0.02} and {Y = 0.55, Z = 0.45}, respectively. 
See also Fig. 7 in Crowther (2007). The dashe d line corresp o nds to Galactic O-supergiants, 
and has been taken from the WLR derived by IMokiem et all j2007bl see Fig [13}, assuming 
prototypical values for Voo(2000 kms^') and /f*(15 Rq). 



loss rate by a factor 2 — 4 (e.g.. lHamann and Koesterkell 19981) . Because on the one 
hand mass-loss rates have been down-revised, whilst on the other hand luminosi- 
ties have been revised in an upwards direction, the burden for radiative acceler- 
ation to drive WR mass loss to values exceeding the single-scattering limit has 
become less severe. Nevertheless, compared to 0-stars, WRs display significantly 
larger wind densities and mass-loss rates (factors ^ 10 for the latter, see Fig. 117b. 
In contrast, the terminal velocities of WR stars are comparable to those from O- 
stars, of the order of 1000 km s~' (lat e WNs and WCs ) to 2500 kms^'and higher 
(early WNs and WCs), see Table 2 in lCrowtheJ (|2007|) . 

It is relevant to note that despite the fact that WR temperatures are now cal- 
culated to be larger than a decade ago thanks to the inclusion of metal-line blan- 
keting, the inner core radii (which are a factor of several smaller than the "photo- 
spheric" radii) computed from the most sophisticated NLTE models are still a fac- 
tor of several (^^3) larger than those predicted for the He-main sequence by stellar 
evolution models. In other words, the problem of how WR winds are connected to 
their s tellar core is still lingering, despite the impressive progress recently reported 
(ICrow ther 2007). 

One might now question why there is a difference between the winds from 
WR and 0-stars ;/ also WR winds were accelerated by radiative line-driving, as 
supported by a number of investigations (for references, see Sect. 13.4.3b . This 
problem has been extensively discussed in the literature, but cannot yet be conclu- 
sively answered. However, there are some facts and findings which might serve as 
a guideline. 

(/) Since for a given luminosity the initiated mass-loss rate is considerably 
larger than in 0-stars, whereas the terminal velocities are quite similar, the mo- 
mentum transfer in WR winds must be more efficient. This efficiency is supposed 
to be related to very effic ient multi-line scat tering even ts, possible because of pho- 
ton trapping dLucv and A bbott 1993, Spring mann and Puis 1998): In 0-star winds 
the ionization equilibrium is almost frozen-in, leading to gaps in the frequency 
distribution of the lines and thus to a "premature" photon escape when such gaps 
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are encountered. In WR stars with a high core temperature and a larger wind- 
density, the ionization equilibrium is stratified, and the previous gaps are closed 
due to the emergence of new lines from a new ionization stage. Thus, photons 
become trapped instead of escaping, and can transfer more momentum, (n) More- 
over, again due to the larger wind-densities, the ionization equilibrium in the outer 
WR winds is shifted towards lower stages than present in the corresponding O- 
star winds, and the number of lines is larger, leading to an efficient acceleration 
of the outer wind. (Hi) Following the ionization stratification inwards, there is a 
(hot) iron opacity peak a round 160 kK in the deeper atmospheres of WRstarS 



dNugis and Lamersl[2002h which according to iGrafener and HamannI (l2005h pro- 



vides the required opacity for initiating WR winds (as shown for a WC5-type 
object), (iv) Finally, it should be mentioned th at theoretical model s of WR stars 
display the strongest strange-mode iristabili ties dGautschy and Glatzel 1990; for a 



review, seelOlatzel and Chemigovski 200 ih found so far dGlatzel et aL.1993.). Re- 



cent non-linear simulations ( Wende et al I2OO8I) have shown that Wolf-Rayet mod 



els can reach radial velocity amplitudes which amount to up to 30% of their escape 
velocity. The corresponding mechanical energy transferred by shock waves to the 
stellar atmosphere might easily help or even dominate the ini tiation of mass-loss 
and the acceleration of the lower WR wind. But note also that 'Moffat et al| (|2008|) 
found no coherent oscillations in the hot WC-star WR 111 after photometrically 
monitoring this object with the microsatellite MOST for three weeks. The derived 
upper limit for coherent Fourier components turned out to be almost two orders of 
magnitude below corresponding predictions from strange-mode pulsation simula- 
tions. 

Re turning to the (more conventional) WR wind model by lOrafener and HamannI 
d2005h . let us stress that there are two crucial aspects which had to be adopted in 
order to enable a consistent wind driving: a rather large clumping in the outer 
wind (/ci= 50, but see Sect. 16.3.2b . and the proximity of WRs to the Eddington- 
limit. The latter was suggested by these authors as the most critical parameter for 
initiating a strong mass flux (see Eq. [T4| and Sect. 13.3b . Note that there are no 
present self-consistent models of WR winds with the underlying star being (far) 
away from this limit ! 

We now turn to the empirical dependence of WR winds on initial metallicity. 
Until a few years ago, there were both indications for as well as against the exis- 
tence of a mass-loss-Z dependence for WR winds, however quantitative numbers 
on the exponent m of the mass-loss Z exponent were not provided u ntil theoretical 
results of a M-Z dependence for WR stars became available ( Vin k and de Koteil 
I2OO5I and Sect. 13.4.^ . and empirical exponents have been quoted for both WN (m 
^ 0.8) and WC (m ^ 0.6) stars - similar to those predicted. 

An oft-u sed formula in stell a r evol ution computation for WR winds is the one 
provided by iNugis and LamersI d2000l see also Fig. \Tt} in which the mass-loss 
rate depends on both the He (Y) and the total metallicity (Z). We note that the Z 

in this formula represents the total Z, w hilst, theoretically, Fe is thought to be the 

dominant wind driving ion for WR stars dVink and de KoteJ2005LlGrafener and HamannI 
IIOO8). For this reason we s hould not attach too much physical significance to the 
Z dependencies quoted by Nugis and LamersI ((2000). Nevertheless, the formula 
may well provide a proper representation of empirical mass-loss rates (cf. Fig. 7 

Additionally, there is also a weaker, (cool) opacity peak around 40 kK. 
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in lCrowthe3l2007h . and since a general theoretical description of mass-loss rates 
from WR stars is no yet available, they probably present the most reliable formu- 
lae for use in stellar evolution calculations of He stars. 

Particularly relevant for the issue of LGRBs are the rotational properties of 
their progenitor stars. Owing to the broad emission lines in their spectra, tradi- 
tional V sin i studies are not feasible for WR stars. Luckily, there is a powerful al- 
ternative in the form of linear spectropolarimetry. Apart from its role in constrain- 
ing wind clumping properties (Sect. l63^ , large scale asymmetries r esulting from 
rotatio nall y-flattened w inds may also be probed with the technique. Harries et all 
(Il998h and lVinkI (l2007h studied samples of resp. Galactic and LMC WR stars and 
found equal incidences of polarization "line effects" in WR stars at Z© and '-^0.4 
Z0. As there exists no evidence for significant differences in rotational properties 
of the presumed progenitors of LGRBs, these results suggest that one needs to 
study lower than LMC Z environments to detect a significant population of rapid 
rotators at the end of their lives, if there is a low Z bias for the creation of LGRBs 
via the collapsar model. 



5 The weak wind problem 

First statistically relevant evidence for the presence of 0-stars with theoretically 
unexp ected thin, weak winds was provided already by IChlebowski and Garmanyl 
(Il99lh . From studies of the wind and X-ray properties of 0-stars, they found that 
the mass-loss rates derived for the 0-dwarfs of later spectral type were signifi- 
cantly lower (up to a factor of 10) than expected from the M^U power-law rela- 
tion (Sects. inl and [3]l known for brighter 0-stars. In the same year.l Kudritzki et all 
([l991) analyzed the UV spectrum of the B-type giant j3 CMa (Bill) and con- 
cluded that its theoretically predicted mass-loss rate (3 • lO~^A/0yr~') had to be 
reduced by a factor of 5 to reproduce th eSilv A1397, Cll l A 1 175, and Civ A 1550 
profiles. This result was confirmed by iDrew et all d 19941) from an analysis of the 
ROSAT X-ray emission of /3 CMa, who claimed a similar reduction for the B2II 
star e CMa and presented two alternative explanations for this behaviour. The first 
one invoked the X-ray emission from an "outer zone of hot shocked gas" to pen- 
etrate far enough inwards to increase the ionization around the critical point, thus 
reducing the radiative acceleration and the mass-loss rate. Alternatively, the source 
of wi nd heating could be due to ion-drag effect^ (Springmann and Pauldracl^ 



19921). Usi ng corresponding scaling relations from Gavlev and Owockl (.1994) . 



Drew et al argued that the winds of these two B -giants were not significantly af 



fected, leaving wind-shocked X-rays as the most pl ausible heat sour ce. 

Using Ha as a more reliable mass-loss indicator, IPuls et all ([19961 showed that 
the observed WLR for 0-type dwarfs exhibited a severe curvature towards very 
low wind-momenta for luminosities lower than 10^^ Lq, indicating a break-down 
of the WLR in that domain. It should be stressed, however, that their sample was 
biased towards higher luminosity stars for which H^ could be optimally used. In 
fact, their H^t -diagnostics could provide only upper M-limits for the lower lumi- 

frictional heating accompanying the collisionial momentum transfer between radiatively 
accelerated metals and the bulk wind matter, see page|6| 
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Luminosity [Lj log(L/LJ 

Fig. 18 The weak wind problem. Left: Mass-loss rates for 0-dwarfs in NGC 346 iBouret et aj| 
12003. LMC, filled circles), compai'ed to "normal" winds from Galactic (squares) and other SMC 
stars (triangles) as analyzed by Puis et al ( 1996). Right: Modified wind momenta as a function 
of stellar luminosity for O-stars (adapted from Martins et al 2004). Filled and open symbols 
correspond to Galactic and LMC (small red)/SMC (blue) stars, respectively. Luminosity class 
V, III and I are represented by asterisks, squares and circles. Note the low wind-momenta of the 
SMC objects and the Galactic star 10 Lac (large triangle). 



nosity objects, and subsequently the M-determination in weak-winded stars was 
deferred until unified atmosphere models were able to do so. 

From a theoretical point of vi ew, physical possibilities to explain this be- 
haviour continued to be explored (Pui s et allll998l lOwocki and Pulslll999l iBabell 
ll995Lll99"^lKrticka and KubaL2000) . where thre e effects were considered a s par- 
ticularly promising and have been reviewed by iKudritzki and Pulsl (|2000|) . The 
first one relates to the potential decoupling of metal ions from the bulk plasma, 
which should happen below a critical density when the Coulomb-collisions be- 
come ineffective (Sect. ITTI and references therein). The second one refers to the 
shadowing of the wind-driving lines by the photospheric ones, which consider- 
ably lowers the line forc e, resulting in a reduction of M, especially in the case 
of B-dwarfs (lBabellll996 ). Finally, for low mass-loss rates where the continuum 
is thin throughout the transonic region, curvature terms of the velocity field lead- 
ing to source function gradients may reduce the line-accelerations to values much 
smaller th an in the standard computati ons, resulting again in reduced mass-loss 
rates (Puis et aliri998L lOwocki and Puis 19991 . Unfortunately, due to the lack of 
reliable M estimates it remained open to which extent the above processes could 
account for the weak wind problem. 

First res ults with modern un ified codes using both UV and optical lines were 
obtained bv iHerrero et all (l2002h . For the Galactic 09V star 10 Lac they derived 
an upper limit of M< lO~^M0yr~^ from the optical and lO~"'M0yr~' by 
including additional constraints from the UV, which is m ore than one dex b elow 
the theoretical pred ictions (see also FigfT9]l. Subsequently. iBouret et all(l2003h and 
iMartins et all ( |2004|) obtained similar results for a larger sample of SMC 0-dwarfs 
(see Fig.llSll. The latter authors studied in detail the reliability of their mass-loss 
rate estimates derived from UV spectra and concluded that, at most, the true M 
could be underestimated by a factor of six. With these determined mass-loss rates 
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Fig. 19 The role of X-rays in the O-dwarf weak wind problem. Model fits to C iv 1 1550 from 
the 09V star 10 Lac, using different mass-loss rates. Without X-rays, an M-value as low as 
2.5-10^"' (dashed blue, M inMeyr^') is required, whilst a value of 10^'" (dotted red) produces 
a much too weak profile, and a value of 10^' (dashed-dotted orange) results in a too strong, 
P Cygni type profile. If X-rays are included (Lj/L,~ 10^*-^, dashed-dotted-dotted green curve), 
however, the observations can be reproduced with M = 10^' Moyr^'. Adapted from Najarro, 
Puis & Herrero (in prep, for A&A). 



iMartins et all revisited in detail the physical possibilities discussed above and con- 
cluded that none of them could account for such a strong mismatch between ob- 
served and theoretically expected M. They speculated that although the physical 
mechanism leading to such weak winds remained unknown, the low mass-loss 
rates were intrinsically related to the youth of the stars, possibly testifying the 
phase of wind onset in young 0-stars shortly after their formation. 

iMartins et all (l2005bh extended this study to Galactic objects by performing a 
combined UV and optical analysis of a sample of 12 0-dwarfs. Again, a mismatch 
as high as a factor of 100 was found for later spectral types with log(L/L0) < 
5.2. A major step was the investigation of the effects of X-rays, magnetic fields, 
adv ection and adi abatic cooling in the determination of M from UV spectra. This 
lead iMartins et all to conclude that the inclusion of X-ray emission (likely related 
to magnetic mechanisms, see Sects. IZ2!2l and l6!5. lb in models with low density is 
crucial to derive accurate mass-loss rates from UV lines, while it was found to be 
unimportant for high density winds. This is shown in Fig{T9] where models with 
and without X-rays and different mass-loss rates are compared for C IV A 1550 
from 10 Lac. 

One important conclusion from the above studies is that the discrepancy in the 
derived mass-loss rates seems not to be related to metallicity, as it is present in both 
the Galactic and SMC stars. Interestingly, however, the break-down of the WLR 
starts at earlier spectral types (higher luminosity) for the lower metallicity envi- 
ronment of the SMC ( around 06-6.5 ViBouret et alll2003l [Martins et alll2004l) than 



for the Galaxy (09V, iMartins et alll2005bl) . Thus, although the principal problem 
appears to be unrelated to metallicity, its onset (with respect to limiinosity) could 
be linked to the metal content. 
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Fig. 20 The M-X-ray degeneracy. Model profiles computed for different mass-loss rates and 
scaled X-ray fluxes for strategic UV lines from 10 Lac, with M- values in units of M0yr^'. 
Adapted from Najarro, Puis & Herrero (in prep, for A&A). 



Not only the mass-loss rates, but also the termin al velocities seem to behave 
unexpectedly in weak winds. The values derived bv iMartins et all ( l2005 b) lie far 
below the expected and observed average value of Voo/vesc~ 3 for stars above the 
bi-stability jump, rather close to a ratio of unity. Given the low wind densities, 
only upper limits could be obtained from C IV A 1550, which constitutes the line 
with the highest sensitivity to M and v«, for these late spectral type s. Inverting 
the scaling relations for radiation driven winds (Eq. [13]|, ' Martins etall argued that 
a value of a ^ 0.3 was required to reproduce the observed ratios, again quite 
different from the expected a ^ 0.6 value. Inte restingly, the mechanism related to 
X-ray emission proposed bv lDrew et all d 19941 . see above) was recalled to explain 
such a low value. 

A crucial point regarding the reliability of UV M-determinations for weak- 
winded stars is posed in Fig. |20| Here, we show the diagnostic UV and pro- 
files for a set of models where the mass-loss rate is varied by almost two orders 
of magnitude whilst the X-ray luminosity is increased in parallel to keep the ion- 
ization structure and thus the UV lines at the observed level as far as possible. 
Note that the changes imposed on the wind do not reach the photospheric lev- 
els, which is obvious from the fact that the underlying UV iron forest does not 
change. Apart from the cores of Si IV /1 1397 and O V A 1371, the observations can 
be equally well reproduced with any M combined with an appropriate Lx. Thus, 
no independent UV "measurement" of mass-loss rates may be obtained unless the 
X-ray properties of the star are accurately known. 
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Fig. 21 Left: Synthetic Ha profiles for z Sco, at different mass-loss rates M = 2 • 10 ',2-10 ^ 
and 6- 10^** M,.,yr^' (dotted, bold, and dashed, respectively). Note the insensitivity of Hq at 
lowest mass-loss rates. Right: Synthetic Bra profiles, for the same models. The line is a sensitive 
diagnostics for M even at lowest values. 




Fig. 22 Synthetic Bra profiles for 10 Lac for the model sequence displayed in Figurj20l Note 
the sensitivity of Bra on M even at lowest values. At larger M, the He ll-component becomes 
sensitive on the amount of X-ray emission as well (see text). 



Consequently, a reliable M-diagnostics is required that does not suffer from 
effects such as X-rays, advection, or adiabatic cooling. This is provided by the 
BTg line in the infrared. The relevance of Bra has been pointed out already by 
lAuer and MihalasI (Il969l) . who predicted that even for hydrostatic atmospheres the 
(narrow) Doppler-cores should be in emission, superimposed on rather shallow 
Stark-wings. The origin of this emission is a combination of various effects in 
the upper photosphere and the transonic region that depopulate the lower level of 
Bra, « = 4, stronger than the upper one, « = 5. As detailed by Najarro, Hanson & 
Puis (in prep, for A&A), this is the consequence of a typical nebula-like situation, 
due to the competition between recombinations and downwards transitions, where 
the strong decay channel (4-^3) (stronger than 5^4 itself) is of particular 
importance. 

As shown in Fig. |2l] (right), the emission component of Bra (an d the line 
wings !) react strongly on M, particularly for very weak winds (see also .Najarro et all 
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Il998h . In fact, if we now investigate the same M-Lx combinations as in Fig.l20lus- 
ing Bra, we obtain a really promising result. From Fig.|22]we note that Bra reacts 
to changes in M even at lowest values. Evidently, for such low values (with a deep- 
seated line formation region), the line profile is basically not affected by X-rays. 
With increasing mass-loss rate, however, the line starts to react on changes both in 
M and X-rays, due to an outwards moving line formation zone, reaching the lay- 
ers to which X-rays can penetrate. This is clearly seen in Fig. |22] from the strong 
reaction of the Hell emission component of Br^, which is caused by the large 
sensitivity of the He ll/He III ionization equilibrium to X-rays in the region where 
the emission arises. Our models indicate that under "normal" X-ray luminosities 
(Lx/Lboi < 10~^) the (hydrogen) core of Br^ is not affected though. 

We thus expect that Br^ will become the primary diagnostic tool to mea- 
sure very low mass-loss rates at unprecedented accuracy, thus allowing to clar- 
ify whether they actually lie well below the theoretical predictions and to which 
degree. 



6 Inhomogeneous winds 

In our discussion of diagnostical methods (Sect. 14. lb and recent results (Sect. 14.21) 
we have encountered a multitude of findings which are difficult or even impossi- 
ble to reconcile with the standard model of a smooth and stationary wind. Even 
the inclusion of rotation and magnetic fields does not help to explain additional 
features, at least not at first glance. In the following, we will summarize various 
observational findings (ordered from low to high energies) which point towards a 
structured and variable wind. Further evidence will be outlined in Sect. |6.3| where 
we will concentrate on wind-clumping, maybe the most important aspect with 
respect to deviations from the standard model. 



6.1 Observational background 

Non-thermal radio emission has been discussed on page|40l According to recent 
investigations, the chances are high that this emission might be related to colliding 
winds alone. Insofar, we refrain here from considering non-thermal emission as an 
indication for shocks being present in the wind, in contrast to previous reasoning 
by various authors (including ourselves). 

Electron-scattering wings of recombination lines in Wolf-Rayet stars are weaker 
than predicted by smooth models, indicating the presence of clumping dHilliej 
[1991). Due to a routine inclusion of the latter process into recent spectroscopic 
analyses of WR-stars, perfect fits have been obtained for (almost) all lines, and 
the derived mass-loss rates are a factor of two to three lower than if derived by 
means of unclumped models. Note that this diagnostics cannot be used in OB- 
stars, due to the much lower wind-density and the corresponding lower optical 
depth in electron-scattering. 

Optical line profile v ariability has been detected in the wind lines of OB-stars (Hg: 
Markova et all2005l and references therein; He II 4686: Grady et aJl"983l . lHenrictii 
199lh and BA-supergiants (.Kaufer et al 19961 Typical time-scales range from few 
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hours/days (wind-flow times) to several months. Direct evidence for outward ac- 
celerating wind inhomogeneities (clumps) is implied by the "moving bum ps" su- 
perimposed on the emission lines of Wolf-Rayet stars (e.g.. iRoberd 1 1 994h . whilst 
similar evidence for 0-stars has been reported for only few objects (see Sect. l6.3b . 

Not only wind lines, but also lines which for m to a significant or m ajor part 
in the photosphere (e.g.. He I A5876) are variable. iFullerton et all (1 19961 . and ref- 
erences therein) studied a large sample of 0-stars and found that almost 80 % of 
the investigated stars (all supergiants and a few late-type dwarfs) displayed photo- 
spheric line profile variability, with an amplitude increasing with stellar radius and 
luminosity. These findings were interpreted as being likely related to strange-mode 
oscillations (pagel59l), linking photospheric variability and structure formation in 
the wind. 

Spectropolarimetric variations were detected bv lLupie and Nordsieckl (Il987l) in a 
sample of 10 OB supergiants which pointed to the presence of clumps. Given that 
the position angles of the polarization, and thus the orientations of the electron 
scattering plane, were found to vary randomly, this implied that the polarization 
changes were not attributable to an orderly interaction with a companion, but the 
data provided direct evidence for the presence of an inhomogeneous wind struc- 
ture. 

Black troughs in saturated UV P Cygni profiles strongly point towards the pres- 
ence of multiply non-monotonic flows, and can be simulated in stationary models 
only if allowing for a supersonic velocity dispersion (increasing with distance 
from the star) as the dominating line-broadening agent (page[36]l. 

UV line profile variability is evident in UV P Cygni lines, with large variations 
of the blue edge, contrasted to an almost stationary red emission part (Henrich^ 
ll99lLlPrinialll992l) . implying a significant fraction of small-scale structure. 

Discrete absorption components (DACs) detected in the absorption troughs of un- 
saturated P-Cygni profiles are the most intensively studied indicators of wind 
structure and variability. To date, they are believed to be formed in co-rotating 
interaction regions (CIRs) and will be discussed in Sect. |6.4[ together with the 
so-called "modulation features". 

X-ray emission from hot stars has been measured mostly at soft energies (0.1 to 
^2 keV), either at low resolution in the form of a quasi-continuum or at high 
resolution allowing to see individual lines. The latter will be covered in the last 
part of this section ( 16.5b . 

Already the first X-ray satellite observatory, EINSTEIN, revealed that 0-stars 
are soft X-ray sources. Follow-up investigations, mostly by ROSAT, have con- 
firmed this notion and al so allowed for quantifying the X-ray properties for a mul- 
titude of OB-stars (see .Kudritzki and Puis 2000 and references therein). Roughly, 
the (soft) X-ray luminosity scales with w 10~^Lboi. First there appeared to be 
a rather large scatter indicating the influence of additional parameters, however a 
more recent XMM-Newton study by Sana et al 2006 of the young open cluster 
NGC 6231 showed that the majority of 0-stars display a very limited dispersion 
(of 40 or only 20 per cent depending on whether cool dwarfs are included or not). 
This suggests that the intrinsic X-ray emission of normal 0-stars is rather constant 
for any given object and that the level of X-ray emission is accurately related to 
the basic stellar and wind parameters. 
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The source of the X-ray emission is widely believed to be shocks embedded in 
the stellar wind and related (though not necessarily directly, see Sect. 16.2b to the 
line-driven instabilit y. In term s of a stationary descriptiorPH, a simple model (e.g., 
iHillier et~a l 1993, C assinelli e t al 1994) assumes randomly distributed shocks where 
the hot shocked gas (with temperatures of a few million Kelvin and a volume fill- 
ing factor of the order of 10^^. . .10^^ in the case of denser 0-type winds) is col- 
lisionally ionized/excited and emits X-ray photons due to spontaneous decay, ra- 
diative recombinations and bremsstrahlung. The ambient "cool" stellar wind (with 
a kinetic temperature of the order of T^ff) can re-absorb part of the emission due 
to K- and L-shell processes if the corresponding optical depths are large. This 
simple model (sometimes extended to account for the different post-sh ock cool- 
ing zones of radiative and adiabatic shocks, see Feldmeier et all ll997|) has been 
used to analyze large samples of 0-stars observed with ROSAT TKudritzki et all 
Il996h . and is still used in the context of NLTE diagnostics to account for the 
influence of X-ray/EUV emission on the photo-ionization rates (see below). A 
si mple scaling analysis of the involved X-ray emission and absorption processes 
by lOwocki and CohenI d 19991) showed that the "natural" scaling for optically thin 
(with respect to X-ray absorption) winds is given by °^ (M/voo)^ and for opti- 
cally thick winds by Lx °^ {M/vooY^^ if the volume filling factor follows a radial 
stratification / r\ This would imply, for 0-stars, i w —0.4 to explain Lx °^ ^boi 
by means of the scaling relations for line-driven winds (Eqs. 1131114b . 

Whereas ROSAT observations of OB-star winds showed correlations between 
Lx and Lboh WC stars were not found to emit in the X-ray regime, unless the star 
is part of a binary system and the X-ray emission can be attributed to colliding 
winds. The fact that single WC stars (for WN stars the situation is less clear) are no 
signific ant X-ray emitters doe s not im ply that shocks do not develop in their winds. 
In fact, 'ig nace and Oskinoval ( 11999 ) could explain the lack of X-ray emission in 
single WR winds by realizing that the wind is optically thick to the hot X-ray 
gas. The absence of X-ray detections in single WC stars was reinforced by more 
sensitive XMM-Newton and CHANDRA data, where the lack of X-ray emission 
was attributed to photoelectric absorption by the wind and where the large opacity 
of the WR winds was thought to put the radius of optica l depth unity at hund reds 
or thousands of stellar radii for much of the X-ray band (l Oskinova et alll2003h . 

Since the X-ray and associated EUV luminosity emitted by the shocks is quite 
strong, it can severely affect the degree of ionization of highly ionized species 
such as Cv, Nv and Ovi, by Auger-i onization CMacfarlarie et al|[r993) and even 
more by direct ionization in the EUV (IPauldrach et all 1994 ). In fact, the so-called 



superionization in the winds of O supergiants (exemplarily, the strength of the 
Ovi resonance line in ^ Pup, see Sect. 14. lb can be explained only by means of this 
additional source of photons. A systematic investigation of these effects on the 
complete FUV spectrum, as a function of stellar parameters, mass loss and X-ray 
luminosity has been performed by Garcia (2005). 

Corre lated variability. Finally, let us note that both H„ and U V variability dKaper et all 
Il997i) and and X-ray variability ('B erghoefer et al|[r996l) have been found to be 



correlated, which indicates the propagation of disturbances throughout the entire 
wind. 

Generally, the X-ray lumi nosity of O and early B-stars is constant to within 10 to 20%, see 



uenerally, me A-ray lurm 
iBerehoefer and Schmitlj 1994 
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6.2 Theoretical background: The line-driven instability and time-dependent 
models 



Historically, most of the theoretical effort to explain the above observational fea- 
tures has concentrated on the line-acceleration itself. Later on, however, it turned 
out that at least one of the primary indicators of time-dependence and structure, 
the DACs, can be understood within the concept of CIRs, rather (but not com- 
pletely) independent on the specific physics of the wind acceleration (for details, 
see Sect. 16.4b . 



Linea r stability analysis. Already in their pioneering paper on line-driven winds, 
ILucv an d Solomon (1970) pointed out that the accelerating mechanism of these 
winds should be subject to a strong instability, and subsequent investigations tried 
to quantify its properties, at first based on linear stability analyses. By linearizing 
the time-dependent equations for density and momentum in a frame comoving 
with the mean fluid and using a planar approximation with height co-ordinate z, it 
is easy to show that the conventional ansatz 

P{z,t)\ po(z)+5p(z,0 ] 

v(z,0>= 5v(z,0 Y^Akz-cot) (48) 

5i-ad(z,f) J grad,o(z)+ 5grad(z,0 J 

for the density, velocity and line force (quantities with subscript '0' refer to the 
mean flow, vq = 0, and 5p, 5v, 5grad are the corresponding amplitudes of a sinu- 
soidal disturbance with wavenumber k and frequency ft)) results in 



4ft,5p + poifc5v = ^ Tp=V,j\ S^^^Sg^ 

-ico5v = ^ ft) = i%i I 

ov 



(49) 



In the la tter equation, we have negle cted the pressure terms (for a complete anal- 
ysis, see lOwocki and Rybicki|[l984h . Thus, the phase angle between density and 
velocity disturbance, cos(p, the growth-rate of the instability, Q, and its phase 
speed with respect to the mean flow, V(p, 

cos(p= ~^^^^7/y , X2 = 3(ft)), v^ = ^{(o)/k (50) 

are completely determined by the response of the line-acceleration on the veloc- 
ity disturbance, at least in the supersonic regime where the neglect of the pres- 
sure force is justified (9t(jc) and 3 [x] are the real and the imaginary part of the x 
quantity, respectively). Thus, a correct description of the line force is of highest 
importance for time-dependent simulations (see below). 

Within the Sobolev approximatior0 (page|7]l it turns out that 5g^°^ / 5v ik, 
i.e., it is purely imaginary. Consequently cos 9 = — 1 and density and velocity 
are completely anti-correlated. There is no growth of the instability, = 0, 

but only an oscillatory behaviour, and the phase speed (independent of k, i.e., the 

Remember that most stationary wind models are based on this approach. 
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medium is non-dispersive in this approximation) is inwards directed. The inclu- 
sion of pressure terms leads to slight modifications: the inwards directed phase- 
speed becomes marginally larger (roughly by /V(p), and a slow, outwards prop- 
agating mode becomes possible, at a speed w /V(p < a ( Abbott 1980). 

The (absolute) value of V(p (inwards mode) at the critical point (page|9ll of sta- 
tionary models based on a Sobolev line-force turns out to be identical with the 
corresponding flow speed, much larger than the speed of sound. Ins ofar, this crit i- 
cal point may be considered as a critical point in the "usual" sense (iHolzeJ 19771) : 
It is the outermost point in the flow where a communication with lower layers 
is possible, in this case via the above radiative-acoustic waves which have been 
named as Abbott-waves according to their "inventor". Following Holzer's "wind- 
laws", M then depends solely on processes in the sub-critical region, whilst v«, is 
controlled by the acceleration in the super-critical part. For a detailed discussion 
of Abbott-waves and the corresponding modifications of the above wind-laws for 
line-driven wind we refer the reader to Feldmeier and Shlosman (2002). In any 
case, however, remember that the use of the Sobolev approximation is justified 
only for wavelengths which are large compared with the Sobolev length (page 
[T]), L^°^ = Vth/(dv/dr), i.e., for small k. A discussion of critical points in models 
based on non-Sobolev line-forces is beyond the scope of this r eview. 

A somewhat opposite approach dMacGregor et al 1979, , .Carlbergll 19801) results 
in a purely real response, 5grad = A5v, with 3 (A) = and A > 0, which is valid as 
long as the disturbances remain optically thin. In this case, the finite sound-speed 
is of importance and disturbances can propagate only for large k, in both directions 
and becoming strongly amplified, with a phase-speed of the order of a. 

A combination and generalization of these two alternative approaches has been 
obtained by Owocki and Rybicki ( 1984) in the form of a "bridging-law", which re- 
sults from a perturbation analysis of the line-force. In this case, 5grad/ has both 
a real and an imaginary component, and includes the two limiting cases for small 
and large k from aboveH The resulting phase relation between density and veloc- 
ity is predominantly anti-correlated for all k, and there is a fast, exponential growth 
of the inwards propagating modes, within a dispersive med ium. Important addi- 
tions t o this analysis have been provided by lLucyl (Il984) and lOwocki and Rybickil 
(Il985h . by considering the effects of the disturbed diffuse radiation field and as- 
sociated line-force (the so-called "line drag'EB) which had been neglected in the 
previous approaches. 

Time -dependent models. After these initial studies, effort has concentrated on 
a direct, hydrodynamic al modeling of t he ti me-dependent str ucture of line-driven 
winds (for reviews, seel Owockilll994cl and iFeldmeiedl 19991) . As already pointed 
out, the key aspect concerns a consistent description of the line force which fol- 
lows the local and non-local conditions (due to the coupling with the radiation 
field) at all scales of structure formation. In order to keep the problem com- 
putationally feasible, an "exact" calculation of the line-force is prohibitive, and 

The critical length scale is here the "bridging length" ' , with Xb the line opacity required 
to reach unit optical depth. 

" The disturbed diffuse radiation field in resonance lines leads to a damping effect which, 
at the foot-point of the wind, exactly cancels and otherwise diminishes the instability resulting 
from pure absorption lines. 
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complex integral forms of vari ous degrees of approximation have been developed 
dOwocki and Pulsll 1 996L 1 1 999t) to allow for this objective. 

After incorporating the stabilizing effect of the line drag by the diffuse, scat- 
tered radiation field (see above), hydrodynamical models have proven to be signif- 
icantly more stable in the lower wind than the first ones based on pure absorption 
dOwocki et aJl988h . The outer wind (r ^ 1.3/?*), however, still develops extensive 
structure that consists of strong revers^^ shocks separating slower, dense mate- 
rial from high-speed rarefied regions in between. This is a direct consequence 
of the strong amplification of the inwards propagating mode with a predominant 
anti-correlation of density and velocity, as derived from the linear analysis. Such 
structure is the most prominent and robust result from time-dependent modeling, 
and the basis for our interpretation and description of wind-clumping. Moreover, 
the multiply non-monotonic nature of the velocity field explains the b lack absorp- 
tion troug hs observed i n saturated UV P Cygni lines in a natural way (lLucylll982L 
[^83, Pu is et allll993aLlOwocki|[l9943) . 

iFeldmeied (11995 ^ extended these models by accounting for the energy trans- 
port including radiative cooling. By investigating the impact of various photo- 
spheric disturbances which c an affect the on-set of structure formation (via excit- 
ing the line-force instabilitv). lFeldmeier et all (Il997l) concluded that clump-clump 
collisions (and not the cooling of the reverse shocks themselves, as previously 
suspected, since this results in much too low filling factors) are the origin of the 
observed X-ray emission. 

Recent effort has focused on more-dimensional simulations and on the study 
of the outermost wind. By means of an efficient 1-D "pseudo-periodic box for- 
malism" to investigate the e volution a t distances very far from the central star, 
iRunacres and Owockil (l2005l . see also iRunacres and Owockill2002h showed that 
structure can survive out to distances of more than 1000 stellar radii, due to super- 
sonic clump-clump collisions which counteract the pressure expansion. 

Initial results of 2-D simulations have been obtained using different approxi- 
mations for the calculation of the line-force. In the simplest case (radiatively iso- 
lated azimuthal zones and neglect of lateral line-dra g^ , the shells arising in 1-D 
models are broken up by Rayleigh-Tayl or or thin-shell instabilities , resulting in a 
completely in coherent lateral structure (Dessart and Owockil l2003h . A follow-up 
investigation (jPessart and Owockill2005[) showed that by accounting for the lat- 
eral line-drag and the lateral mixing in a more self-consistent manner (three rays 
at different impact parameters, as already used by Owocki. 1 999) . such models - at 
least at highest resolution - show a much larger lateral coherence than correspond- 
ing one-ray models. Quantitative results are still lacking though. Further effort is 
needed to answer the important question about the lateral and radial length scale 
of clumps, required to check and improve present assumptions when modeling the 
spectra from the IR to the X-ray regime. 



Time-dependent versus stationary approach. Obviously, time-dependent mod- 
els based on a non-local, unstable line-force are able to explain a large number 



Reverse shocks travel backwards in the comoving frame (though outwards in the stellar 
frame). In the pre-shock region (starwards from the shock), the velocity is high and the density 
low, and vice versa in the post-shock region. 

see lRvbicki etail ( 119901) . 
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of observational findings, at least qualitatively. On the other hand, they appear to 
be in stark contrast with our assumptions within the standard model (Sect. 12.1b . 
both with respect to modeling and diagnostics. However, when viewed with re- 
spect to the mass distribution of density and velocity, the stationary and time- 
dependent models are quite similar (e.g., Fig.[32l left panel). Given the intrinsic 
mass-weighting of spectral formation, this explains why detailed line-synthesis 
calculations assuming line-opacities proportional to the local density a re able to 
genera te (time-dependent) line profiles quite similar to observations ( jPuls et all 
Il993ah . 

Let us also emphasize that the gross wind properties derived from time-depen- 
dent models like the terminal flow speed and time-averaged mass loss rate agree 
well with those following from a stationary approach. D eviations from this rule 
due to source function gradients in the transonic region (iOwockiandPulslll999L 
see also Sect. |5]l are not a consequence of the line-driven instability. Instead, they 
originate from an incorrect use of the local, Sobolev approach in the stationary 
models in a region where this approach is no longer justified, mainly because of 
the strong curvature of the velocity fieldF^ 



6.3 Wind clumping 

In our review of "observed" wind parameters (Sect. |4ll, we have frequently re- 
ferred to the process of wind-clumping, and particularly to its impact on the de- 
rived mass-loss rates, without providing further details. Indeed, the investigation 
of this process and its consequences constitutes a major fraction of studies per- 
formed in the last decade, and various contributions regarding the present status- 
quo can be found in t he proceedings of a recent workshop dedicated to this topic 
(iHamann et alll2008bh . 

In the previous two sections, we have summarized important observational and 
theoretical findings which unambiguously indicate that radiation driven winds are 
structured and time-variable, and we have outlined that the line-driven instability 
is a primary candidate for explaining at least part of these findings. Thus, we 
are now in a position to consider wind-clumping in more detail, starting with a 
discussion of more specific evidence. 

6.3.1 Observational evidence 

Half a century of clumping. The presence of cl umping as "extreme density fl uc- 
tuations" around hot stars was first invoked by lOsterbrock and FlatheJ ( Il959l) to 
explain the discrepancy between observ ed and p redicted fluxes in the H II region 
in the Orion nebula. Two decades later. iLucy an d White ( 1980) proposed a phe- 
nomenological model to explain the observed X-ray luminosity of Pup, in which 
the stellar wind would break up into a population of radiatively driven blobs and 
where the radiation should originate from the bow shocks preceding the blobs (se e 
also Sect. 16. 5. Zt . Applying the same formalism as lOsterbrock and Flatheil(ll959t) . 



Remember that standard Sobolev theory requires a vanishing curvature over the coiTespond- 
ing Sobolev length, which is small in the transonic region. 
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lAbbott etail (Il981bh investigated the effects of dense clumps on the thermal free- 
free radio emission of stellar winds as a function of both the volume filling factor 
(/voi<l) and the density ratio between clumped and inter-clump matter. For the 
standard assumption of vanishing inter-clump density (see below), they showed 
that the radio flux is a factor of fm\~^^^ larger than that from a homogeneous wind 
with the same average mass-loss rate. Thus, using Eq.|45] radio mass-loss rates de- 
rived from clumped winds are lower by a factor of fyoi^^ compared to unclumped 
winds. 

Wind-clumping in dense winds. In the mid-80s, two more observational diag- 
nostics suggested that the standard, and convenient, assumption of homogeneous 
stellar winds might fail due to clumping. The first line of evidence came from 
studies of the lig ht-curve eclipses as a funct ion of wavelength in the WR + O bi- 
nary V444 Cyg (ICherepashchuk et allll984l) . The differences in shape and depth 
of eclipses from the UV to the IR could not be simultaneously interpreted using 
a smooth electron density distribution in the wind. The authors required the pres- 
ence of individual density condensations becoming optically thick to reconcile the 
electron scattering (°= p) dominated UV and optical fluxes with the free-free and 
bound-free (oc p^) infrared values. 

The second line of evidence relied on studies of continuum and line variabil- 
ity. For continuum investigations, linear polarization variability is particularly ad- 
vantageous. This process, being directly related to the electron density through 
electron scattering processes, provides information on the presence of localized 
inhomogeneities and is more reliable than pure continuum photometric variabil- 
ity studies, where other processes such as rotational mod ulation or pulsations 
may come into play. In addition to the OB-star work by iLupie and NordsiecS 
dips?), mentioned already in Sect.l6.1l extensive an d systematic monitoring of lin- 
ear po larization in WR-stars dSt.-Louis etalll 19871) and LBVs (Taylor et al"l 99lL 
iDavies et al,2005) revealed stochastic variability on short time scales, indicating 
that the material had to be confined in clumps. More recent results helping to 

quantify the clumping properties are discussed below. 

In the case of lines, time -resolved spectroscopic monitoring of WR-stars dMoffat et all 

I1988L Robert 1994) revealed the presence of narrow emission features on top of 
the broad emission lines that appeared to propagate from the line centers to the 
line wings on timescales of the wind flow-time. These moving sub-peaks were 
interpreted as evidence for blob ejections into the wind and are now considered to 
be one of the most clear-cut e vidences for clumping in WR winds. 

iLepine and MoffatI d 19991) monitored the line-profile variations in the Hell 
emission lines of WN stars and in C III emission lines of WC stars, and studied 
these variations by means of a detailed wavelet analysis. As well, these authors 
performed simulations of line-profile variability (Ipv) assuming a wind consisting 
of a large number 10'*) of discrete wind-emission elements (DWEEs). Obser- 
vations were best reproduced by assuming a large ratio of the radial to the lateral 
velocity dispersion for the DWEEs, consistent with the nature of the line-driven 
instability and the corresponding lateral line-drag (pagelTOll. Physical models of 
such Ipv interpreted in terms of clumps generated by the line-driven instability 
have been provided by Dessart and O wockil (l200 2a.b). where the latter study ex- 
tends the iLepine and Moffat (.1999) wavelet formalism to modeling Ipv from hy- 
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drodynamic simulations. A key result is that one can fit most of the characteristics 
if one assumes a lateral coherence scale of ^3 degrees. 



Wind-clumping in OB-star winds. For 0-stars, observational Ipv-studies are 
more difficult due to the ir roughly factor of te n thinner winds. First investiga- 
tions were carried out by lEversberg et all d 19981) on the He ll 4686 line in Pup, 
which revealed outward-moving inhomogeneities in its wind that started near the 
stellar surface and became invisible at ^ 2/?* . From the wind acceleration, the 
authors derived a value of j3 = 1, consistent with results from independent diag- 
nostics assuming unclumped winds. Thus, it seems that the clumping properties of 
O- and WR stars are s i milar and this "universality" has been recently confirmed 
by iLepine and Moffaj (120081) who extended the study to five hot stars in differ- 
ent evolutionary stages. They show that the moving sub-peaks identified on top 
of the broad emission lines in WR-stars appear also in the Of stars Pup and 
HD 93 129 A, in the more evolved hydrogen-rich and luminous Of-like WN stars 
HD93131 and HD 93162, and in the more mass-depleted WC star in Vel. The 
authors conclude that stochastic wind clumping is a universal phenomenon in the 
ra diation-driven, hot win ds from all massive stars. 

iMarkova et all (l2005h presented a study on the Ipv of Ha for a sample of 15 
0-type supergiants, ranging from 04 to 09.7. All stars showed evidence of signif- 
icant Ipv, mostly dominated by processes in the wind, where the variations were 
found to occur from the wind base on out to ^ 0.3 Voo. From profile simulations, 
they concluded that the mass-loss rate is only marginally variable, resulting in 
insignificant variations of the wind momenta and hence of the WLR. 

Of course, best direct evidence for the presence and behaviour of clumps in 
stellar winds would be reached using imaging devices with enough spatial res- 
olution to "see" the individual clumps close to the star and track them as they 
expand with the wind flow. By such imaging, also ge ometric cons traints on, for 
instance, the sizes of the clumps might be obtained (e.g. lChesneau et aL200Q') . Un- 
fortunately, micro-arcsec resolution would be required to carry out such studies, 
even for the nearby 0-stars. Present observing facilities allow to see such clumps 
only at much larger scales, such as the outer win d regions where the thermal radio 
emission originates from d Williams et all [19971) . However, although these obser- 
vations indicate that clumps seem to survive over the entire stellar wind (see also 
iRunacres and O wocki 2005), no information is provided on their formation region 
nor how they evolve throughout the wind. It should be noted that even after the 
possible commissioning of 100-m class telescopes, one would still not be able to 
spatially resolve the innermost wind regions where the clumps originate, and it is 
thus pivotal that alternative methods to constrain clumping properties are further 
developed, such as linear polarimetry (see below). 

At least for 0-stars, however, we presently have to rely mostly on indirect 
methods to address the behaviour of clumping throughout the stellar wind. Such 
indirec t evidence for the presence of clumping was found in parallel by iMarkova et all 
(l2004l) and iRepolust et al (l2004l) . By means of a comparison between the ob- 
served and predicted WLR for large samples of Galactic 0-stars, these authors 
showed that if the stellar sample was plotted as a function of H^ profile type (stars 
with Ha in emission and those with Ha in absorption, partly filled in by wind 
emission), two different WLRs appeared, inconsistent with the basic theoretical 
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prediction that the wind-momentum rates of 0-stars should depend on luminosity 
alone (Sect.|2]and Fig.[5]l. It was suggeste43 that this difference is related to wind- 
clumping (being "visible" only in denser winds, i.e., in those with in emission) 
and shown that a unification of both "observed" WLRs with the theoretical one is 
possible when the average clumping factor, /ci(=/voi ^ see below) of the clumped 
matter in stars with in emission is of the order of /ci ~ 5, implying a net reduc- 
tion of M by a factor of 2.3. Note that these values differ significantly from those 
derived by means of UV stud ies with net M reductions up to factors of 10 (e.g., 
iBouret et alll2003l . lFullerton et al 2006 . see below). 

Therefore, at this stage, several fundamental questions need to be addressed: 

- Which are the upper and lower limits for the clumping factor obtained from 
observationally driven studies and how reliable are they? 

- What is the stratification of the clumping factor throughout the wind, and what 
is the shape of the clumps? 

- How realistic is the treatment of clumping assumed so far and how well can 
modern, sophisticated atmospheric codes account for it? 

- How much mass-loss reduction can evolutionary models of massive stars cope 
with to still be consistent with independent constraints such as the observed 
numbers of WR-stars? 

Before we can consider these points in detail, we have to describe the basic con- 
cept underlying most of the following analyses. 

"Micro-clumping" . This concept bases on the hypothesis that the wind con- 
sists of small-scale density inhomogeneities, which redistribute matter into over- 
dense clumps and underdense inter-clump matter. Motivated by the principal re- 
sults from hydrodynamical simulations including the line-driven instability, the 
inter-clump matter is assumed to be void0 In this case then, the average density 
can be described by 

(p)=/volP + , (p')=/vol(p + )' (51) 

where p+ is the density inside the overdense clumps, and ( p^) is the mean of th e 
squared density. Thus, the clumping factor as introduced bv lOwocki et all ( ri988l) . 

/ci = (p')/(p>' ^ U=Ur' and p+=fci{p), (52) 

corresponds to the inverse of the volume filling factor and measures the overden- 
sity of the clumps. Since we assume void inter-clump material, matter is present 
only inside the clumps, with density p+, and the corresponding opacity is given 
by K = Kc{fc\{p)), where C denotes the process evaluated inside the clump. Op- 
tical depths have to be calculated via T = / Kc{fc\{p))fvo\'ir, with a reduced path 
length (/voidr) to correct for the volume where a clump is "seen". Note that this 

following a previous idea by Puis etaH J20031) . 

Alternativ es from t his assumption are rare, and not used in standard diagnostical tools until 
recently when IZsargo e t al (2008a) re-investigated the "superionization" of O vi (page|67) in 
^ Pup. Their simulations show that clumped wind models that assume a void interclump medium 
cannot reproduce the observed O vi profiles, though enough O vi can be produced if the voids 
are filled by a low-density gas. 
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expression is valid only if all clumps are optically thin. Extensions for optically 
thick clumps are given in Sect. 16.3.31 Combining these results, optical depths for 
a clumped medium can be expressed by a mean opacity, K", 

T= /)fdr with )f=Kc(/ci(p))/voi = ^K<:(/ci(p))- (53) 

■I Jcl 

Thus, for processes that are linear in density, the mean opacity of a clumped 
medium is the same as in the homogeneous case, whereas for processes scal- 
ing with the square of density, mean opacities (and emissivities) are effectively 
enhanced by a factor /d. Moreover, radiative processes described by the micro- 
clumping approach do not depend on the size, the geometry or the distribution 
of the clumps, but on the radial stratification of the clumping factor alone. The 
enhanced opacity for dependent processes has the consequence that mass-loss 
rates derived by such diagnostics are lower by a factor of vTci than corresponding 
mass-loss rates assuming no clumping. Consequently, the optical depth invariant, 
Q (see Eq.l39b, for such processes might be extended according to 

as long as there is no strong radial variation of /d or /voi. 



6.3.2 Quantification of wind-clumping 

Let us now review how observations and models have evolved during the last 
decade and summarize the most important results concerning the "chase" of the 
clumping factor. 

One of the first ( and stil l most important) attempts to quantify wind-clumping 
was carried out by iHillierl (Il99lh who analyzed the corresponding effects on 
the p^-dependent emission components and the p-dependent electron scattering 
wings in WR emission lines. Noting that standard, homogeneous m odels severely 
overestimated the electron scattering wings of the observed profiles. [Hillieil d 1 99 ih 
showed how both components could be simultaneously reproduced when clump- 
ing is accounted for and the mass-loss rate is reduced by a factor of ^ 2. 

For the method to work, the electron scattering wings have to be strong and 
clearly defined, i.e., a good continuum rectification is required. The method also 
assumes that the strength of the p^-dependent emission can be described by micro- 
clumping alone. Thus, th is metho d may also be applied to "well behaved" LBVs 
like P Cygni, for which iNajarrol (l200lh found /voi « 0.5 by analyzing all the 
Balmer lines. 

Results using the new generation of unified model atmospheres accounting 
for blanketing and clumping started to appear with the ne w millennium (se e Sect. 
14.2.1b . From quantitative IR studies in the Galactic Center. lFiger et all ( 1200 2') found 
that /voi-values around 0.1 (M-reductions of ^^3 with respect to homogeneous 
models) were required to best reproduce the spectra of WNL and Of stars in the 
Arches Cluster. Anal yses carried out for O supergiants in the Magellanic Clouds 
dCrowther et alll200"2l using UV and optical data) revealed that /voi= 0.1 could re- 
produce the observed PV/l/11118- 1 128 doublet wi thout requiring a reduced phos- 
phorus abundance as suggested by IPauldrach et all (fl994 .200 1.) . Further analyses 
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Fig. 23 Effects of clumping and a correspondingly reduced mass-loss rate in the UV O v A 1371 
and NlV/11718 lines of HD 96715 (04((f))). The derived /voi=0.02 value implies a factor of 7 
reduction in M. Dotted: Profiles from an un clumped model wh ich reproduces other, dumping- 
insensitive lines equally well. Adapted from lBouret et all ( 120051) . 



in the SMC found likewise /voi values around 0. 1 dHillier et al 2003h and lowe r 
0.1 - 0.01 in four O stars of the SMC cluster NGC346. ' Bouret et~aill2003h . 
A similar s tudy was car r ied ou t using UV (FUSE -i- lUE) spectra of two Galactic 
04 stars bv lBouret et all (l2005h . In this case they found that values as low as /voi= 
0.02 and 0.04 were required to satisfactorily reproduce not only the P V A A 1118- 
1128 and the O V A1371 lines, but also NlV A1718, as shown in Fig.[23] Though 
indicating an enormous degree of clumping (overdensities of factors between 25 
to 50!), these numbers are consistent with the alternative explanati on in terms of a 
strong oxygen depletion (factor 50) in early O stars, suggested by IPauldrach et all 
(1200 ih as a result of analyses using homogeneous models. 

Interestingly, the lower /yoi -values found for the SMC objects seem to cor- 
respond to the Galactic ones down-scaled by the abundance ratio. Major uncer- 
tainties in these UV-based methods are due to the role of X-rays and the validity 
of micro-clumping (see below) assumed for the models. Other effects such as 
changes in the ionization structure driven by clumping are included, however. 



The P V problem. Although the problem wi th the O V A 1371 li ne had b een exten- 
sively debated in the literature (Ha ser et aJl 998. Paul drach et a l 2001, B ouret et all 
l2003h . it was the advent of the FUSE mission providing the cleaner P V AA1118- 
1128 mass-loss diagnostics that gave rise to the so called P V problem (see lFullerton et ai 
|2008| for a comprehensive review of this issue). Due to the very low cosmic abun- 
dance of phosphorus, the P V XX 1 1 18-1 128 doublet basically never becomes sat- 
urated, not even when P+^ is the dominant ion. This enables to obtain a direct es- 
timate of the product M{q), where (q) is a spatial average of the ion fraction of the 
species (cf. page[36l). Unfortunately, estimates of (q) for any given resonance line 
is problematic unless detailed NUTE modeling is performed, and even then, prob- 
lems such as X-rays or subtle blanketing/blocking effects may lead to incorrect 
predictions. On the other hand, an empirical determination of ionization fractions 
is not feasible, since access to resonance lines from consecutive ionization stages 
are not available in the far-UV and UV regions. In the case of P V, some insight 
is gained from the available FUSE atlases of Galactic and Magellanic OB spectra 
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Fig. 24 roSE spectra of HD 14947 (05If^) and HD 203064 (0 7 5 III) for the P v 1 118- 
1128 doublet, together with the best tits from the simulations by IFullerton et all i2006l) with 
a low M (solid) and the M derived from H^ using homogeneous models (dashed). To unify 
both analyses and assuming (g)(P+'*) = 1, values as low as /vi » 0.03 (HD 14947) and 
/vol ~ 0.003 (HD 203064) are required, resulting in a down-scaling of the "unclumped" mass- 
loss rates by factors of 6 and 18, respectively. 



(IPellerin et all200llWalbom et all2002h . These show that the doublet is very weak 
for the earliest 0-subtypes, exhibits a stable maximum over a wide range of spec- 
tral types, and then starts to disappear rather quickly for stars between 09.5 and 
BO. Such morphology seems to secure a maximum of the ionization fraction 
for mid 0-type stars, a behaviour confirmed by model atmosphere calculations. 
Therefore, for those 0-stars within a safe (q) w 1 region, the P V AA 1 1 18-1 128 
lin e should provide a rob ust and independent estimate of M (pure p -diagnostic). 

IFullerton et all ( l2006l) carefully selected a sample of 40 Galactic 0-stars, for 
which also reliable M-estimates (Ha, radio) were available, and carried out 
a detailed comparison between both methods. Surprisingly (see Figs. I24land 1251) 
they found a severe discrepancy by a factor of ^10 in M for mid-range 0-stars 
and by substantially larger factors for earlier and later spectral types. A more re- 
strictive consideration, taking into account only those stars for which (^)(P"'"^) is 
expected to be 1, i.e., the "safer" mid-range O supergiants, resulted in a median 
discrepancy M(p^)/(M(P V)(^)) = 20, which would imply 400 if interpreted 
in terms of micro-clumping, i.e., extreme clumping, much larger than derived by 
other methods. This discrepancy is the P V problem and is nicely illustrated in 
Fig.|25] At the end of this section we present a possible solution of this problem. 

Polarimetric studies. Soon after the detection of moving subpeaks on top of the 
broad emission lines, other indications for clumping of WR winds were found. Be- 
sides the diagnostics provided by the electron scattering wings, a straightforward 
comparison of mass-loss rates based on p^-diagnostics with rates determined us- 
ing o ther methods, such as those from polarization variations in WR-hO binariefl 
(e.g. ISt.-Louis et suggested that the mass-loss rates determined from 

emission-lines are overestimated by about a factor three (lRobert|[l994l) . which 



The basic principle is that the light from the O-star acts as a probe of the WR-wind via elec- 
tron scattering. O-star photons scatter off the abundant free electrons in the WR wind, which 
polarizes the light received by the observer, where the level of polarization depends on the 
orbital-phase. 
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Fig. 25 Comparison of M (der ived from p^-d i agnos tics) with M{q)(P+*) for a large sample of 
Galactic O-stars. Adapted from lFuUerton et al ( I2006t) . see text. 



agre es with the factors of 2-4 derived from NLTE analyses using clumped models 
(e.g. lHamann and Koester ke 1998). 

Extending the above technique, ISt.-Louis and Moffa3 (l2008h studied the po- 
larimetric variability of O-stars winds in six massive O-l-O binary systems. They 
reported clumping independent mass-loss rate estimates that suggest only small 
clumping corrections. However, the O-stars analyzed by this method had relatively 
weak winds and, as a consequence, the electron scattering optical depth of the 
wind is low. Thus, a significant fraction of the polarization signal originates from 
scattering off the photosphere, implying that the derived mass-loss rates should be 
lowered. 

In their investigation of the polarization variability in WR- winds, [Robert et"ail 

(119891) found an anti-correlation between the terminal wind velocity and the de- 
gree of random, intrinsic scatter in polarization, and interpreted this finding as 
due to the p resence of blobs that survive or grow more effectively in slow winds. 
iDavies etall ( i2005[) found this trend to continue into the regime of the LBVs, 
which have much lower v„ than WR-winds and are an ideal testbed for constrain- 
ing clumping properties, due to their long wind-flow times. 

iDavies et~ail (l2005h performed a spectropolarimetric survey of LBVs in the 
Galaxy and the Magellanic Clouds and found more than 50% of them to be in- 
trinsically polarized. As the polarization angle was found to vary irregularly with 
time, the polarization line effects were attributed to wind clumping. Monte Carlo 
models for scattering off wind clumps were developed bv iCode and Whitneyl(ll995h . 
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Fig. 26 Time-averaged polarization over a broad range of ejection rates per wind tlow-time, 
At ^ ~ 20 , tiie optical depth per clump exceeds unity and the overall polarization falls off (see 
iDavies et"ail !2007 for details). The observed level of polarization for P Cygni is marked by the 
dash-dotted line. Note that there are two ejection-rate regimes where the observed polarization 
level can be achieved. 



iRodrigues and Magalhaesl(l2000l) . lHarriesl(l2000h . whilst analy tical models to pro- 
duce the variability of the linear polarization were presented by Davies et all (l2007h . 

An example of an analytical model that predicts the time-averaged polariza- 
tion for the LBV P Cygni is presented in Fig. |26] The clump ejection rate per wind 
flow-time o/K is defined as ^ = Ntf[ = NR-i,/voo, where the clump ejection rate, 
N, is related to the mass-loss rate via M = NNeHeinn, with A^,, the number of elec- 
trons in each clump and /ip the mean mass per electron. There are two ejection- 
rate regimes where the observed polarization level can be achieved. One regime 
is where the ejection rate is low and a few optically thick clumps are expelled, 
and the other regime is that of optically thin clumps where the number of clumps 
is very large. These two models can be distinguished via time resolved polarime- 
try. Given the rel atively short timescale of the observed polarization variability, 
iDavies et allboOTi favor the latter scenario, suggesting that LBV winds consist of 
thousands of optically thin clumps in close proximity to the photosphere. 



Stratification of the clumping factor. Most of the studies considered so far as- 
sume that clumping starts just above the sonic point and remains constant through- 
out the wind. We now address whether this approximation is feasible or needs to 
be relaxed by requiring stratified clumping. Obviously, this can be achieved by 
identifying reliable clumping diagnostics at different distances from the star. First 
indications that /yoi (or /d) vary as a function of radial distance from the star were 
obtained by Nugis et al (1998) by means of IR to mm and radio continuum stud- 
ies of WRs. To explain a spectral index from mm to cm wavelengths of the order 
of ^0.76 ( s teeper than a value of 0.6 expected for smooth winds, see Eq. |46] |. 
iNugis et"ai (Il998h proposed a stratification of the clumping factor such that /d is 
unity at the wind base, reaches a maximum at 5-10 7?*, and then recovers unity in 
the outer wind parts. With respect to mass-loss rates derived from homogeneous 
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Fig. 27 Radial str atification of the clump i ng fac tor /d for an O-supergiant wind. Solid: Theoret- 
ical predictions bv lRunacres and Owockil ( [2002 h from hydrodynamical models, with self-excited 
line-driven instability. The formation regions of different diagnostics are indicated. Dotted: Av- 
erage clumping factors, as derived for the wind from ^ Pup, assuming an unclumped outer wind 
jPuls et al2006D . Dashed: as dotted, but assuming an outer wind with clumping as theoretically 
predicted. The corresponding mass-loss rate is a factor of (Z™')"'' smaller than in the dotted 
case. See text. 



models. iNugis et ai ( Il998h concluded that these were overestimated by ^^0.2 dex 
for WN stars and by -0.6 dex for WC star0 

A similar parameterization of the clumping factor was used by iFiger et all 
(j2002) in their analysis of WNL and Olf + stars in the Arches Cluster. Only a van- 
ishing clumping (/ci= 1) in the outer wind zones of these objects could account 
for the observed Pa^ equivalent widths and radio fluxes of the stars. These authors 
used Pace, Br^, and the radio fluxes as diagnostics for the radial dependence of /d. 

A further st ep towards the characterization of wind clumping was taken by 
IPuls etall(l2006h . who were able to derive constraints on the radial stratification of 
fci by simultaneously modeling and the IR/mm/radio emission from a sample 
of 19 0-stars with well-known parameters. Actual constraints could be obtained in 
the (1 - 5 /?*) lower and intermediate wind region where Ha and the IR form, and 
in the outer (10 - 50 /?*) wind region where the mm/radio emission arises. Partic- 
ularly for the best constrained object, ^ Pup, the derived s tratification is in strong 
contrast to theoretical predictions (see Fig.|27ll. Actually. iPuls et all ( |2006[) found 
considerable clumping already close to the stellar surface for all denser winds, 
consistent with the findings from other diagnostics covered in previous sections. 
Clumping would then remain rather constant over a large volume before decreas- 
ing in the outer wind, with ratios of clumping factors between the inner and outer 
regions ranging from 3 to 6 (M(Ha) = 1.7 - 2.5 MrMm), which is very similar to 
the values suggested by Markova et al ('2004') and lRepolust et all (l2004l) on a com- 
pletely different basis (page|73]i. For weaker winds, on the other hand. IPuls et all 
(l2006h found similar clumping factors for the inner and outermost regions, i.e., 

M(H„) W A/radio- 



In a subsequent studv. lNugis a nd Lamers ( 2000) derived a n empir ical mass loss formula for 
Galactic WR stars accounting for these reduced rates, see Sect. l4.2.3l 
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Fig. 28 Examples for the radial structure of the volume tilling factor, /voi(r), for different values 
of CL2,3,4 according to Eq.|55] and CLj = 0.1. For models with an unclumped outer wind, CL4 = 
1, whilst for models with constant clumping in the outer part, CL4 = CLj . The velocity structure 
(solid magenta, in units of v^^) is displayed to illustrate the onset of the clumping variations. 



A major shortcoming of this study is caused by the dependence of all con- 
sidered diagnostics, which lead the authors to derive only relative clumping fac- 
tors. Note, however, that mass-loss rates based on the assumption that /ci(radio) 
= 1, i.e., an unclumped outer wind, would give rise to a unique WLR (depen- 
dent only on L), being in very good agreement with the theoretical predictions by 
IVink etall(l2000h . 

The next natural step is then to perform multi-wavelength analyses including 
bot h p and diagnost ics. These have been recently presented by Naiarro (2008) 
and lNajarro et all (l2008l) by means of UV to radio observations of OB stars and IR 



observations of LB Vs. These authors identify key diagnostics lines to obtain the 
clumping structure throughout the stellar wind, and provide guidelines to constrain 

the degree of clumping in stellar winds for different stellar types^ 

To investigate the clumping stratification in detail. Naiarro et all (l2008t) utilized 
the OSIf"*" star CyOB2#7 as a work bench and suggested the following clumping 
law (see Fig.| 



/™i (r) = CLi + ( 1 - CLi ) e + (CL4 - CLj ) e , (55) 

where CLi and CL4 are volume filling factors and CL2 and CL3 are velocity terms 
defining locations in the stellar wind where the clumping stratification changes. 
CLi sets the maximum degree of clumping reached in the stellar wind (provided 
CL4 > CLi) whilst CL2 determines the velocity of the onset of clumping. CL3 and 
CL4 control the clumping structure in the outer wind. This is illustrated in Fig. 1281 
which displays the behavior of clumping in a stellar wind for different values of 
CL2, CL3 and CL4. The clumping law Eq.|55]has been constructed in such a way 
that as the wind velocity approaches Voo, the volume filling factor starts to change 
from CLi towards CL4. If CL4 is set to unity, the wind will be unclumped in 
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the outer most region (as sugg ested bv lNugis et allll998l and adopted bv lFiger etail 
120021 and iNaiarro et alll2004 ' for the case of WR-winds, see ab ove). From Eg. 155] 
we no te that for CL3 ^ we recover the law as proposed by iHillier and Milled 
(Il999h . 

All clumping effects considered so far refer to a situation where the ionization 
fraction of the ion responsible for the considered process is not or only weakly af- 
fected by the clumping itself (at least beyond the conventional V/voi-scaling 
Eg. 154b. Thus, for a resonance line (e.g., P V) it must be secured that the ionization 
stage the line belongs to (P^"*) clearly remains as the dominant stage, so that by 
varying the clumping alone no changes in the line profile are produced. On the 
other hand, for a recombination line, the ionization stage the line belongs to must 
be clearly less populated than the next higher ionization stage (i.e., HlI^Hl for 
the hydrogen lines), so that an V/voi-scaling preserves the number of recom- 
binations. We consider the corresponding regime of the ionization equilibrium to 
be on the "safe" side. 

Most worrying is the case where two adjacent ionization stages have similar 
populations. Then clumping will cause a net reduction of the mean ionization, 
due to an enhanced number of recombinations (<== p^) compared to ionizations 
(roughly °= p). This will result in a line-profile behaviour similar to a reduction of 
the stellar temperature. We denote such a regime of the ionization equilibrium by 
the term "changing". This indirect effect of clumping has been warned about from 
early times on (e.g.. Puis et al 1993b), and Bouret et al (2005) have investigated 
this effect in detail. As an example. Fig. |29] shows how the oxygen ionization 
equilibrium is significantly altered due to wind clumping. 

One may, therefore, distinguish between lines formed in a "safe" region and 
those arising from a "changing" region. Figure [30l shows that the N IV /14058 and 
N V A4600 lines in CygOB2 #7 are clearly formed in such a "changing" region, 
whilst the rest of the optical lines are less affected by the large change of clump- 
ing and seem to follow the "safe" region behaviour. On the other hand, most of 
the infrared lines react strongly to changes in the clumping law. Since, as out- 
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Wavelength 

Fig. 30 Best model fits to UV, optical, and iiifrai-ed profiles of CygOB2 #7, for two values of 
CLi , denoted in the figure by / and corresponding to the maximum degree of clumping in the 
wind (= minimum volume filling factor). The weakly clumped model (dashed) with CLi,2,3,4 = 
(0.4, 350 kms^', 25 kms^', 1.) has a mass-loss rate of 6-10~^M(.-i yr~', whereas the s trongl y 
clumped model (dashed-dotted) with a clumping law according to IHillier and Milleri ( 119991) , 
CLi,2 = (0.025, 120 kms-i), has M = 1.4-lO-'*M0yr-'. 



lined on page|42] the continuum adjacent to NIR-lines is already formed in the 
wind, not only the lines, but also the continuum will depend on clumping (via 
bound-free and free-free processes p^), resulting in a high sensitivity of the 
continuum-re ctified line profiles on the clumping stratification, particularly on 
CLi and CL7. lNaiarro et all ( l2008l) found that if clumping starts with large clump- 
ing factors relatively close to the photosphere of CygOB2 #7 (with^j_=^5 reached 
at ^^200 km s^ ^ , corresponding to 1 . 1 similar to the results by lPuls et al 200^ 
for this object), consistent simultaneous fits to the UV, optical, and IR observations 
are possible (see Fig.[30]l. 

Whilst the optical and IR spectra of CygOB2 #7 provide strong constraints 
on CLi and CL2, the UV, submillimeter, and radio observations constitute crucial 
diagnostics to determine CLi and CL3. Inde ed, UV and s ubmillimeter data support 
the presence of constant clumping (Naiarr o et alll2008l) . at least up to mid/outer 
wind regions where t he millimeter cq ntinua of CygOB2 #7 are formed. However, 
radio observations bv lPuls et all (120061) with a formation region at much larger radii 
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showed that clumping may begin to vanish in the outermost wind regions, which 
is consistent with the fact that the radio emission from such models with constant 
clumping severely overestimates the upper limits provided by the observations by 
|Puls et al (2006) of CygOB2 #7. This demonstrates the need of multi-wavelength 
observations to constrain the run of the complete clumping stratification. 

Finally, for quantitative analyses one would wish to have a nice toolkit pro- 
viding diagnostic lines that, depending on the stellar type of the object, could be 
used to constrain the absolute degree of wind clumping. At temperatures around 
30 kK, there is the possibility to use He II 4686 and Ha in parallel, as outlined on 
page|5T] Otherwise, one has particularly to rely on the "changing" region situation 
discussed above. No strong UV lines should be included as they may be signifi- 
cantly affected by X-rays and corresponding EUV emission. Note, however, that 
even UV-lines formed in the lower wind could be affected by very hot material po- 
tentially located close to the wind base (see Sect. 16.5 1) . For early O supergiants, the 
N IV A4058 and N V A4600 lines, and for late O supergiants, the He II, N III and 
C III lines, turn into important clumping diagnostic lines. Valuable clum ping in- 
forma tion may be obtained from cool LB Vs for which the He I and Fe II] (iNajarrd 
l2Q08h lines will arise from "changing" regions. 



6.3.3 Treatment of clumping - Micro- and macro-clumping 

With studies yielding clumping factors ranging from /d = 2 to ^ 100 (correspond- 
ing to volume filling factors /yoi = 0.5 to ^ 0.01), one may question whether our 
present treatment of clumping is physically correct. If related to a certain type 
of instability (e.g., the line-driven one), we may expect this treatment to be in- 
adequate in those regions where the instability is not fully grown but still in its 
linear or only weakly non-linear phase. In these (lowermost) wind regions, the 
assumption of a void inter-clump medium is certainly questionable. Most of the 
atmospheric codes take into account only variations in the density, but the hy- 
drodynamic simulations also reveal strong changes of the velocity field inside 
the clumps, sometimes even changing the sign of the velocity gradient. Most im- 
portant, however, is the recent finding that the prerequisite of the standard micro- 
clumping approach, namely that the clumps are optically thin, needs to be relaxed. 

Within the micro-clumping approach (page l74l) . an optically thin clump is 
identified as a clump with a size smaller than the photon mean free path of the 
relevant matter-light interaction. In an optically thick clump ("macro-clumping"), 
photons will interact with matter many times before eventually being destroyed or 
scattered off the clump and escaping through the inter-clump matter. Of course, 
whether a clump is optically thin or thick will depend on the abundance, ion- 
ization fraction, and cross-section of the involved transition. Thus, we will have 
clumps that are optically thick for some processes while being optically thin for 
others. For a r ecent r eview of micro- vs. macro-clu mping, we refer the reader to 
iHamann et ail (l2008al see also lOskinova et alll2007l) . 

Micro-clump ing has been firstly implemented into NLTE model atmospheres 
by ISchmutj (fl995), and the corresponding formalism was outlined on page |74] 
For our discussion of macro-clumpin g we follow the s ame notation (see also 
lOwocki et aill2004l and the review by .Hamann et alll2008ah . 
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For optically thick clumps, photons do care about the distribution, the size and 
the geometry of the clump (see Sect. 16.5b they will encounter and multi-interact 
with. The conventional description of macro-clumping bases on a uniform clump 
size, Z(r), and an average separation among the statistically distributed clumps, 
L(r), which are related by the volume filling factor, 

/vol = (7)' = ^. (56) 

^ /cl 

Following the nomenclature of Eq. |53j the optical depth across a clump of size / 
and constant opacity Kc is given by 

Xc=Kcl= Kf^ll =^J2= ^^'^^ 

with mean opacity k (Eq. |53] | and porosity length h = I? jf', as introduced in 
Eq.|47] Remember that the porosity length constitutes a key parameter to quantify 
the characteristics of a structured medium and corresponds to the photon mean 
free path in a medium consisting of optically thick clumps (page[56]l. 

The effective cross section of the clumps, i.e., their geometrical cross section 
(surface) diminished by the fraction of transmitted radiation, is given by 

ac = /2(l-e-^c) (58) 

such that the effective opacity that needs to be used in the models results in 

/2(l-e-^c) (l-e-^c) (1-e-^c) 

V h Tc 

when ric is the number density of the clumps. Note that this equation holds for 
clumps of arbitrary optical thickness. When taking the optically thin limit, we im- 
mediately recover the micro-clumping approximation, jQff = k", which depends 
on the volume filling factor alone but not on clump sizes and distribution, whereas 
in the optically thick case, the effective opacity becomes strongly reduced, K^g = 
k/tc = h~^ and depends only on the porosity length, i.e., it is grey. Further com- 
ments on effects resulting from specific clump geometries are given in Sect. 16.51 

Effects on wind structure. iKrticka et all (12008') have studied the influence of 
micro- and macro-clumping on the predicted wind structure of 0-type stars, by 
artificially including clumping into their stationary wind models which have been 
described in Sect. 13. II 

If they assumed all clumps to be optically thin, the radiative line force in- 
creased compared to corresponding unclumped models, due to an increase of the 
electron density dependent 5-term in the force-multiplier (page|9ll, representing a 
changed ionization equilibrium with more lines from lower stages. This increase 
of the line-force leads to an increase of either the mass-loss rate or the terminal 
velocity, depending on the onset of wind -clumping (be low or above the critical 
point). Similar results have been found by lde Koter et ai 12008). 

If, on the other hand, the clumps were considered as optically thick. iKrticka et all 
(I2OO8) showed that a net reduction of the line-force may occur. Note, however, that 
this finding is in contrast to hydrodynamical simulation^, and originates from as- 

which do not find changes of the gross wind properties compared to a stationary solution, 
though the medium is clumpy (paeel71t. 
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Fig. 31 Porosity as a possible solution for the PV problem. Adapted from lOskinova et all i2007l) . 
see text. 



suming that the velocity law remains unchanged, in contrast to what would happen 
in a consistent approach. Finally, an inclusion of wind porosity into the continuum 
opacity (with the onset of porosity below the critical point) might also lead to a 
reduced M, due to an increase of the wind-ionization. 

In conclusion, a further, consistent investigation of the impact of wind-clumping 
(with its various facets) on the wind structure of massive O stars is crucial, since 
this might influence the predicted wind parameters, to which the observed ones 
have to be compared. 

Effects on line profiles. lOskinova et all (l2007h investigated the effects of poros- 
ity on the line profiles of 11, Pup, by means of the above formalism (Eq. [59] ) and 
assuming a "high" value for the mass-loss rate, M= 2.5 10~^M(T)y r~\ which is 
only mildly lower than the upper limit derived by IPuls etall (l2006l) . Using over- 
densities, /ci, inferred from diagnostic lines, and the equation of continuity to 
parameterize the radial dependence of the clump separation, L{r), they were able 
to account for (part o¥% the radiative transfer effects of clumping. 

In particular, only changes in the "formal integral" were considered (by using 
the effective opacity), whereas the feedback of macro-clumping on the NLTE- 
populations was neglected by arguing that most transition rates are not affected. 
From Fig. [3T1 it becomes clear that the most important effect regards strong res- 
onance lines like P V AA 11 18-1 128, which can be reproduced by this approach 
without relying on (very) low mass-loss rates, due to the reduction of the effective 
opacity when the clumps become optically thick. On the other hand. Ha is not 
affected by wind porosity, remaining optically thin inside the clumps and being 
reproduced in parallel with P V. Insofar, these results might open a way out to the 
"P V problem". 

Clumping from the other side: "Vorosity". As indicated above, line transitions 
should be strongly affected by a strongly non-monotonic velocity field as result- 
ing from the line-driven instability. Aiming to examine the potential role of such a 

Note that the original porosity concept had been developed for continuum processes; line 
processes are additionally affected by changes in the velocity field, which are simulated here by 
adopting a velocity dispersion inside the clumps (cf. page[36t, but see below. 
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Fig. 32 Vorosity. Multiply non-monotonic velocity structure arising in a 1-D hydrodynamic 
simulation with self-excited line-driven instability, plotted versus a mass-coordinate (left), and 
schematic simplification (right). The smooth wind flow is represented by the straight line. 
Adapted from Ow ocki ( 2008), see text. 



flow structure in reducing the observed strength of wind absorption lines. lOwockil 
f2008) performed dynamical simulations of the line-driven instability in which the 
line strength can be described in terms of a "velocity clumping" factor. Interest- 
ingly, this velocity clumping factor is insensitive to spatial scales. 

The left panel of Fig. 1321 shows the velocity structure arising in a 1-D simu- 
lation. The intrinsic instability of line-driving leads to substantial velocity struc- 
ture, with narrow peaks corresponding to spatially extended but tenuous regions 
of high-speed flow. Such regions are followed by dense, spatially narrow clumps/ 
shells which, when plotted as a function of mass-coordinate, give rise to extended 
velocity plateaus. The right panel illustrates the situation by means of a simplified 
model, where the velocity clumping is represented by a simple staircase structure, 
quantified by a velocity clumping factor. This value is set by the ratio between the 
internal velocity width 8v to th e veloc ity separation Av of the clumps. 

Using this model, lOwockil ( |2008|) computed dynamic spectra directly from 
instability simulations, which indeed exhibited a net reduction in absorption, due 
to the porosity in velocity space, an effect he denoted by "vorosity". The key 
problem in explaining a "reduced" P V line-strength in terms of vorosity is that 
one needs to have a large number of big velocity gaps, which is not what comes 
"naturally" from the instability simulations (Fig. [32]- left). Thus, only a modest 
reduction at a 10-20% level was found, which is well short of the factor 10 (or 
more) required in terms of the P V problem. 

Obviously, it remains a future challenge to investigate scenarios including both 
porosity and vorosity and their possible inter-relation. 

6.3.4 Evolutionary limits 

If mass-loss rates would need to be revised drastically (e.g., by factors of ^ 10 or 
more), this would have dramatic consequences for the evolution of and feedback 
from massive stars. Recalling a corollary statement in massive star evolution, that 
a change of the mass-loss rates o f massive stars by even a factor of two has a 
dramatic effect on their evolution dMeynet et all 1 1994^ . one immediately realizes 
the importance of quantifying the wind clumping phenomenon. This statement 
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illuminates a potential alternative to constrain the "allowed" range of mass-loss 
red uctions, namel y by evolutionary considerations themselves. 

iHirschil (120081) has recently explored the evolutionary implications. Taking for 
example a 120 Mq star with a lifetime of 2.5 million years and an average 
2.5 • lO~^M0yr~\ the star will have lost arou nd 50 Mfv^ at the end o f the main 
sequence using the theoretical prescriptions of I Vink et alll2000l 1200 ll whilst the 
star would lose only 5 Mq if the mass-loss rate was reduced by a factor of 10. 
The two big questions regarding evolution are then how to produce WR-stars with 
such low mass-loss rates and whether one would end up with too many critically 
rotating stars. 

Concerning the WR issue, one way out would be to shift the bulk of mass 
loss on other evolu tionary phases such as the LBV (e.g. Nugis and Lamers 200(1 
ISmith and OwockillMjaD and the red supergiant (RSG) stages. However, although 
mass-loss rates are more challenging to determine in these two stages, it seems 
unlikely that they could cope with the bulk of it. 

A seco nd, provocative , possibility is to place all massive stars in close binary 
systems CKobulnickv and Frver.2007,) . This scenario, however, would not only fail 
to produce the many RSG stars observed, but there is also observational evidence 
against all WR-stars being in close binary systems (only 30-40% in the Magel- 
lanic Clouds, Foellmi et al 2003a.b). Another possibility is to assume that WRs 
are produced by strong rotational mixing rather than by mass loss. Since fast rota- 
tors constitute only a small fraction of the whole population, this scenario seems 
very unlikely though. 

Combining the above points. IHirschil 120081) concluded that evolutionary mod- 
els could s urvive with M re ductions of at most a factor of ^2 in comparison to the 
rates from lVink et all (l2000l) . whilst factors around 10 are strongly disfavored. We 
note that a reduction of the Vink et a 1 (|2000L 120011) rates by a factor of two may 
correspond to an "allowed" reduction of the empirical mass-loss rates of a factor 
of about four - at most. 



6.4 Co-rotating interaction regions (CIRs) 

In the previous section, we have extensively discussed features attributed to micro- 
structure in the wind, expressed in terms of micro- and macro-clumping, and con- 
sequences of this. However, there is also large-scale structure present, inferred 
most directly from the so-called discrete absorption components (DACs). 

Discrete absorption components. These are optical depth enhancements ob- 
served in the absorption troughs of unsaturated UV P Cygni profiles and are 
present in most 0-/early B-star winds (Ho warth and Prinjalll98 9V but also in late 
B-superg iants (Bates and Gilheanv 1990) and at least one WN7 star (Prinja and S niidil 
[l992fFr 

Typically, these DACs accelerate to the blue wing of the profile on ti mescales 

of a fe w da ys, becoming narro wer as they approach an asymptotic velocitv. lHowarth and Prinja! 
( il989h and jPrinja et all ( Il990l) have used this asymptotic value to determine Voo for 

Absorption components have been also observed in the Balmer lines of the LBV P Cygni 
dMarkovalligsg) . 
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objects with no black trough in CiV (see Sect. 14. it . The acceleration and recur- 
rence time scales of DACs were found to be correlated with v sin /, i.e., with the 
rotational period (|Prinialll988llHenrichs et alll 19881 IXaper et alUlQgg l). The accel- 
eration of most strong DACs is much slower than the mean wind acceleration 
(IPrinja et allll"992L IPrinjalll994h . suggesting that DACs might arise from a slowly 
evolving perturbation through which the wind material flows. Such a perturba- 
tion could consist of a higher dens ity or a lower vel ocity gradient (e.g., a velocity 
plateau, as speculated alr eady by Lame rs et all 19821) . or result from a combination 
of both dFullerton and Owo cki 1992,). 

Though the lack of variability in the emission compo nent of P Cygn i lines 
(|Prinjall992 ') and the lack of significant infrared variability ('H owarthlll99 2') seems 
to rule out a spherically-symmetric disturbance, the structure must be large enough 
to cover a substantial fraction of the stellar disk in order to produce the observed 
strong absorption features. The latter constraint strongly argues against any micro- 
structure related to the wind instability as the source of the DACs, because (i) 
micro -structure would be averaged out and (ii) the presence of the lateral line- 
dras (iRvbicki etallll990h should suppress any lateral velocity disturbances and 
maintain lateral length scales as set by a potential base perturbation. 

Dynamical models of CIRs. The above arguments suggest that DACs should 
originate from coherent structures of significant lateral extent, with a clock related 
to stellar rotation and an increased optical line depth due to density and/or veloc- 
ity field effects. Already in 1984, i.e., well before these constraints were known, 
Mullan had sug gested c o-rotating interaction region^ (CIRs) as a potential can- 
didate (see also lMullan|[l986l) . Over the following years, this suggestion became 
more and more likely, due to an amount of empirical knowledge and since CIRs 
are compatible with all constrain ts as outlined above. 

ICranmer and Owockil (Il99 (j) were the first to investigate the CIR scenario in 
detail, by means of 2-D time-dependent hydrodynamic modeling of the wind of 
a rotating 0-star. Photospheric disturbances in the form of azimuthal variations 
were induced by a local increase of the line forcS due to a bright stellar spot 
in the equatorial plane. Such a stellar spot should be regarded as a representative 
for various kinds of photospheric disturbances (e.g., localized magnetic fields or 
non-radial pulsations), since the induced wind structure turned out to be rather 
insensitive to the particular way by which the photospheric conditions were dis- 
turbed. 

A locally enhanced line acceleration generates streams of higher density and 
lower velocity (compared to the undisturbed mean flow), and CIRs of enhanced 
density form where the faster mean flow collides supersonically with the slow 
material at the trailing edge of the stream. When the unperturbed wind collides 
with the CIR, a non-linear signal is sent back towa rds the star which forms a sharp , 
propa gating discontinuity ("Abbott-kink", see also lFeldmeier and Shlosmanl2000L 
l2002h in the radial velocity gradient. The inwards directed propagation velocity 
(with respect to a comoving frame) of this kink is just the characteristic speed of 

well known and studied in the solar wind. 

Actually, the authors investigated also the consequences of a local decrease in the line force 
due to dark spots. Such a model, however, could be ruled out, since the resulting synthetic 
profiles behave quite different from the observed variability pattern. 
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the radiative-acoustic Abbott waves considered in Sect. 16.21 In the stellar frame, 
this feature travels slowly outwards and a velocity plateau-like structure is formed 
between the trailing Abbott-kink and the CIR compression. 

By calculating the corres ponding line optical depths and time-dependent ("dy- 
namic") synthetic profiles. ICranmer and Owockil (Il996h showed that these slowly 
moving kinks, together with their low velocity gradients (T °= (dv/dr)~M), are the 
likely origin of the observed DACs rather than the CIRs themselves (but see be- 
low). A fit to the temporal wavelength drift of the synthesized DACs resulted in a 
velocity field parameter j3 w 2 — 4, i.e., in a slow apparent acceleration, in accor- 
dance with observations. Though the recurrence time scale of the DACs is strictly 
correlated with Vrot, due to the link of CIRs and the rotating photosphere, no cor- 
relation between acceleration time scale and rotational speed was found in the 
synthetic profiles, in contrast to the observational indications as outlined above. 

This "missing" correlation is the central topic of an interesting kinematical in- 
vestigation by Hamann et al (2001) , who stressed the fact that the low drift rate is 
not a consequence of the CIRs themselves. Features formed within the CIR would 
display a drift even faster than features formed in a non-rotating wind. Instead, the 
slow drift is a consequence of the difference between mean flow and the velocity 
field of the pattern in which the f eatures form. A pattern wit h upstream propaga- 
tion (as the kink in the model bv iCranmer and Owocki|[l995) inevitably results in 
a wavelength drift with a slower apparent acceleration than displayed by features 
formed in the bulk or within the CIR itself. Any of these drifts, however, were 
shown to be independent on the rotation rate, leaving the observed correlation (if 
actually present) still unexplained. 

An implicit assumption in the work by iHamann et all (l200lh is that the CIRs 
are induc ed by disturbances lock ed to the stellar surface, in accordance with the 
model by ICranmer and Owockil ([l996). In a detailed study, iLobel and Blomm3 
(|2008") relaxed this assumption in order to enable a quantitative analysis of the 
time evolution of DACs for the fast-rotating BO. 5 lb supergiant HD 64760, one 
of the best observed objects (see below). The authors use their own 3-D radiative 
transfer code, coupled with 3-D hydrodynamics restricted to the equatorial plane. 
A large grid of models and dynamic spectra for different spot parameters (bright- 
ness, opening angle and velocity) has been computed, and the effects of these 
parameters on the wind structure and DAC evolution studied. A best fit could be 
obtained with a model with two spots of unequal brightness and size on opposite 
sides of the equato r, with spot velocities being five times slower than Vrot- All ba- 
sic conclusions of ICranmer and Owockil ( IT996) could be confirmed, particularly 
the importance of velocity plateaus in between kinks and CIRs. For their models 
of spots which are not locked to the stellar surface, however, the independence of 
the DAC acceleration on Viot can no longer be warranted. In this case, the internal 
clock is given by the rotational speed of the spots, Vjp, whereas a change in the 
stellar rotational speed modifies the underlying bulk flow (densities and veloci- 
ties) via a different centrifugal acceleration (Eqs. l20ll21b such that the time-scales 
are no longer conserved when Vrot is changed. In other words, two models with 
Vsp = Viot (locked) and Vsp < Vrot (non-locked) which both reproduce the DAC 
recurrence time behave differently with respect to the acceleration time scale. 

The physic al scenario of different spot and stellar rotation velocities follows 
a suggestion bv lKauferetal (120061 These authors detected non radial pulsations 
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in the photospheric lines of HD 64760, with closely spaced periods. In order to 
explain the observed Ha variability, they invoked the beat period between two of 
these periods (lower than the rotational one!) to be responsible for the CIRs in the 
wind. Interestingly, howeve r, this beat period is not con sistent with the spot veloc- 
ity derived from the UV by iLobel and Blommd (l2008l) . For further discussion we 
refer to the latter publication. 

Reality or artefact? As should be clear from our discussion of the onset of 
Abbott-waves, the dynamic models discussed above rely on a Sobolev line-force, 
because of computational feasibility. Thus, two questions need to be answered. 

(i) In how far do the kinks still occur when the Sobolev approximation is relaxed? 

(ii) In how far can the line-driven instability disrupt the CIR structure? A prelim- 
inary answer might be found in Fig. 2c from lOwockil (Il999f) . which displays the 
radial velocity calculated for a CIR model based on a non-local line force cal- 
culated from a three-ray treatment (cf. page 1701) . Also here, somewhat "rounded" 
kinks followed by a plateau are present ! Unfortunately, however, the line-driven 
instability effectively destroys all macro-structure, but this might be related to a 
still inadequate three-ray treatment and/or insufficient resolution, leading to too 
low a lateral damping. To cite Stan Owocki (2008, priv. comm.), "there is a strong 
potential that such kinks and plateaus are real physical effects in a line-driven 
wind, and not just an artefact of the Sobolev approximation." 

Future work using computationally expensive methods to calculate a consis- 
tent more-D line force will hopefully prove this expectation. 



Modulation features. The observed correlat ions of the DAC properties with 
V sin i inspired the "WE MEGA Campaign " (iMassa et all Il995h . during which 
three prototypical massive stars (i^ Pup (0411), HD 64760 (BO.SIb) and HD 50896 
(WN5)) were monitored almost c ontinuous l y over 16 days. Results of this cam- 
paign have been summarized by Fullerton' (1999). The most important finding 
was the detection of a new type of variability, namely periodic modulations in the 
UV wind lines, derived from the dynamic spectra of HD 6 4760 dPrinia et all 19951) 
which most likely is observed equator-on. As shown by iFullertoii et all (119971) . 
these modulations result from two quasi-sinusoidal fluctuations with periods of 
1 .2 and 2.4 days, which propagate, as a function of phase and starting at roughly 
750 kms~^, simultaneously towards lower and higher line-of-sight velocities, un- 
til a certain minimum and maximum is reached and the pattern begins to repeat 
itself. Because of its peculiar shape, this effect has been called "phase bowing". 
Similar features (th ough not as pronounced) have be en detected in only two other 
stars, C Pup (Mass a et alii 19951 iHowarth et allll995h and ^ Per (07.5 III (n)((f)), 
iKaper et allll999() . Note that these features occur in parallel with the conventional 
DACs. 

The peculiar phase properties of the modulation features indicate the presence 
of absorbing material with the same phase at two different projected velocities. By 
means of a kinematical rnodel an d corresponding synthetic spectra, Owocki ejLaJ 
([19951) and lFullerton et all d 19971) suggested that these features are formed in az- 
imuthally extended, co-rotating and spirally-shaped wind structures linked to a 
surface density which is modulated by non-radial pulsations. Due to rotation, a 
certain spiral exits from the lines of sight in front of the stellar disk first at an 
intermediate projected velocity (750 kms~4n case of HD 64760). At later times 
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and because of its curvature, this structure exits from the absorption region si- 
muhaneously at lower and higher velocities. Since the spirals are assumed to be 
linked to the surface density modulation, different locations along the spiral have 
the same phase, because they arise from the same "spot" on the surface. In concert 
with the kinematics of the spirals, this explains the modulation features and their 
phase-bowing, verified in the aforementioned synthetic spectra. Consistent hydro- 
simulations of this scenario are still missing, but these features are a strong hint 
on the presence of azimuthally extended macro- structure and obse rvational "ev- 
idenc e for co-rotating wind streams rooted in surface variations" dOwocki et all 
[HH). 

6.5 X-ray line emission 

After the advent of the X-ray satellites xmm-Newton and CHANDRA, with their 
capability to perform high-resolution spectroscopy, detailed knowledge about the 
X-ray line emission fro m hot stars ha s been accumulated (for recent reviews, see 
lOskinova et alll2008l and ICoherill2008l) . Two "classes" of spectra have been found. 
On the one side, the prototypical spectrum of d ' Ori C, one of few massive stars 
with a strong magnetic field (Sect. \2.2.2i . is rather hard (implying high plasma 
temperatures) and shows lines with small widths. On the other side, the spectrum 
of the prototypical "normal" supergiant ^ Pup is softer and displays much broader, 
asymmetric lines. Moreover, the lines from hydrogen-like ions are stronger in 
the former class, and vice versa for the lines from helium-like ions, which im- 
mediately points to differen ces in the ionization ba l ance a nd thus again to dif- 
ferent plasma temperatures. jWojdowski and Schulzl (120051) derive a peak value 
r = 30 • 10* K for 0i Ori C and 1-2- 10* K for six other 0-stars (incl. ^ Pup), 
where the latter values are of t he same order as tho se derived from earlier, low 
resolution data (e.g., by ROSAT. lKudritzki et all! 1 9961 see Sect.lej]). 

6.5.7 Magnetically confined winds 

By means of a tailored MHD modeQ for the wind of 6 ^ Ori C. lGagne et al ( l2005l) 
showed that the magnetically confined wind shock (MCWS) model (discussed in 
Sect. l2.2.2] l fits quite well all the data, predicting the temperature, luminosity, and 
occultation of the X-ray emitting plasma with rotation phase. Particularly, the bulk 
of the shock-heated plasma should be located in the magnetically confined region 
close to r w 2/?*, where the speed and line-of-sight velocity of this material is low, 
giving rise to rather narrow X-ray emission lines. Similar arguments might hold 
for the X-ray spectrum of T Sco which also shows narrow lines, consistent with 
the finding that also T Sco has a strong magnetic field. 

6.5.2 Non-magnetic stars 

In accordance with previous results derived from lower resolution observations 
(Sect. 16. Il l, also the gross properties of highly resolved X-ray spectra from "nor- 
mal" 0-stars may be well explained by the line-driven instability scenario (either 



extending the work bv lud-Doula and Owockil l i2002h with respect to the energy equation. 
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Fig. 33 Schematic Grotrian diagram for Helium-like ions, showing levels and transitions in- 
volved in the formation of the fir-complex in X-ray spectra of hot stars. Solid arrows: collisional 
excitations. Dashed arrows; radiative transitions. The f,i,r transitions are denoted by correspond- 
ing letters. The transition s betw een the ^S and the levels are located in the (F)UV. Adapted 
from iGabriel and Jordai] ^9691) . 



self-excited or triggered, with dominating clump-clump collisions). This refers 
particularly to the softness of the spectra and the large line widths, resulting from 
the high velocity of the shock-heated wind. Remember, e.g., that the line-driven 
instability should be fully grown above 1.3-1.5 R^,and that the predicted shock 
temperatures are of the order of "only" a few million Kelvin, due to jump veloci- 
ties of a few hundred km s^ ' . 

fir-diagnostics of Helium-like ions. A key diagnostics to observationally con- 
strain the onset of the X-ray emission is provided by the ratio of line fluxes from 
Helium-like ions s uch as Ovil, NelX, MgXI and SiXIII. These ions show charac- 
teristic^?- triplets taabriel and Jordarilll969h of a forbidden (f), an intercombina- 
tion (i) and a resonance (r) component (cf. Fig. [33]). 

In hot stars, the M = f ji line ratio is sensitive to the local mean intensity 
of the UV radiative field, since UV photons can excite the metastable ^Si level 
to the ^Po,i,2 levels. To a good approximation, this mean intensity is given by 
W(r)//v'^°', with dilution factor and photospheric flux H^°^ . If the latter 
is known (from observation^ or model atmospheres), the dilution factor W(r) 
may be inferred from !M, which constrains the formation radius of the X-ray 
emitting plasma. Applying this method, iLeutenegger et all (l2006 m studied the X- 
ray lines of Helium-like ions in four evolved 0-stars (incl. t, Pup) and derived 
minimum radii of X-ray line formation in the range of 1.25/?* < /?fir < 1-67/?*, 
consistent with the the oretical expectation. 

On the other hand, IWaldron and Cassinellil (l2007t) published rather puzzling 
results for a larger sample of OB-stars. Based on similar yrr diagnostics, they con- 
cluded that high-Z ions were predominantly located closer to the photosphere, a 

if the transition is located in the observable part of the UV spectrum, and not contaminated 
by wind-lines. 

see also Leutenegger et al 1 2007) for the effects of resonance scattering in the r-component, 
which give rise to more symmetric fir-complexes. 
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Fig. 34 Two different models that fit the Fexvil line in the Chandra spectrum of ^ Pup. 
Dashed histogram: non-porous model with t, = 2. Solid histogram: porous model with t, = 8. 

= 3.3. The dashed vertical line refers to the transition frequency at rest and the dotted lines 
to frequencies corresponding to ±v„o. See text. (From D. Cohen, priv. comm.) 

problem which they named the "near-star high-ion" problem (see also lWaldron and Cassinellil 
I2O O 1 ) . The derived conditions (temperatures 1 -2 • 1 0^ K close to the star, with /?flr 
as low as 1.1 would be in stark contrast to the standard shocked- wind scenario, 
since, e.g., the implied jump velocities would be higher than the velocity of the un- 
derlying wind. Moreover, and assuming that the average wind velocity increases 
with the radius, lines from lower Z ions should be broader than from ions with 
higher Z, at least in terms of the standard model of X-ray emission (see below) Fl 
Contrad ictory to this expectation, the observed widths of all X-ray lines are simi- 
lar (e .g., rWaldron and Cassinellil200lL [Kramer et aill2003LIWaldron and Cassinellil 
l2007h . More investigations are certainly needed to clarify the situation. 

Analysis of emission line profiles. Further clues on the X-ray emission from hot 
stars are provided from the profile shape of strong, single lines (mostly hydrogen- 
like, e.g., Ovill, Nex, MgXll, SiXlV, but also FeXVII), which is asymmetric and 
skew ed. Note first that such a line shape is consistent with the "standard model" 
(e.g., lHillieretall [l993h which assumes the presence of hot, line-emitting mate- 
rial mixed with a "warm", continuum-absorbing mean wind flow. If the wind 
is optically thick in the continuum, the emission from the receding hemisphere 
will be significantly attenuated, i.e., a deficit of red-shifted photons will be ob- 
served, whereas from the front, less attenuated side there will be comparatively 
more emission, contributing to the blue-shifted part. Taken together, an asym met- 
ric profile should be created, already predicted by iMacfarlane et all (Il99lh and 
indeed observed, as shown in Fig. |34] 

Unfortunately, however, this principal agreement was severely challenged after 
quantitative investigations had been performed. Already from the first analyses on 



ICassinelli et all <2008t) speculate that the "near-star high-ion" problem might be explained 
by infalling clumps and corresponding bow shocks. 
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dOwocki and Cohenll200lLlKramer et aJ2003h it turned out that the red parts of the 
profiles were much less attenuated by the underlying wind than expected. Even 
worse, the continuum opacity, Ky, increases with wavelength, and consequently 
the lines of lower Z ions are expected to be more skewed compared to the higher Z 
ones. Such a dependence, however, has not been found, and the wavelerigth sh ifts 
of line centroids are similar for lines of all ions tWaldron and Cassinellill2007h . 

Lower mass-loss rates or porosity/macro-clumping ? Different hy potheses to ex- 
plain these problems have been proposed and are summarized by ' Oskinova et all 
(12008 ). In the following, we will concentrate on two promising alternatives that at 
least allow for reasonable profile fits, but are still subj ect to a lively debate (see, 
e.g., the discussions recorded in iHamann et alll2008b h. Both approaches rely on 
the notion that a reduced continuum opacity can lead to a better agreement be- 
tween modeled and observed line shapes. This can be achieved either by reducing 
the adopted mass-loss rate or by accounting for porosity effects. The latter have 
the additional advantage that the effective opacity, in the optically thick limit, be- 
comes fr equ ency-independent, i.e., grey (Eg. 1471) . thus potentially explaining the 
similarity^ of the line-shapes from different ionization stages. 

As a typical example, Fig.[34]displays the observed Fexvil pro file from C Pu p 
and two fits. The dashed histogram refers to the b est fit obtained bv' CohenI (|2008|) . 
based on the empirical model bv lOwocki and C ohen (2001), which optimizes 
three parameters: the normalization, the inner radius below which there shall be 
no emission (Ro), and the continuum optical depth of the "warm" wind, parame- 
terized by the quantity T* = KyM / {4nR^v„). The fit is formally good, and the best 
fit values (within 68% confidence limits) are Ro = 1.53;°;}^ and T* =2.0 ±0.4. 

Thus, the Rq derived from the profile is consistent with the expected on- 
set of shock formation and the optical depth is quite small indeed. Translated 
to a mass-loss rate, this would imply M = 1.5 •lO~^A/0yr"', which is (much) 
lower than recent estim ates from NLTE an alyses based on unclumped models 
(Mw 8.8 lO-^M^yr-'. lRepolust etaill2004h . and a factor of 2.8 lower than the 
maximum mass-loss rate {M = 4.2 •lO~^A/0yr~') constrained by independen t p^- 
diagnostics when the outer wind is assumed to be homogeneous dPuls et al.2006l 
see Sect. l6.3b . 

Alternative models were simultaneously presented by lOskinova et all (120061) 
and[Owocki and Cohen (2006). Both publications (which differ with respect to 
important details) suggested that also porosity associated with macro-clumping 
(Sect. 16.3.31) can account for the low degree of asymmetry, without any need to 
worry about red uced mass-loss rate s . The solid histogram in Fig. [34] bases on 
the approach by lOwocki and Cohen! (|2006) and displays the best fit for such a 
porous model, assuming spherical clumps. By fixing the mass-loss rate at M= 
6 •lO~^M0yr~^ (corresponding to T* = 8), the porosity length, h (Eq.|47ll is a free 
fit-parameter in this model. It is evident that the best fitting porous model with the 
"high" mass-loss rate is almost indistinguishable from the best fitting non-porous 
model with the low M. The optimum terminal porosity length (the value of h in 
the outer wind region) is rather large, 2.57?,, < < 4R^ at a 68% confidence 
level. Cohen (2008) argued that this is much larger than any porosity length seen 



The issue of the actual degree of grayness has not been settled though. 
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in state-of-the-art hydrodynamical simulations, and concluded that "there is no 
co mpelling evidence" fo r invoking porosity rather than a reduced mass-loss rate. 

lOskinova et all ( l2006h pointed out that by combining the fits for t* from ^ Pup, 
based on non-porous models (from .Kramer et al 2003) with appropriate absorp- 
tion coefficients, different l ines would imply signi fican tly different mass-loss rates. 
Based on previous work bv lFeldmeier et all ( l2003h and lOskinova et all(l2004l) . they 
describe the wind by a flow of clumps obeying the equation of continuity. As in 
the porosity length formalism, the effective opacity of the clumps becomes grey 
when the clumps become optically thick. They assume that the emission originates 
between two radii ri and r2 (derived fromyzr-lines), and investigate the influence 
of the clump geometry. Comparing spherical clumps with (infinitesimally) thin 
shell fragments (pancakes) oriented perpendicular to the radial direction, they ar- 
gue that the latter geometry leads to more symmetric profiles than the former one, 
as long as the clumps subtend a constant solid angle as they propagate outwards. 
In their subsequent simulations of synthetic X-ray lines, the shape of the clumps 
is simulated by an anisotropic opacity (per default assumed to be pancake-like), 
and the only parameter to be specified (in addition to the wind parameters, M, Voo, 
and j3, adopted from the literature) is the quantity no = n{r)v{r). This quantity 
defines the number of fragments passing through some reference radius per unit 
time (constant due to mass conservation), and was set to mq = v„/R*, in order 
to obtain average clump separations L(r) which are compatible with predictions 
from hydrodynamic simulations, L{r) = /?*v(r) /voo < ^m.. This model has no free 
parameter, but nevertheless provides a satisfactory agreement between a variety 
of synthetic and obser ved lines from d ifferent ions (thus supporting the grey char- 
acter of the opacity). lOskinova et all (12008.) conclude that these results provide 
evidence that the wind clumps are not optically thin, that they are compressed in 
the radial direction, and that they are separated by a few tenths of the stellar radius 
in the wind acceleration zone. 

Let us finish this section with our conclusions. Whilst the X-ray spectra from 
magnetic winds seem to confirm the magnetically channeled wind shock scenario, 
the situation for non-magnetic winds remains unclear. Both non-porous models 
with reduced mass-loss rates and porous winds with "conventional" M are able to 
explain the observed spectra, where the role of non-isotropic opacity needs further 
confirmation by independent diagnostics. Future modeling with a simultaneous 
analysis of processes from the X-ray to the radio domain will certainly uncover 
the most appropriate description and the "real" mass-loss rates and clarify the 
"near-star high-ion" problem. Such modeling should include not only the X-ray 
line-features, but also the SED in the complete X-ray band including the wealth 
of the ROSAT data. 



7 Summary and further implications 

We reviewed the current status of both theoretical and observational aspects of 
mass loss from hot massive stars, starting with the standard model and the wind- 
momentum luminosity relation (WLR). Mass-loss predictions from line statistics 
were discussed with particular emphasis on the temperature and metallicity de- 
pendence of the wind momenta, emphasizing that the spectral-type dependence 
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of the WLR, especially with respect to the bi-stability jump, is still a most active 
aspect of present-day hot-star wind research. 

We continued with the impact of rotation. Wind compression and the predic- 
tions of predominantly prolate winds due to gravity darkening were reviewed, 
whilst we concluded that the rotationally-induced bi-stability concept for explain- 
ing the equatorially outflowing winds from B[e] supergiants requires future 2D 
modeling. With respect to magnetic fields, we made a clear distinction between 
internal and surface magnetic fields and note that firm observational constraints 
on the strength and geometry of the fields have become within reach of current 
instrumentation (e.g., using ESPADONS or NARVAL). We reviewed the wind struc- 
ture as a function of the magnetic confinement parameter, and pointed out that 
weak magnetic fields below present detection limits could have a large impact on 
weak winds. It remains a possibility that the co-rotating interaction regions that 
are the likely root of ultraviolet line profile variability of 0-stars might be related 
to an anchored field on the surface. 

A key topic concerned stationary models of radiation-driven winds. We de- 
scribed two basic methods in use to predicting the mass-loss rates of massive stars 
and their pros and cons: hydrodynamical methods based on the improved CAK 
approach and those based on Monte Carlo radiative transfer methods. Predic- 
tions were provided for O and B supergiants, followed by those for stars in closer 
proximity to the Eddington limit, including Luminous Blue Variables (LBVs) and 
Wolf-Rayet (WR) stars. Furthermore, we provided some basic insights into the 
possible existence of stellar winds at very low and zero metallicity content (Popu- 
lation III) involving issues such as the potential self-enrichment of metals through 
the evolution of the first generations of massive stars themselves, and the potential 
importance of continuum driven winds from super-Eddington stars, which do not 
require the presence of metals in their atmospheres. Particularly, we considered 
the predicted metallicity dependence of WR winds since this has relevant impli- 
cations for the angular momentum evolution of massive stars at low metal content 
and the progenitor stars of long GRBs. 

As far as observational wind diagnostics are concerned, we discussed the tradi- 
tional UV, Ha, and continuum bound-free/free-free emission, in addition to mod- 
ern non-LTE atmosphere analyses and emphasized the future potential for NIR 
spectroscopy. One of the major steps in the last decade has been the inclusion of 
metal line-blanketing, which has led to a new temperature scale of OB-stars, with 
generally lower temperatures and luminosities. Due to these lower luminosities 
the 0-star mass-loss predictions per object have been down-revised as a conse- 
quence of the changes in underlying stellar parameters - without actually chang- 
ing the predictions themselves. Empirical mass-loss rates have been determined 
in parallel with stellar parameters, by profile fitting methods (using the conven- 
tional "by eye" method or optimization via genetic algorithms), for a variety of 
spectral types, luminosity classes and environments, particularly within the VLT- 
FLAMES survey of massive stars. In many cases, the derived mass-loss rates of 
0-supergiants were a factor of ~2 above predicted values. This could either mean 
that the theoretical values are not large enough, or that the empirical values are too 
high (or that both need revision). For mid and late type B-supergiants, on the other 
hand, empirical M-values were found to be lower than the theoretical ones. Com- 
paring the mass-loss rates for the Galaxy, the LMC and the SMC, and accounting 
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for an average clumping correction, the empirical dependence of O-star mass-loss 
rates on metallicity has been derived, which agrees quite well with the theoretical 
predictions. 

In addition to the minor problems outlined in the last paragraph, our general 
understanding of radiation driven winds is challenged by two urgent problems, 
the weak wind and the clumping problem. Regarding the former, we have dis- 
cussed that the presence of X-rays can lead to a degeneracy of mass-loss rates 
derived from the UV, and suggested the NIR Bra-line as a promising tool for fur- 
ther progress on quantifying mass-loss rates for thin winds. We have also outlined 
that just these X-rays might be the origin of the weak wind problem, by modifying 
the ionization equilibrium and thus the efficiency of radiative driving in the lower 
wind. Future work has to provide tight constraints on the domain of occurrence 
of weak winds, and particularly on the question which parameters discriminate 
objects with "normal" from those with weak winds. To be even more provocative: 
Can we be sure that there are any late O-dwarfs with "normal" winds at all, or 
did we over-interpret present mass-loss diagnostics, and the weak wind case is the 
normal one? 

In the last part of this review we discussed important issues related to time- 
dependence and structure, featuring the clumping problem. We reported that dif- 
ferent diagnostics seem to suggest also different clumping properties (with clump- 
ing factors differing by factors up to 20-100), and summarized recent attempts to 
unify these results by invoking optical depth effects (micro- vs. macro-clumping) 
and clumping in velocity space. One particular aspect concerns the radial stratifi- 
cation of the clumping. Most analyses seem to have at least one result in common, 
namely the indication of considerable clumping close to the wind-base. This is in 
strong contrast to the conventional interpretation of wind-clumping as caused by 
the line-driven instability (either self-excited or triggered), since all correspond- 
ing simulations predict the onset of wind-clumping only at rather large velocities 
(some hundreds of kms~^). Thus, future effort should be dedicated to improve our 
understanding of the physical cause of the clumping, but also to continue with em- 
pirical studies, by combining evidences from different diagnostics that are formed 
in different positions in the wind. We have provided a promising "toolkit" for 
this objective. Of course, the final goal is not only to "measure" clumping/ volume 
filling factors and porosity lengths, but also to accumulate information about the 
shape of the clumps. Most spectral diagnostics (UV, H^, Br^, IR-/radio contin- 
uum) do not provide such geometric constraints, but we have pointed out that 
linear polarimetry could help, as well as the diagnostics of X-ray line emission. 

Given the overwhelming observational evidence for wind clumping, it is ev- 
ident that those empirical rates which are based on diagnostics and homoge- 
neous wind models yield maximum values and that the real mass-loss rates must 
be lower. Still, the most relevant question is by how much these empirical rates 
must be down-revised. If the clumping factor is only a factor of '^4, mass-loss 
rates would need to be reduced by a corresponding factor of ~2. Then, at least 
for O- and early B-stars, empirical rates would be in full agreement with present 
theoretical rates (still assumed to be uncontaminated by clumping effects), which 
are used in current models of massive star evolution. 

If, on the other hand, the empirical rates needed to be down-revised by very 
large factors, of 10 or more, they would fall significantly below theoretical rates. 
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with some devastating effects on current massive star evolution models. This would 
tie in with the downwards revision of WR mass-loss rates about a decade ago, with 
implications involving that black holes have larger masses than previously appre- 
ciated, and the possibility that a significant fraction of the total mass loss would 
need to be lost in super outbursts of LBVs. We note that for a variety of reasons 
including stellar rotation rates and the origin of WR stars, present day stellar evo- 
lution models could not survive with theoretical rates down-revised by more than 
a factor of two, or, equivalently, empirical rates reduced by more than a factor of 
four. 

Only the future will tell! 
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